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Summary 


Photoelectric observations are described, intended to detect or set an 
upper limit to anomalous absorption in the violet wing of the solar K line 
following flares. Such absorption would be direct evidence of a cloud of 
approaching particles, the presumed source of great geomagnetic storms. 

During 1949-50 numerous photoelectric observations failed to reveal in 
most cases any absorption exceeding 0-5 per cent, corresponding to particles 
with approach velocities between 650 km/sec and 1500 km/sec, for several 
hours after flares, even when followed by magnetic disturbance. A few cases 
of doubtful effects up to one per cent are noted. 

Since the strongest absorption would be expected shortly after expulsion, 
the negative results imply a low abundance of Ca II ions in the cloud. This 
has been suggested by theoretical work, though in conflict with earlier photo- 
graphic observations of absorption during magnetic storms. 





1. Introduction.—It is now well established that most (but not all) great 
geomagnetic storms are significantly associated with antecedent solar flares of high 
intensity (1, 2, 3). Newton has compiled lists of great flares and the associated 
magnetic storms, going back to Carrington’s celebrated naked-eye observation in 
1859 (4). It appears inter alia that the mean flare-storm interval varies from 22 
for very great storms to 34" for small storms, and that flares followed by storms 
are usually (but with notable exceptions) within 45° of the visible disk centre 
(r/R<o-7). 

The usual explanation of the association is that these magnetic storms are caused 
by clouds of solar particles emitted simultaneously with the flare within a limited 
conical angle. The mean Sun—Earth velocities corresponding to the above 
flare—storm intervals are from 1900 km/sec to 1200 km/sec. Chapman and others 
have discussed the geometry of electrically neutral but ionized streams of corpuscles 
emitted from the Sun, and the mechanism whereby such corpuscles may cause 
magnetic storms (5). In particular, Chapman has suggested that approaching 
solar corpuscles might be detected spectroscopically, from weak Doppler-shifted 
absorptions shortward of strong Fraunhofer lines (6, 7). 

Some evidence for such displaced absorption during great magnetic storms has 
been found by Richardson (8, g) and by Briick and Rutllant (10). They observed 
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possible absorption bands about one per cent deep in the violet wings of the H and 
K lines of ionized calcium; the corresponding approach velocities were between 
600 and 1400km/sec, and between 700 and 800 km/sec, in the two cases. Such 
small absorption is very close to the limit of detection by photographic photometry, 
and these observers stressed the need for further confirmation. 

Kahn (11), ina theoretical study of the expansion of a corpuscular cloud emitted 
from a model solar flare, has shown that the H and K lines are the most likely to 
show the anomalous absorption. Even for these the absorption after 24" or so is 
not expected to reach one per cent except in the case of a rather dense cloud 
(initial density ~10! Ca atoms per cm‘), and in this case the absorption should be 
very much stronger, and easily observed, immediately after the expulsion of the 
cloud. There is also the possibility of displaced H« absorption during the early 
stages. This is a probable explanation of the frequently observed asymmetry 
of the flare H« profile (12, 13) ; the depression of the violet wing can sometimes be 
traced out to a distance corresponding to an approach velocity of 400 km/sec, and 
may represent an early stage of the cloud’s acceleration. 
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Fic. 1.—Direct-intensity tracing of violet wing of K line. 


(Reproduced from the Utrecht Photometric Atlas of the Solar Spectrum by permission 
of Professor Minnaert.) 


The crowded nature of the solar spectrum in the neighbourhood of H and K 
makes the detection of a small displaced absorption here very difficult. This is 
emphasized by Fig. 1, a tracing from the direct-intensity Utrecht Photometric 
Atlas of the Solar Spectrum (14). At Cambridge the possibility of detecting 
absorption in the violet wings of the D lines, by photographic photometry of 
integrated sky-light, was first explored. However, although this region is less 
crowded with lines, the great variability of the atmospheric water-vapour lines, 
which were not resolved with low dispersion, masked any possible effect due to 
approaching corpuscles. Next the infra-red subordinate triplet of Ca 11 (AA 8498, 
8542, 8662) was investigated. Here the easily identified continuum with few 
absorption lines simplified photometry; but it soon became clear that accidental 
errors due to the rather coarse-grained infra-red sensitive plates would make the 
detection of absorption as small as one per cent quite uncertain in routine 
observations. Attention was therefore turned to photoelectric methods, combining 
high accuracy and simplicity of reduction with the possibility of integrating over a 
broad spectral region. 

In this and a following paper there are described two series of photoelectric 
observations, carried out at Cambridge and Edinburgh respectively, intended to 
detect or set an upper limit to anomalous absorption in the violet wing of the K 
line. Particular attention is given to measurements immediately following flares ; 
it is important to observe a possible corpuscular cloud before lateral expansion has 
diminished its absorption (7). 
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2. Method of observation.—The photomultiplier tube offers great advantages 
in direct photoelectric photometry of the solar spectrum. An amplifier is usually 
unnecessary, the output is linear over a very extensive range, and the signal/noise 
ratio is high. The rather poor 931—A tube available for these investigations 
performed satisfactorily, with adequate light. When requirements are more 
critical, selection from among a number of 931—A’s, which show wide individual 
variations in sensitivity and noisiness, will usually produce a tube with 
characteristics equivalent to the 1P21 (15). 

Direct-recording photoelectric photometers for the solar spectrum, using a 
moving analysing slit and photocell in the focal plane of a spectrograph, have been 
constructed at Cambridge (16) and Kodaikanal (17). Unlike a photographic 
plate the photocell “sees” different spectral regions at different instants, and its 
intrinsic accuracy is not attained unless the atmospheric transparency is exception- 
ally steady during the run. Experience has shown that in the 4000 A region the 
intensity is rarely constant within a few per cent. Consequently some form of 
compensation is necessary; for example, a second fixed photocell may be used to 
check variations of transparency. 

When this is done, the increased precision of photoelectric intensity measure- 
ments must be matched by increased precision in setting on the wave-length under 
examination. ‘This is especially important in a chopped-up region such as that of 
H and K (Fig. 1), where a small wave-length shift would produce a large intensity 
variation. Since a relatively broad shallow absorption was sought, it was con- 
venient in this investigation to use a set of fixed exit-slits in the plane of the 
spectrum, light from which was successively admitted to the photomultiplier. 
Interpolation was then required to allow for the wave-length variation of atmo- 
spheric extinction at different solar altitudes. Thus the necessity for estimating 
the position of the continuous spectrum, an important source of error in such a 
region, was avoided. Further, by integrating the intensity over several angstroms 
and using a wide primary slit, wave-length centering was made less critical without 
sacrificing sensitivity for detecting a broad shallow absorption. The principal 
new source of error was slow changes in the relative effective transmissions of the 
fixed exit-slits arising from accumulation of dust or from small shifts of the photo- 
multiplier relative to the slits. 

3. Apparatus.—A schematic diagram of the apparatus used at Cambridge 
during 1949-50 is given in Fig. 2(a)._ A 12-inch photovisual object-glass, fed by 
an 18-inch heliostat, formed a 55 mm diameter solar image on the slit S, which was 
prevented from overheating by a heat-filter F of Chance ON2o0 glass. The image 
was slightly defocused to avoid local seeing fluctuations. A 132mm x 88mm 
aluminized plane grating G, together with a 9-inch f/9 aluminized mirror M used 
in autocollimation, formed a first-order spectrum at Q with a dispersion of 
8-46A/mm. Fixed exit-slits at Q passed light from selected spectral regions to the 
931-A photomultiplier via a field-lens which imaged the collimating mirror on to 
the photocathode. Preliminary wave-length setting was made by rotating the 
grating, and a line-shifter L provided fine adjustment. 

The exit-slits at Q are shown in Fig. 2(d). A pair of centering slits was used for 
wave-length adjustment, the line-shifter being operated until these slits gave equal 
photocurrents. The slits V, X and R were used for the actual spectrum measure- 
ments, as described below. These five slits were produced photographically on a 
single glass plate, removing the risk of changes of relative position but not of those 
of relative transmission. ‘They were opened and closed by separate light metal 
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shutters. It will be noted that the part of the spectrum examined represented a 
region slightly off-centre on the solar disk; the coordinates of this region, but not 
its angular distance from the disk centre, changed as the heliostat field rotated. 

Fig. 2(c) shows the photomultiplier circuit. Dry batteries totalling 720 V fed 
the dynodes from tappings on a resistance chain, the gain being controlled by 
varying the current with series resistors as indicated. ‘The galvanometer was a 
Tinsley taut-suspension instrument giving a full-scale deflection of 50 cm at I metre 
with a current of about 5A. The small but detectable dark-current was always 
allowed for. ‘The surprisingly high voltage required indicated the low sensitivity 
of the particular 931—A used. 
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Fic. 2.—Cambridge spectrophotometer. 


In addition, a ?-inch diameter selenium barrier-layer cell was later mounted 
redward of the main exit-slits. This was used for compensation for seeing, and its 
output was measured by a Pye taut-suspension galvanometer indicating on the 
same scale as the main galvanometer. 
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4. Reduction of observations.—The approximate wave-length regions covered 
by the exit-slits were as follows : 


R: A3984-°6 to A3994°4 

X: A3915°0 to A39g25-0 

V: A3838-0 to A3861-3 

A3927:0 to A3931°8 

A3935°6 to A3940°3 

The primary slit-width was 0-4 mm, corresponding to 3-4 A. The X slit covered 
the region in which anomalous absorption in the violet wing of the K line might be 
expected. A fairly wide region was chosen, corresponding to a range of approach 
velocities between 650km/sec and 1500 km/sec, at the risk of dilution of the 
absorption if this occurred within a narrower wave-length range. The V and R 
slits, well clear of H and K, were used for comparison, as described later; V was 
made wider than X, and R was covered with a weak neutral filter, in order approxi- 
mately to equalize the photocurrents from each slit. 

For each “observation”’, the K line was first set accurately between the 
centering slits, by adjusting the line-shifter so that when the sliding shutters 
covering these slits were opened alternately, equal galvanometer deflections were 
obtained. Wave-length adjustment was necessary for every “‘observation”’, 
on account of temperature changes in the spectrograph. Several readings on each 
of the three main slits were then taken, and these were averaged, so as to allow as 
far as possible for variations in atmospheric transparency. ‘Thus each “ obser- 
vation ’’ yielded mean relative intensities for each of the spectral regions V, X, R. 

Compensation for the variation of atmospheric extinction (“slope of the 
spectrum ’’) with solar altitude was made by assuming a A~* law, which has been 
shown by Fowle (18) to be valid in the region A 3600-A 4500. In any case, the 
exact form of the extinction law is unimportant over the short wave-length range 
here involved. 
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Fic. 3.—Variation of measured line-prefile with solar altitude. 
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Denote by AV, AX, AR the changes in observed intensity at wave-lengths A,, 
A2, Ag between two observations, as shown in Fig. 3. Let 
J, = relative intensity in solar spectrum at wave-length A, 
S, = relative sensitivity of spectrophotometer at wave-length A, 
x, (t)=total atmospheric extinction (in magnitudes) at time of observation t. 
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Then the photomultiplier current 7, is given by 
log, =log (I, S;) — 2 x, (t) + const., 
and if we consider two observations at times ¢, ¢’, 
A logi, = log, —logi’, = 2 x, (t')— 3 x; (t). 
Now the atmospheric extinction x, (t) is of the form 
x, ()=r-*f(0), 
where f(t) depends on the solar altitude, etc. Thus 
Alogi;=24-* {f(t')—f(0)}. 

So for the two sets of intensity measurements V, X, R, 

A log (X/R)=3 (As *—As) {f(¢')—f(0)}, 

A log (V/R)=3 (y-*—As-*) {f(t’)—F (0)}, 


whence 


Alog(X/R)_Ay#—A*# 
Nlog(V/R) A, *#—A, 4°97 





using the wave-lengths given earlier. Thus any departure of the observed value of 
A log (X/R) from that given by this formula would indicate an anomalous effect. 
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FIG. 4. 
(a) Variation of slope of spectrum during course of day. 
(b) Typical reduction curve. 
(V, X, R represent measured intensities in three fixed spectral regions.) 


The measured intensities V, X, R were reduced to a magnitude scale, using 
specially computed tables incorporating non-linearity calibrations of the galvano- 
meters, and thus the quantities $log(X/R) and $log(V/R) were formed. 
Fig. 4 (a) shows a typical example of the variation of $ log(V/R) during the course 
of a day. At Cambridge the Sun was usually bluer in the morning than in the 
afterncon; consequently the variation of $log(V/R) was asymmetrical about 
noon, and no correction involving sec z only could have been used. 

Now according to the above theory a plot of $ log (X/R) against $ log (V/R) 
should have been a straight line with slope 0-572. In fact, a slight curvature was 
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found when points referring to low solar altitudes were included, as in Fig. 4 (6). 
This effect was presumably due to slight wave-length shifts in the regions measured, 
caused by the method of K-line centering used ; at low altitudes, the violet wing of 
the K line was depressed relative to the red wing, necessitating a slight shift of the 
centering slits to provide balance. ‘The usual procedure therefore was to draw an 
empirical curve, departing as little as possible from the predicted straight line, to 
represent the “normal” variation of § log (X/R) with $ log (V/R), and this 
curve was used so long as the apparatus remained undisturbed. Finally, residuals 
in § log (X/R) were calculated, and these, apart from observational errors, were 
taken to represent any variations in the region X in the violet wing of the K line. 
Atmospheric lines should not cause trouble here; Rowland (19) lists none short- 
ward of A 5292. 

5. Compensation for seeing.—Observations at Cambridge commenced in 1949 
June, and it soon appeared that the chief source of error was the variation in atmo- 
spheric transparency between readings on different slits. This variation, rarely 
less than one per cent and often of the order of ten per cent or more in the presence 
of very slight cloud or haze, could not always be compensated by averaging a 
number of readings. Accordingly, in 1949 September, a selenium cell and galvano- 
meter were arranged to indicate variations of intensity in a fixed spectral region 
several hundred angstroms redward of the V, X, R exit-slits. 

Denote by S the output of the selenium cell, measured by its separate galvano- 
meter indicating on the same scale as the main galvanometer: V and S were read 
simultaneously, then X and S, then R and S. Provided that no considerable 
changes in the “slope” of the spectrum occurred during the few minutes 
of observation, the observed ratios V/S, X/S, R/S were independent of trans- 
parency changes. From these ratios, $log(V/R) and $log(X/R) were found 
immediately. 

Automatic optical or electronic methods of indicating the ratio of two varying 
light-fluxes are possible, and have since been applied in stellar spectrophotometry 
by Hiltner and Code (20). In the present investigation it was simpler and more 
accurate to read the separate photocurrents and to find their ratios; this involved 
differencing pairs of numbers, since a logarithmic calibration of the galvanometers 
was used. 

This method of compensation for seeing enabled the number of readings per 
exit-slit to be reduced to four or so, with a final accuracy in $log(V/R) or 
$ log (X/R) of 0™-003 or better. 

6. Observational results: (a) General—The equipment was in use at 
Cambridge from 1949 June to 1950 August. In the reduction of observations, a 
new reduction curve—relating the “ normal” variation of § log (X/R) to that of 
§ log (V/R)—was drawn whenever any mechanical alterations, apart from wave- 
length adjustments, were made to the spectrophotometer. The curve was usually 
drawn to fit the readings for the best (undisturbed) observing day during the 
period. ‘Thus the observations were divided into groups, each with a common 
zero-point. In spite of this, owing to the somewhat makeshift nature of the 
apparatus its long-term stability was not always good; the average readings 
sometimes drifted by one or two per cent over the course of a week or two, 
presumably due to shifts of the photomultiplier relative to the slits, or to gradual 
unequal obscuration of the slits by dust. This drift should not impair the 
efficiency for detecting fairly sudden onsets of absorption following flares. 
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The observations have been fully listed elsewhere (21). Now Newton (2) has 
shown that while about 30 per cent of flares 3 + are followed by magnetic storms, 
this figure drops to 8 per cent for flares 3 and is scarcely above chance expectation 
for flares 2. For present purposes, therefore, attention is primarily directed to 
observations following flares 3+ and3. These are listed in Table I, together with 
associated geomagnetic effects and any absorption effects observed. The flare 
data have been taken from the Quarterly Bulletin on Solar Activity and geo- 
magnetic data from H. W. Newton’s summaries (22) and from the Fournal of 
Geophysical Research. 

‘TABLE I 
Flares 3 and 3+, 1949 Fune—1950 August 
Position Delayed 
on disk geomagnetic Absorption 
Date Class (r/R) Interval effect effect 

1949 Aug. 5%08!-; 3 0°82 (W) SC Aug. 64 118-2 Uptor per cent, 
64—9" after flare 
No effect with- 
in 1" of flare 
<I per cent on 
Sept. 20, 21 


1949 Sept. 134 138-2 or (W) Period of subnor- 


mal geomagnetic { 


1949 Sept. 174 175-3 -28 (W) activity 


1949 Sept. 184 o9!-8 ‘95 (E) 


1949 Oct. 24 14h - ‘87 (E) SC Oct. 4402"-1 <0 per cent, 
just after flare 
and on Oct. 3, 4 


1949 Nov. 194 1o0!-5 -94 (W) No obs. 


1950 Feb. 174 01-6 ‘66 (E) SC Feb. 194 23"-7_ No obs. 


Great storm 
(H 420 y) 
1950 Feb. 214 238-7 ‘62 (W) 34h- SC Feb. 23 10"-7 No obs. 
Small storm 
(H 155 y) 
1950 Apr. 124 15)-1 ‘43 (E) sit 
1950 May 3%09"-8 ; -31 (CM) 


No obs. 
_ oak No obs. 
( 59)-3 SC Aug. 74 102-g No obs. 
Moderate storm 
1950 Aug. 442356 - ‘25 (CM) 4 (H 240 y) 
44-6 Abrupt movement 
L Aug. 64 20-2 


After the maximum in 1947, solar activity was still relatively high in 1949, but 
declined rapidly during 1950. Great magnetic storms followed two of the seven 
flares 3 + but none of the flares 3 during 1949; the one flare 3 + and two of the four 
flares 3 during 1950 were also followed by storms. ‘The decline in flare occurrence 
after 1950 August almost simulated solar minimum conditions and accompanied 
the onset of a well-marked ‘‘ M-sequence”’. Ofthe eleven flares 3 + or 3 during the 
period of these investigations, listed in Table I, five were near the Sun’s limb 
(r/R >0o-7), and produced no magnetic storms; three of the remainder preceded 
moderate or small storms, but none of these was observed at Cambridge. Thus no 
really favourable opportunity arose for detecting absorption by approaching 
corpuscles. However, we might still expect absorption, at least in the early 
stages, from corpuscular clouds which did not reach the Earth. At present little 
is known about the relative importance of flare intensity and disk position in 
determining the intensity of the subsequent magnetic disturbance. 
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(6) Flares 3, 3 +.—The graph (Fig. 5 a) of observations at the time of the flare 
3+ of 1949 August 5 shows a general rising trend, of instrumental origin; allowing 
for this, we see that there was no absorption exceeding 0-6 per cent (0™-006) 
between 1"-1 and 4" after the flare peak, but an apparent absorption, I per cent deep 
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Fic. 5.—Intensity variations in the violet wing of the K line following flares 3,3+-. Absorption 
(in magnitudes, o™-01=0-9 per cent) is indicated by decreasing ordinate. Single and double small 
arrows denote flares 1, 2 respectively. 


and gradually diminishing, between 6" and 9! after the peak. ‘There is thus some 
evidence of a corpuscular cloud emitted in the region of the flare (r/R = 0-82) and 
covering the Sun’s disk within 6". The only delayed geomagnetic effect was a 
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sudden commencement (SC) after27"-1. This is typical of several cases in Tables 
I (flares 3, 3+) and II (flares 2) which seem to indicate that a flare unfavourably 
placed or not intense enough to produce a magnetic storm may nevertheless 
produce an SC within 48"; no definite conclusion can be drawn from these limited 
data. The west-limb flare of 1949 August 5 in fact occurred during a small 
geomagnetic storm (SC Aug. 3402"-1) which may have been associated with the 
central meridian passage (CMP) of the active region on August 1, but with no 
outstanding observed flare. 

The only readings obtained on 1949 September 13 (Fig. 5 5), within 0"-5 of the 
west-limb flare 3, show no evidence of absorption. Similarly, scattered readings 
on September 20, 21 show no absorption exceeding I per cent following the flare 3 
of September 17 and the flare 3 + of September18. None of these flares produced 
the slightest geomagnetic activity (apart from simultaneous crochets), although 
that of September 17 was close to the disk centre (r/R = 0-28). 

More extensive readings were obtained for 1949 October 1-5 (Fig. 5c). The 
east-limb flare of October 2 was placed by Newton and Jackson (3) in class 3 
because of simultaneous crochet, radio noise and ionospheric effects. Readings 
for October 2 show an apparent drop of about 1 per cent before this flare, but no 
more than 0-6 per cent at the time of the flare and no effect after 15-5; there is no 
evidence of absorption exceeding 0-5 per cent on the two following days. Again 
the only delayed geomagnetic effect of this excentral flare was an SC after 36%. 
The cause of the low readings preceding the flare is not obvious, but it seems 
difficult to associate them with emitted corpuscles. 

Owing to poor observing conditions during 1950, readings were not obtained 


for any of the remaining flares listed in Table I. It is remarkable that two of the 
three flares near the disk centre produced no geomagnetic effect. 


TABLE II 
Flares 2 during observing periods 


Position Delayed 
Date Class __ on disk Interval geomagnetic 
(r/R) effect 
1949 June 2% 151-8 0-42 (W) 300-1 SC June 34212-g <2 per cent, 
Moderate storm 2-3 after flare 
(H 230y) 


1949 June 22417"-3 24 ‘84 (E) <0°5 percenton 
239 145-9 ‘77 (E) June 22, 23; up 
2770296 2 -30 (W) oe ~ to 1 per cent, 12 
279 16-1 72 (W) after flare on 
289 o9h-5 “50 (W) June 27 and 65 
after flare on 
June 28 


Absorption 
effect 


1949 July 3406-5 -38 (E) is sis <0°5 per cent, 
7» after flare 


1949 Sept. 24 10!-7 
34 08-9 
44 1ol-s 
54 o2h. 
5a o8h.- 
5d 12h. 
69 ogh. 


‘99 (W) o°5 per cent, 
‘70 (E) within 1" of 
‘53 (E) flares on Sept. 3, 
‘44 (E) SC Sept. 84 10"+1 5,6, orbefore SC 
"31 (E) 
‘85 (E) 
‘75 (E) 


N NNN NN WN 
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TABLE II (cont.) 


Position Delayed 
Date on disk Interval geomagnetic 
(r/R) effect 


1949 Sept. 124 06)-5 76 (W) 
124 73)-3 ‘80 (W) on +. <I per cent at 
124 y5h-5 76 (W) 142 40m 
124 2oh-2 ‘82 (W) 


Absorption 
effect 


1949 Sept. 214 065-8 ‘60 (E) _ an <I per cent, 
2h_oh after flare 


1949 Sept. 264 112-2 ‘69 (W) es mS <0°’5 per cent, 
just after flare 


1949 Oct. 34 03!-3 75 (E) <0°5 per cent, 
34 11)-9 75 (E) just after flare on 
44 082-9 ‘62 (E) SC Oct.6%07"-3 Oct. 3; up to 
44 130-4 -58 (E) Small storm 0°5 per cent, 2" 
49 16-4 “50 (E) (H 140y) after flare on Oct. 
(442096  ? 2 43< 0°5 percent, 
54 1-3 -51 (E) 3" after flare on 
Oct. 5 


114 92!-2 1-2 “49 (E) 

114 97h. ‘81 (W) 

rr yyzh- “88 (E) SC Oct. 13420"-2 <1 per cent on 
rd yyh. ‘47 (E) (325-8) | Moderate storm Oct. 12 

rr zsh. 2 ‘41 (E) (H 220y) 

114 ,6h- “56 (E) 
(124 yr hg -28 (E)) 


1949 Oct. 244 og!-2 -78 (W) SC Oct. 27404)-8 <1 per cent 
Small storm during storm 
(H 1507) 


1949 Oct. 294 15)- . SC Nov. 1409-9 <1 per cent on 
Small storm Oct. 31 and 
(H 160y) during storm 


1950 Mar. 84 og!- +58 (W) ies ae Up to 1°5 per 
cent,3"after flare 


1950 May 1041-2 24 ‘45 (E) SC May 11% 175-2 < 05 per cent 
rh_sh after flare 


1950 May 23% 08!-8 ‘66 (E) a ne <0°5 per cent, 
just after flare 


7950 June 19 06-6 “11 (CM) <0°5 per cent, 
14 7gh.3 -11 (CM) ? SC June 67004-3 22 after flare on 
34 07-4 ‘34 (W) June 1; 
59 o7h-7 77 (W) < per cent, 2)- 
45 after flare on 
June 3 or 1—4h 
after flare on 
June 5 
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(c) Flares 2.—In order to make full use of the observational data it is worth 
looking also at suspected geomagnetic effects of flares 2. While Newton (2) has 
shown that on average there is no significant rise of activity following flares 2, 
there are isolated cases of a probably genuine association with subsequent storms. 
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Fic. 6.—IJntensity variations in the violet wing of the K line following flares 2, 1949. Absorp- 
tion (in magnitudes, o™-01=0-9 per cent) is indicated by decreasing ordinate. Single and doubl: 
small arrows denote flares 1, 2 respectively. 

Accordingly, in Table II there are listed all recorded flares 2 during observing 
periods with the present apparatus, together with suspected geomagnetic and 
absorption effects. The classification of flares is at present rather arbitrary, and 
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undoubtedly many flares of minor importance are included here. It is, however, 
possible in some instances to trace the connection of magnetic disturbances with 
specific flares 2, or with the CMP of flare-active regions. 

The fairly central flare 2 of 1949 June 2 preceded by 30"-1 a moderate magnetic 
storm with SC. The readings before the flare (Fig. 6a) are fairly steady: the 
single reading 2"-3 after the flare is uncertain, but definitely limits absorption to 
less than 2 per cent. 

None of the numerous smaller flares between 1949 June 22—28 produced any 
geomagnetic disturbance. ‘The continuous series of observations obtained 
(Fig. 65) fails to show any absorption exceeding 0-5 per cent, except for possible 
drops of the order of 1 per cent about 1 after a flare 1-2 on June 27 and 6" after a 
flare 1-2 on June 28. A single reading on July 3, 7" after a central flare 1-2, 
shows no absorption greater than 0-5 per cent. 

The SC of 1949 September 8 is an example of a geomagnetic disturbance 
preceded by the CMP ofa flare-active region but not obviously associated with one 
particular flare. The readings at this time (Fig. 6 c) show a falling trend, presum- 
ably of instrumental origin; apart from this there is no evidence of absorption 
exceeding 0-5 per cent within 1" of a flare 2 on September 3 or within 0-5 of flares 
2 on September 5, 6. 

A single reading on 1949 September 12 shows no evidence of absorption within 
1"5 ofaflare2. The readings for September 21, already discussed in connection 
with a flare 3 + (Fig. 5 5), indicate also the absence of absorption greater than 1 per 
cent between 2-9? after a flare 1-2. Similarly on September 26 (Fig. 6d) there 
was no absorption exceeding 0-5 per cent during the hour following a flare 2. 
Readings on that day appeared, however, to recover at about 10" from a depression 
of rather more than I percent. It would be tempting to associate this effect with a 
central flare 1 (r/R=0-28) at 06%, although corpuscular emission from a flare 1 
seems unlikely. ‘There was no geomagnetic effect except an SC recorded at one 
station only, 24-7 after this flare. 

The flares listed between 1949 October 3—5 occurred in an active region which 
crossed the central meridian on October 7. This region had already produced on 
October 2, while near the east limb, a flare 3 followed by an SC (see Table I). A 
small storm commenced with a further SC on October 6; we might associate this 
with the observed flare 2 + on October 4 at 13"-4, or perhaps with an (unobserved) 
flare at 20"-6 deduced from a geomagnetic crochet and a sudden ionospheric 
disturbance, giving an interval of 345-7. Readings already discussed (Fig. 5c) 
show a slight drop of notmorethano:5 per cent commencing about 2" after the flare 
2+ on October 4, but no absorption exceeding 0-5 per cent during 3" following the 
flare 2 on October 3; nor was any effect found 3"-5 after the flare 1-2 on October 5. 

Although there were several fairly central (r/R=o0-5) flares 2 on 1949 October 
II, it seems more probable that the SC and moderate magnetic storm of October 
13 were associated with the central (r/R = 0-28) flare of October 12, 32"-8 previously. 
Classified as only 1+ at Greenwich, this flare was accompanied by marked iono- 
spheric disturbance, and its importance may have been underestimated. A few 
readings before the flare on October 12 show no evidence of absorption. 

A flare 2—, accompanied by crochet, on 1949 October 24, is probably the 
only flare which we can associate with the SC and small magnetic storm of October 
27. ‘There is no evidence of absorption exceeding I per cent during the storm. 
The SC and small storm of November 1, if not of the recurrent type, might be 
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associated with the east-limb flare of October 29, variously classified as 1, 2 and 3. 
There is little evidence of absorption greater than 1 per cent on October 31 or 
during the storm. 
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Fic. 7.—Intensity variations in the violet wing of the K line, 1950. Absorption (in 
magnitudes, oO™-01=0°9 per cent) is indicated by decreasing ordinate. Single and double small 
arrows denote flares 1, 2 respectively. 


Even flares 2 in 1950 were infrequent. That of 1950 March 8 produced no 
magnetic disturbance; the few readings (Fig. 7a) show a falling trend of little 
more than I per cent up to 3" after the flare 1-2, or 1 after the flare 1, which pro- 
duced a stronger ionospheric effect. An SC occurred 30" after the flare 2+ of 
1950 May 10, but there was no absorption exceeding 0-5 per cent or so between 
15 after the flare or on the following day (Fig. 75). A pair of readings just 
before and 5™ after the peak of the flare 1-2 on 1953 May 23 show no absorption of 
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as much as 0°5 per cent; there was no geomagnetic disturbance. Finally, an SC 
followed by slight disturbance on 1950 June 6 may have been associated with an 
active region which crossed the central meridian on June I. Scattered readings 
(Fig. 7c) show no evidence of absorption greater than 0-5 per cent within 2" of a 
central flare 2 on June I, or greater than I per cent between 2"—4» after a flare 1-2 on 
June 3 or between 1"—4 after a flare 1-2 + on June 5; low readings on June 3 are 
accidental, due to poor observing conditions. 
Tas_e III 


Magnetic Preceding Absorption 
Date disturbance Interval solar event effect 


1950 Mar. 19% o5!-7 SC 32h-7 crochet <2 per cent on 
storm (H 245) (flare-active 18th 
region near CM) 


1950 Apr. 224 og!-2 SC vee i <I per cent on 
20th; 

234 058-8 SC <o'’5 per cent 
on 21st; possible 
effect 1-5 per 
cent on 22nd 


1950 May 2% 09-8 SC 40!l-0 crochet <2 per cent on 
storm ~13!" (central flare 1?) Mayr 
(H 1757) 
1950 Aug. 184 155-6 SC ?25h-7 flare 1 <2 per cent on 
great storm (r/R=0'17); 16th; <1 per 
(H 310y) crochet cent, 7" before 
storm 

(d) Other magnetic disturbances.—We have seen that geomagnetic distur- 
bances may in some cases follow less intense flares. It is noteworthy that the 
second greatest storm during the period (H 385y, SC 1949 Oct. 154 08-1) followed 
445-4 after a flare 2+ on the central meridian. Also, the flare responsible for a 
storm may occasionally be unobserved. Thus a few storms associated with flares 
may be omitted from Tables I, II; in Table III, therefore, there are listed geo- 
magnetic disturbances, during observing periods, for which no intense preceding 
flare was observed. 

The moderate magnetic storm of 1950 March 19 occurred two days after the 
CMP of a flare-active region; a geomagnetic crochet, such as often accompanies 
more intense flares (23), preceded the storm SC by 32"-7. A few readings on 
March 18 definitely limit absorption then to less than 2 per cent. 

Doubtful SC’s on 1950 April 22, 23 were preceded by no outstanding solar 
event. Readings for April 20-22 (Fig. 7d) show no evidence of absorption 
exceeding I per cent on April 20 or 0-5 per cent on April 21; the drop of 1-5 per 
cent on April 22 is not obviously explained. 

A crochet, probably accompanying a central flare 1, occurred 40" before the 
SC of 1950 May 2; asmall magnetic storm developed 4" later. A few readings on 
May 1 limit absorption to less than 2 per cent. 

The great magnetic storm of 1950 August 18 followed 25"-7 after a crochet 
accompanying a central flare 1 (r/R=0-17). Readings on August 16 show no 
absorption exceeding 2 per cent; there is no evidence of absorption on August 18, 
7» before the SC. 
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7. Conclusion.—The observations described have failed in most cases to 
reveal any absorption of as much as one-half per cent, either immediately after 
flares or during magnetic storms, due to approaching Ca tlions. A few instances 
of apparent absorption, for example following the flare 3 + of 1949 August 5, may 
represent a genuine effect. In view of possible observational errors, however, 
these should be regarded with caution. 

The negative results do not fulfil expectations aroused by earlier isolated 
photographic observations, but confirm Kahn’s prediction (11) of the difficulty of 
detecting a corpuscular cloud spectroscopically. Using Kahn’s tables, we see 
that absorption of less than one per cent within a quarter of an hour of the flare 
implies an initial density for the cloud of not more than 10! particles per cm%, 
assuming Stroémgren’s value of 6 x 105 for the ratio H: Ca. A greater proportion 
of Ca would cause heavy absorption unless the cloud were so much less dense as to 
prevent appreciable recombination of Cam to Cam. The upper limit of 101 
particles per cm? is considerably lower than Svestka’s estimate (24) of 10! for the 
particle density in an intense flare, but evidently there is uncertainty in both figures. 
Again, on Kahn’s model the upper limit to the density of this cloud on reaching the 
Earth would be 400 particles per cm’, which is still considerably greater than the 
usual estimates of the density required to produce a magnetic storm (5). 

In relation to the earlier photographic observations of Richardson (8, g) and of 
Briick and Rutllant (10), however, the present results indicate lower absorption 
than might have been expected. The one per cent absorption found by them 
during magnetic storms would imply much stronger absorption soon after the flare, 
which the present observations have failed to reveal. The possibility of dilution 
of the absorption, owing to the wide spectral band used, cannot be overlooked, but 
it could scarcely reduce strong absorption to an undetectable value. Similarly, 
the range of approach velocities covered (650 km/sec to 1500km/sec) could not 
exclude a major part of the absorption in every case. And although the spectrum 
of only one part of the solar disk was examined, we should expect a corpuscular 
stream to cover the disk fairly soon after emission; the comparative frequency of 
storms following flares shows that the stream cannot usually be limited within a 
conical angle of much less than 45° (1). 

Thus, although no case of a central flare 3+ followed by a great storm was 
observed, the present results throw some doubt on the reality of the effect found 
by the previous observers. ‘Those observations were advanced with caution; 
internal evidence, such as the dissimilarity between the absorption profiles in the 
wings of H and K, indicated considerable uncertainty. 

It remains to be investigated whether the small effects occasionally found here, 
of the order of one per cent, are genuine. It may be necessary to look for local 
absorption on the Sun’s disk, if the crucial density of Ca 11 for detection is reached 
only for a short time during the corpuscular cloud’s expansion. ‘There may also 
be a stage at which there is detectable displaced H« absorption, distinct from the 
asymmetry of the flare profile itself. 
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ELECTRON TEMPERATURES AND ELECTRON DENSITIES 
IN PLANETARY NEBULAE 


M. F. Seaton 
(Received 1953 September 29) 


Summary 


Electron temperatures T, and electron densities N, are estimated from 
a study of the relative intensities of the forbidden lines in the following 
nebulae : NGC 7027, 2440, 7662, 6572, IC 418 and Orion. 

For certain ions (N 11, O 111, O11 and S11) it is possible to measure the 
relative intensities of two multiplets originating from different metastable terms. 
Each intensity ratio gives a relation between T, and N,. A further relation 
may be obtained if it is assumed that N,(O 1)/N,(N )=N,(O 11)/N,(N 01), 
where N,(X) is the number of ions X in the ground state per cm*. Recently 
calculated electron excitation cross-sections enable these various relations to 
be obtained in numerical form. 

The observed relative intensities of the forbidden lines are corrected for 
absorption and systematic calibration errors by comparing observed and 
theoretical relative intensities in the hydrogen spectrum. Simultaneous 
solutions for T,, and N, are then obtained from the forbidden line intensity 
ratios. The values of 7, range from 1:3 x10‘ to 2-0 x 10 deg. K, and the 
values of N, from 8 x 10° to 5 x 10 per cm*®. The consistency of the various 
T., N, relations is particularly satisfactory for NGC 7027, NGC 2440 and 
IC 418; for the other objects considered the observed intensities tend to be 
either less reliable or less complete. 

With the exception of IC 418, the densities obtained from the forbidden 
lines are greater than those obtained from the surface brightness by a factor 
of order 3. These differences are probably a consequence of the densities 
being non-uniform. ‘The surface brightness method gives essentially a 
r.m.s. density for the total geometrical volume, while the forbidden line 
methods give average densities for the active emitting regions (clouds and 
filaments). Such differences would not be expected for IC 418, which is 
particularly uniform. 

In a final section the abundance ratios N(O)/N(N) are estimated to be 1-6, 
0°43 and 2:0 for NGC 7027, NGC 2440 and IC 418 respectively. It is 
concluded that the results reveal a real difference in composition between 
NGC 2440 and the other two objects. 





’ 1. Introduction 


In order to make a trustworthy determination of the chemical compositions 
of planetary nebulae, and to solve a number of other outstanding problems, 
it is essential to have reliable estimates of the electron densities N, and the 
electron temperatures T,,. 

The earliest electron temperature estimates were obtained from studies of 
the relative intensities in the hydrogen recombination spectrum (1, 2). Since 
the relative intensities are insensitive to 7,, this method is not capable of giving 
much better than order-of-magnitude estimates. 

All previous density estimates (3, 4) have been obtained from estimates of 
the absolute emission rate in the hydrogen recombination spectrum. This 
method gives N, to be of order 104 cm-* in typical planetaries. To apply this 
method it is necessary to determine the absolute surface brightness and the 
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absolute nebular dimensions, both of which are difficult to measure to a high 
precision. 

In principle it should be possible to obtain more reliable temperature and 
density estimates from the relative intensities of the forbidden lines. For a 
given ion the relative intensities of the forbidden lines are determined by the rates 
of electron excitation, electron deactivation and spontaneous emission. If the 
electron excitation and deactivation cross-sections and the radiative transition 
probabilities are known, the relative intensities may be expressed in terms of 
known functions of 7, and N,. 

The ions giving the forbidden emissions contain two metastable terms. When 
deactivation is negligible for both of these, which requires that the transition 
probabilities should be sufficiently high and the electron density sufficiently 
low, the ratio of the emission rates from the two levels will be proportional to the 
ratio of the excitation rates. The latter ratio may be expressed as a function of 
T, alone; hence if the intensity ratio is known, 7, may be calculated. 

To a good approximation deactivation is unimportant for Oil. For most 
other ions it must be taken into account for the lower metastable level but not 
for the higher one. For these ions the intensity ratio depends on both 7, and 
N,, and hence may be used to determine N, if T, is known. 

In practice it is desirable to consider the relative intensities for several ions 
and to solve simultaneously for N, and 7,; if a unique solution for N, and T, 
is obtained we can have some confidence in the correctness of the results. 

Temperature estimates were first obtained from the [OIII] intensity ratios 
by Menzel, Aller and Hebb (5), who used the cross-sections calculated by Hebb 
and Menzel (6) and the intensity ratios measured by Aller (7). The tempera- 
tures obtained were of order 10‘deg. K and varied comparatively little from 
one nebula to another. 

Some doubt was cast on the correctness of these results when it was shown by 
Bates, Fundaminsky, Leech and Massey (8) that the Hebb and Menzel cross- 
sections were much too large, since they violated a general theorem on charge 
conservation. A different method of calculating the cross-sections has recently 
been developed by the author (9, 10, 11) and cross-sections have been obtained 
for a number of transitions of importance in planetary nebulae. The results 
should be sufficiently reliable. 

In the present paper the new cross-sections are used to calculate temperatures 
and densities in a number of planetary nebulae and also in the Orion nebula, 
for which the same methods may be used. 


2. The hydrogen spectrum 


It is convenient to summarize certain formulae* used in discussing the 
hydrogen spectrum (Menzel (12), corrected for a numerical error (4)). The 
density of atoms in the mth quantum state of hydrogen (ionization potential x,,) 
is written in the form 


apa he 72, ae 
N,(H)=by | pe | nPexp(xalT INN (2) 





* Certain formulae are included which will not be required in the present paper but which will 
be required in future papers which are in course of preparation. It was thought that a compact 
summary would be the most convenient. 


I1* 
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which defines 6,(6,=1 for thermodynamic equilibrium). The coefficients 5, 
have been calculated as functions of T, by Baker and Menzel (13). WN, is 
the proton density and the other symbols have their usual significance. 

The energy in ergs radiated per cm® per sec in transitions from state n to 
state n’ is given by 





I(n, n')=1-42 x 10-18 foe, =P (xn/®T.) N,N. (2) 


(nx')> * T 3? 

where g,,,, is the Kramers-Gaunt factor (13, 14). The intensities of a few 
Balmer lines, relative to H8 = 100, are quoted in Table I from Baker and Menzel 
(13) (Case B). 

TABLE I 


Theoretical Balmer intensities 
I 


An. 


A , 
(A) T,=1 X 10 T,=2 xX 108 T.=4 10 
6563 250 259 271 
4861 100 100 ‘100 
4340 51 50 49 
4102 31 30 29 





The intensity per angstrom in the continuum resulting from recombination 
to quantum state m is given by 


1a8xn I exp[— (hv—-x,, WRT ed) yy Nz 


I(x, n) =6°47 x 1077 3 Ts (3) 





where A is the wave-length in A and « is a quantum number for continuous states 
defined by 


hy= x, (= + “): (4) 


K 


The intensity per angstrom resulting from free—free transitions in the fields of 
protons is given by 
h 
1 f—f)=2-05 x 10° J ame pa NN, (5) 


where (Menzel and Pekeris (14)) 
a — hv\18 2kT, 6 
g=1+0172 - ade 9 (6) 


3. Theory of the forbidden line spectrum 


3.1. Notation.—Most of the forbidden lines observed in planetaries have 
been identified with transitions between the three terms contained in the ground 
configurations of various ions. These configurations are of the form 2p? or 3p% 
with g=2, 3 or 4. The three terms in these configurations will be denoted by 
n=1, 2 and 3 in order of increasing excitation energy. The relation between 
this notation and the usual spectroscopic notation is given in Table II. 


TasLe II 
Terms in 2p% and 3p% configurations 
q n 


3 
3 
3 
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The energy differences between the J-levels of the terms are much smaller 
than the mean kinetic energy of the electrons and may therefore be neglected. 
We further assume that there is a Boltzmann distribution among the J-levels 
(i.e. within a given term the populations of the J-levels are proportional to 
(2/+1)). This will be further discussed in Section 3.5. 

The statistical weight of term n is denoted by w,, and the term energy by E,, 
(centre of gravity). For £,,>E,, it is convenient to introduce 


Enm = (E,—E,,). (7) 
The radiative transition probabilities A,,,, (sec~!) are averaged over initial 
J-levels and summed over final J-levels, i.e. 


I , 
stl (2L,, + I)(2S,, + ha. (2J,, + 1I)A(L,, Syn» J, Sia Lin» SinvJ m)- (8) 


The number per cm* of ions of a given type in state nm is denoted by N,. If it is 
necessary to specify the type of ion this is written N,(X). The energy radiated 
per cm® per sec in the transition n—> m is given by 
I(n, m)= EnmA nmN n- (9) 
The electron excitation and deactivation cross-sections Q(m — n’) are con- 
veniently expressed in terms of the dimensionless parameters 


A 





Q(n, n') = (7) 2s O(n—n’), (10) 
7 
introduced by Hebb and Menzel (6). v,, is the velocity of the electron incident 
on state # and the other symbols have their usual significance. Seaton (10) has 
suggested that the Q(n, n’), which are symmetrical in m and n’, should be termed 
the collision strengths. 
The probability per sec of electron impact producing the transition n — n’ 
is denoted by 
Inn’ =2nwN,- (11) 
Assuming a Maxwell distribution of electron velocities (temperature T, deg. K) 
the deactivation coefficients «,,, are related to the excitation coefficients «,,, by 
the relation 


Cyn = —" Oy EXP(—Eqn/RT,), (En >Ey)> (12) 
which ensures detailed balancing for collisional processes. The deactivation 
coefficients for positive ions are given by 


Q2(m, n) 


Onm = 8°54 x 10-6 — Tis? (13) 


it being assumed that the collision strengths are independent of electron energy 
(10). For neutral atoms (2(m, n) in (13) must be replaced by 
~ 4s 
Y(T,;n,m)= | Q(m, n)exp(—me,3/2kT,)d( = ) (14) 
Jo 2k1,, 

3.2. The atomic data.—Table III gives the transition probabilities, collision 
strengths and energies for the ions with which we will be concerned in the present 
paper. 

For Nu, Our and O1 the transition probabilities of Garstang (15) have 
been used. For Ni, Orr and Si we have used the A,, of Garstang (16) and 
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the A, and A3, of Pasternack (17) after correcting for the improved s, integrals 
of Garstang (16). 

The collision strengths are quoted from Seaton (10,11). They are believed 
to be correct to within + 40 per cent in absolute value ; the relative values such as 
Q(1, 2)/Q(1, 3) are probably of a higher accuracy than this. 

Table IV gives the deactivation coefficients for the neutral atoms O1 and NI 
(Seaton (9, 11)). The values for NI are estimated and are therefore of a 
slightly lower accuracy, but it is considered improbable that they are in error 
by a factor as great as 2. 

The accuracy of the atomic data for S11 will be somewhat lower than that 
for the other ions since non-exchange ion wave functions have been used. 


TABLE III 


Atomic data 
ats eure MO). 0G, 9. . 00,39 
0°034 1°08 2°39 0°223 0°46 
0°23 16 1°73 O°195 0-61 
0°0054 0°079 : ~ ae 
0°046 390135 1°44 0218 1°92 
0078 1°28 . ibe be 
0°26 0°44 2°02 0°383 12°7 


TABLE IV 


Deactivation coefficients for O01 and N 1 
(&,m in cm® sec!) 
O1 NI 


aie —<—_? oo 


10°) 10°%3) 10° X39 10°Qo; 10°a3; 








°'9 o'8 o'4 o'5 o'4 
1°6 1°2 o'6 08 06 I°5 
3°5 1°9 26 18 4°5 
5‘! 3°! 3°5 2°7 72 
72 6-1 4°0 3°7 14°3 
6-6 7° 3°5 3°4 16°4 
3-3. The equations of equilibrium.—The equations of equilibrium have been 
discussed by Menzel, Aller and Hebb (5) and by Barbier(18). They are discussed 
from first principles in the present paper in order to examine the validity of the 
various approximations which are made. A simplified general discussion has 
been given elsewhere (19). 
For planetary nebulae we may assume 7", to be of order 104 deg. K and N, 
to be of order 10'cm™*. It is convenient to introduce 
t,=10*T, 
and 
x= 10-*N, /t,", (16) 
both of which will be of order unity. 
Equating the number of ions entering each state to the number leaving we 
obtain* 
No Agi + Ga1 + 923) = Ni9i2 + Na( Ase + 9s2)s (17) 
N3{ Agi + Age + 931 + 952) = N91 + Nogas- (18) 





* It is readily shown that re-absorption is unimportant; cf. Menzel and Aller (20). 
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We consider first the derivation of a simplified expression for N,/N3. 
Eliminating N, from (17) and (18) gives 


Ns — i fem +931) + (Ase + Gsa)lI + 913) al 
Nz 41s \(Ae1 + Gox) + (912928/913)[ 1 + 913/912! 





(19) 
Our first approximation is to assume that 
his <i. (20) 
fiz 
Equations (12) and (13) show that 
dis _ Q(1, 3) 
fie S2(1, 2) 


and since 02(1, 3) <{Q(1, 2) in all cases of interest (Table II), (20) is satisfied for 
all values of 7,. Our second approximation is to neglect deactivation of state 3, 
which is justified if 


exp (— E5,/kT,), (21) 


Ag> 4s, and As2> 452. (22) 
Using (13) it is readily shown that this requires that 
* €10*A5,,093/8-5.4Q(m, 3) (23) 


for m=1and m=2. ‘Table V gives the quantity on the right-hand side of (23) 
for a number of typical ions. It is seen that (22) is satisfied for x of order unity. 


TABLE V 
Ion 10A 3,03/8-54Q(1, 3) 10*A 5,w5/8°542X(2, 3) 
Nu 180 2 700 
Ol I 400 3 100 
N1* 200 I 100 
Ou I 500 49° 
O1* I 500 41 000 
Su 4 800 240 
* For neutral atoms the {2(m, 3) are replaced by the Y(T,; 3, m) (equation (14)), 
calculated for T= 10%. 


Using these two approximations, (19) simplifies to 


Pet tele 
Ns 413 (Aor + 921(931 + 932)/9s1 


For a number of ions we observe two multiplets, the first arising from the 
2-1 transition and the second either from the 3-1 transition or from the 
3-2 transition. For a given ion X the intensity ratio of these two multiplets 
is denoted by 





(24) 


YX) = 1(2, 1)/1(3,m), (25) 
where m=1 or 2. Using (9), (24) and (25) the intensity ratio may be expressed 
finally as 





r Kon ex it, 
Yon( X) = as <P en ) (26) 





M. F. Seaton, 


Electron temperatures and 





Kow= FOG, 3) 


€39 = 10-*E30/k, 


d, = 8-54 x 10-4 


Ayn 





Q(1, 2) (Gs 3) + Q(2, 2) 
Agi. Q(1, 3) . 
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(27) 
(28) 
(29) 


and where ¢, and x are defined by (15) and (16). 
The quantities A;, As» Kin), dy and €go are given in Table VI for the ions of 
principle interest. 
TABLE VI 
Ion An Asm Kim) d, Ese 
Nu 6548 +6584 5755 9°65 O°312 2°50 
Om 4959+ 5007 4363 8°74 0°0438 3°30 
Ou 3727 7319 +7330 17°4 15°7 1°96 
Su 6717+ 6731 4069+ 4076 8-60 5°51 1°39 


We consider next the ratio N,/N,. Eliminating N, from (17), (18) we obtain 


Ns ow 91s( Age + 932)/(Asi + Ase + 9s1 + 952) 
Ny Agi + er + G23( Asi + 951)/(Asi + Ase + 951 + 92)’ 
and assuming (20) and (22) as before this gives 


Nz - iz 
Ny Agy + Goi +923451/(Agi + Az)’ 
which may be simplified to 





(30) 





(31) 


No Qe 
N, ~ Ag+ 4a (32) 

if 
Yar > 923Ag1/(Agi1 + Ase), 


2(2, 3) As, 
Q(1, 2) | Ag, + Aso 


The quantity (Q(2, 3)/Q(1, 2))[As;/(A3; + A32)] is given in Table VII for the ions 
of interest. It is seen that (33) is satisfied, to a sufficient accuracy, for all these 
except SII. 


(33) 





exp (— E3,/kT,) <1. (34) 


TABLE VII 


Q(2, 3) Ax 
Q(1, 2) | Asi t+ Ase 


0°005 9 
0-045 
0079 
0°34 
0005 3 
2°34 

* See note to Table V. 


Ion 





Niu 
Or 
N 1* 
O1 
O 1* 
Siu 


Using (32) the final expression for the intensity may be written as 


1240 


I(2, 1) o En Ay + OV, 


N,N,. 
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3.4. Use of the [O1] and [N1] intensities.—The following lines, all belonging 
to 2—>I transitions, are observed in a number of nebulae* : [0 1], A6300 + 6364 ; 
[N1], A5199; [Om], A3727; and [Nu], A6548+6584. Using equation (35) we 
obtain 1(6300 + 6364) _ Ge,(O1) N,(O1) (36) 

(5199) = Gy (N1) N,(N1) 
- 1(3727) _ _ Gu(Ou) (On) 
(0548 +6584) ~ Ga(Nut) N,(N)’ (37 
where Gy; = Ey1%12A9;/(Ao1+%,N,). A further relation between N, and T, 
may be obtained if it is assumed that 
N\(O1) N,O1) 
a 8 
N,(N1i) ~~ Ny(Nm)’ (38) 
A comparison of the ionization potentials (13-55 eV for O1 and 14-46 for N1), 
and of the photo-ionization cross-sections (Bates and Seaton (21)) suggests that 
the relation (38) should be satisfied approximately. We use (38) to eliminate 
the abundance ratios from (36) and (37). In substituting numerical values for 
the atomic constants we neglect deactivation for state 2 of O1 (which requires 
N,<1-4 x 108 for T,=10*); take «,(O 1)/a»,(N 1) = 48/25 (11) and put 
&9,(N 1) = 3°6 x 10-® cm® sec™ 
(cf. Table IV). We obtain finally 
(6300 + 6364) 1(6548 +6584) _ (1+ 3°4t,'x)(1 + 1°6x) , 
15199) 1(3727) = 0140 “Tr rorosy _ *P (t/t). (39) 

3.5. Relative populations of the J-levels.—For the ions with g=3 (Ou, N1, 
S11) and for the nebulae with which we are concerned, observations of the 
relative intensities of the lines in the various multiplets (?P,—?D,,, ?P,;—>4S3, 
and ?D,->4S,,.) do not show any evidence for systematic departures from the 
Boltzmann distributions among the J-levels of the 7D and ?P termsf (Aller, 
Ufford and Van Vleck (22), Garstang (16)). 

For the ions with g=2 and 4 we may consider the condition for a Boltzmann 
distribution to be established among the J-levels of the ground *P terms. The 
equations of equilibrium are similar to (17), (18) and it is readily shown (6, 
18) that these give a Boltzmann distribution if 

g(*P,—>*P,)> APP, *Py) (40) 
for all J, J’ such that the energy of *P, is greater than that of *P,. Using the 
collision strengths of Seaton (11) and the transition probabilities of Garstang 
(15), it may be shown that the relations (40) will be satisfied for O11 if «>0-6. 
In some nebulae x is of order 0-6 (cf. Section 4.2 below); we may therefore 
expect some departures from the Boltzmann distributions in these nebulae. 

















*(O1] A5577 is also observed in a number of cases, but the measured intensity is probably 
largely due to nightglow radiation. Since (6300+ 6364) and (5577) are comparable in the 
nightglow, the values of [(6300+ 6364) measured for nebulae should be comparatively free from 
nightglow radiation if 1(6300-+6364) > 1(5577). 

So far as the author is aware the [N 1] line at A 3467 has not been reported in nebular spectra. 
It would be of interest to look for this and also, in further work on the far infra-red, to look for the 
[N 1] lines at AA 10407, 10398. 

+ Note added when correcting proofs (1954 July 24).—For nebulae of lower electron density there 
are appreciable departures from the Boltzmann distributions among the J-levels of the *D terms of 
Ou and Su; the relative intensities of the lines in the 7D —‘S multiplets may then be used to 
obtain estimates of electron density (M. J. Seaton, Ann. d’Ap., 17, 74, 1954). 
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The equilibria between the terms will not depend on whether there is a 
Boltzmann distribution among the J-levels if the collision strengths Q(°P,, 'D,) 
and Q(8P,, 1S,) are proportional to (2J +1). It is probable that this condition 
is approximately satisfied (10). It appears safe to conclude that small departures 
from the Boltzmann distributions will not introduce serious uncertainties into 
the equilibria for the terms. 


4. Determination of T, and N, from relative intensities of forbidden lines 


4.1 The observational data.—Relations between N, and T, may be obtained 
from the relative intensities 7,,)(X) for the following ions: Ou, Om, Nu and 
Su. A further relation is obtained from (39) if the [O 1] and [N 1] intensities 
are also known. 

The amount of accurate photometric data available is not as complete as 
could be desired, and it is frequently necessary to combine results obtained by 
different observers. ‘The most complete lists available are those of Wyse (23), 
covering the wave-length range 3705 to 6755 A and including many fainter lines. 
Although these intensities are visual estimates, comparison with photometric 
intensities indicates that they are in general correct to within a factor of 2. The 
Wyse estimates will be used only when photometric intensities are not available. 

In comparing the intensities of lines with large wave-length differences it is 
necessary to make a correction for differential space absorption (together with 
errors arising from atmospheric absorption and instrumental calibration). This 
correction may be made using the theory intensities for the hydrogen lines. 
Although space absorption depends statistically on the path length through the 
stratum of obscuring matter in the galaxy, large individual deviations occur due 
to the uneven distribution of the obscuring clouds. It has been shown by 
Berman (24)* that the average relative hydrogen intensities, corrected for space 
absorption as a function of absorbing path length, are in good agreement with the 
theory values. We may therefore assume that the relative emission intensities 
are given correctly by the theory. The space absorption may then be estimated in 
individual cases by comparing the theoretical and observed hydrogen intensities. 
The dependence of the theory intensities on temperature is so small that tempera- 
ture uncertainties do not introduce any appreciable uncertainty into this 
correction. In all cases we use the theory intensities of Baker and Menzel (13) 
(Case B) calculated for T,=2 x 104deg. K. 

We consider the correction in some detail for NGC 7027, a nebula for which 
the observational data are particularly rich and the space absorption particularly 
large. We denote the theory intensities by J, and the observed intensities by 
Io. The Ip of Aller (7) have been used for the Balmer lines, with the exception 
of H« for which the value of Wyse has been used. H{ (A388) has been omitted 
since it is blended with Het. For the Paschen lines the J, of Aller and Minkowski 
(26) have been used. 

Fig. I gives a plot of log(J,/I,) against the wave number » in reciprocal 
microns (v=10*/A with A in A). The corrected intensity for any line of wave- 
length A is given by 

1(A) = fil oA), (41) 


where log(f;) is the function defined by the smooth curve drawn in Fig. I. 


* See also Baker and Menzel (13) and Aller (25). 
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The adopted forbidden line intensities, relative to HB = 100, are given for a 
number of nebulae* in Table VIII. Smoothed absorption corrections have been 
applied using the method described above for NGC 7027. A slightly different 
correction has been applied for [O 111] 44363. In a number of cases Hy (A 4340) 
appears to be anomalously weak compared with the other hydrogen lines. This is 
probably a consequence of calibration errors, and it is therefore likely that 
1(4363) isin error by acorresponding amount. The procedure used has therefore 
been to take 1(4363) relative to Hy =50 (cf. Table I). 

The following detailed comments may be made: 

NGC 7027. The observational data are particularly complete and there is 
fair general agreement between the results obtained by various observers (Aller 
(7), Wyse (23), Aller and Minkowski (26, 31), White (27), Page (1), Plaskett 
(29) and Berman (28)). The [O11] intensities will be discussed in detail in 
Section 4.4 below. 


+1-Or— 


-OS5 





3.0 20 





_—> 


Fic. 1.—Absorption correction for NGC 7027, obtained from a comparison of theoretical and 
observed hydrogen intensities. I is the theory intensity, Ig the observed intensity and v the wave 
number in reciprocal microns (v=10*/A with A in A). 


NGC 2440 and NGC 7662. Space absorption appears to be small and no 
corrections have been made. For NGC 7662 Wyse states that (6548 + 6584) 
is of a low accuracy. It has therefore not been used. 

NGC6572. This nebula has been studied by several observers, but 
unfortunately there exists a number of discrepancies in the results obtained. 
The Balmer intensities indicate appreciable absorption for A<4861 (H§), but 
the Ha/H§ ratio given by Wyse does not indicate absorption; no absorption 
corrections have been made for A>4861. The value of /(6717 +6731) given by 
Swings and Jose (30) (relative to Hx) is much greater than the value given by 
Wyse. Swings and Jose give the following Paschen intensities relative to 


* These are the only nebulae for which the available observational data are sufficient to allow at 
least two N,, T, relations to be obtained from the forbidden line intensities. 
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Ha=50: P,,=o—1, Py=1, Pys=1. According to theory we would expect 
P,;, Py, and P,, to be of order 0-15 (Ha=50). The Paschen intensities have 
been disregarded. For /(7319 +7330) the intensities of Swings and Jose have 
been used relative to Hx=250. The value of 1(6548+6584) given by White 
(17) (relative to Hx) is much smaller than the value given by Wyse. 

IC 418. The intensity of Hy given by Wyse appears to be particularly low. 
Judging by the other hydrogen intensities we would expect the observed intensity 
of Hy to be about 41. Wyse obtains 25. 

Orion. The intensities of Ha, H8, Hy and Hé given by Wyse indicate appre- 
ciable absorption, but the intensities of the higher Balmer lines are in fair agree- 
ment with theory without making an absorption correction. In consequence 
absorption corrections have been made for all intensities in Table VIII except 
for that of [O 1] A3727. 

The author is indebted to Drs L. H. Aller and R. Minkowski for communi- 
cating to him in advance of publication certain [O 111] intensities (31) obtained at 
Mount Wilson during the years 1945-47. 


TaBLe VIII 


Adopted relative intensities of forbidden lines 
(relative to HB=100 and corrected for absorption) 





NGC 7027} NGC 2440 INGC 7662 INGC6572| IC 418 Orion 

















A Transition 

Int. |Ref.| Int. | Ref. | Int. |Ref.| Int. |Ref.| Int. | Ref.| Int. |Ref. 
3727 Om 21] 45 7 80/7,25| --- |... | 100 | 7|261| 7 | 80] 23 
4069 4076|Sil 31]10°;| 7 8] 23 | 2°, | 23] 36] 7/3°0°] 7] 5 | 23 
43603 Or 32 32] 31 28 7 22] 33 | 12°; | 713°°O} 23 | 4's | 23 
4959 5007 | OIl! 21] 1480] 31 | 1540| 23 | 2250] 313 | 1260] 7 |151]| 31 | 532] 23 
5199 Mi 2212.48 30| 23 <n crs 2| 23 | 2°, | 23 | o'8] 23 
5755 Nil 32 7 7 50] 23 _ sa 12] 23/9 23 | 3°3 | 23 
6300 6364;O1 21] 33 7 | 100} 23 ee ee 23 | 23 | 8-6 | 23 | 4:2 | 23 | 
6548 6584 | Nir 21] 141 |23, 7| 1040 |23, 27) ... | --. | 150] 23 | 220| 23 | 82 | 23 | 
6717 6731 |S 21] 22] 26 12| 23 | 2:9 | 23 5| 23 | 14] 23] 10 | 23 | 
7319 7330|O1m 32] 41] 26] ... ae ost hk oo, 500) OO iid x4 















































The references in Table VIII are: Aller (7), Wyse (23), Aller (25), Aller and 
Minkowski (26, 31), White (27), Swings and Jose (30). When two references are 
given this indicates that both observers obtain practically identical results. The 
intensities of Wyse have occasionally been used to estimate the relative intensities 
within blended features observed by other workers. 


4.2. The solutions for T, and N,.—The relations between 7, and N,, obtained 
from equations (25) and (39) and the intensities of Table VIII, are shown graphic- 
ally in Fig. 2. The curves obtained from (39) are marked ““O/N”’. We also 
include the curves which are obtained from the surface brightness (Menzel and 
Aller (3)). These are marked “ S.B.”’. 

In the surface brightness method the absolute emission rate in the hydrogen 
spectrum is calculated from estimates of the surface brightness and nebular 
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dimensions.* The 7,, N, relations are then obtained using the theory of the 
hydrogen spectrum summarized in Section 2. In previous work the surface 
brightness at the head of the Balmer continuum has usually been used. This 
has the advantage that no knowledge of the calculated 5,, coefficients is required. 
It has the disadvantages that comparatively large errors arise from space absorption 
and that it is necessary to make a rather uncertain correction for the underlying 
“visual continuum” (Menzel and Aller (3), Page (1, 2)). It is considered that 
any uncertainties in the calculation of the b,, coefficients will be smaller than these 
observational uncertainties. We have therefore used the surface brightness of 
Hf, as estimated by Menzel and Aller (3), and have used Baker and Menzel’s 
Case B for the 4,, coefficients. 

It is seen that the 7, N, relations obtained from the forbidden line intensities 
give satisfactorily consistent results for NGC7027, NGC2440 and IC 418. 
In view of the various intensity uncertainties which have already been mentioned, 
it is scarcely surprising that the consistency is less satisfactory for NGC 6572. 
For this object the relations obtained suggest that /(6548 + 6584) may have been 
under-estimated ; the effect of increasing /(6548 + 6584) would be to lower the 
N 11 curve and to move the O/N curve to the right. It also appears possible that 
1(6548 + 6584) may have been under-estimated for the Orion nebula. 

The densities obtained from the forbidden lines will be denoted by N, (F.L.) 
and those obtained from the surface brightness will be denoted by N, (S.B.). 
In Table LX we give what appear to be the most probable values for 7',, N, (F.L.), 


N, (S.B.); we also give x= 10-4N, (F.L.)/1/(10~‘T,). 


TABLE IX 

Nebula wT. 10-*N, (F.L.) 10-*N, (S.B.) x 
NGC 7027 1'9 4 1*4 2°9 
NGC 2440 I-97 4 sen 3°1 
NGC 7662 1°3 5 14 4°4 
NGC 6572 1°3 5 1°4 4°4 
IC 418 20 o'8 o-7 o°6 
Orion 1-3 2 o-1-1°0* 1°8 


* Estimated by Strémgren (32). 


4.3. Discussion of N,.—It is seen that, with the single exception of IC 418, 
N,(F.L.) is considerably greater than V,(S.B.). Despite the various obser- 
vational uncertainties there would appear to be no doubt that real differences 
occur.{ Such differences could result from the densities being non-uniform ; 
the two methods would then give different results due to different methods 
of averaging. 


* The nebular dimensions are calculated from estimates of the angular diameters and of the 
distances of the nebulae. 

+ The relation of the intensities of the Balmer continua relative to the Balmer lines will be 
discussed in a later paper. It appears that Case A may be a better approximation than Case B 
for certain objects. If Case A were used the values of N, obtained from the surface brightness 
should be multiplied by a factor of 1°31. 

} Note added when correcting proofs (1954 July 24).—Further work, in course of publication in 
Monthly Notices, has shown that neglect of space absorption may cause N, (S.B.) to be considerably 
under-estimated. It is found that the surface brightness estimates obtained from photographic 
magnitudes by Menzel and Aller (3), when corrected for space absorption, give values of N, (S.B.) 
which do not differ systematically from the adopted values of N, (F.L.). However, recent photo- 
electric measurements of surface brightness by W. Liller and L. H. Aller (Astr. 7., 59, 186, 1954) 
suggest that the estimates obtained from the photographic magnitudes may be too large. 
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It is readily shown that the quantity determined by the surface brightness 
method is 


N,(S.B.)= EI Ne ar |, (42) 


where N, is the local density and V is the total nebular volume assumed (sphere 
or spherical shell). In this expression we have neglected possible temperature 
variations and assumed that V,=N,. Thus the surface brightness method gives 
a root mean square density for the total nebular volume. 

The forbidden line methods give average densities for the regions where the 
forbidden line emissions are the most intense. The explicit expressions for these 
averages are in general rather complicated. We consider only the simplest case, 
the determination of the density from the [O 1] intensity ratio, assuming T, to 
be known and to be constant throughout the nebula. Ifin addition we assume that 


9a(O 11) > Ag,(O 11) (43) 


threughout the regions giving appreciable emission (which requires x> 0-6), 
it is readily shown that the quantity determined is 


- | N.N,(On) dr 
~ [NOn)ar | 


The differences between N, (S.B.) and N,(F.L.) may be discussed further 
by considering a very simple model. Wesuppose all the matter to be concentrated 
into clouds which occupy a fraction A of the volume V. We also suppose that 
the density, the ionic composition and the temperature are constant within the 
clouds and the same for all clouds. With this model N, (F.L.) would be equal 
to the density within the clouds and (42) would give 

N, (S.B.) = AN, (F.L.). (45) 
Using this model the densities of Table IX would give A~o-1 for NGC 7027, 
6572 and 7662. 

Photographs show that local condensations (clouds and filaments) occur 
in many planetary nebulae (see Curtis (33)). The subject has recently been 
discussed by Zanstra (34), who has suggested a mechanism responsible for the 
formation of condensations. Ina number of cases it appears to be very probable 
that the conditions are such that this mechanism will be operative (19). 

It is probably reasonable to assume that the fraction of the total volume 
occupied by dense clouds is of order one-tenth in certain nebulae, of which 
NGC 7662 is a good example (see the sketch by Curtis reproduced by Wurm 
(35)). It is of interest that IC 418, for which N,(F.L.) is not significantly 
greater than N, (S.B.), is particularly uniform. According to Wilson and Aller 
(36) no fine structure is visible in even the best photographs of this nebula. 

For uniform objects such as IC 418 it would be possible to obtain distance 
estimates by the following procedure: 

(i) Calculate T, and N, (F.L.) from the relative forbidden line intensities. 

(ii) Measure the surface brightness and from the usual theory of the S.B. 
method deduce the nebular dimensions, assuming N, (F.L.)=N, (S.B.). 

(iii) From the measured angular diameter calculate the distance. 





N{On) (44) 


It is of interest to note that the only fundamental length units involved would 
be those required for the measurement of the fundamental constants used in the 
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atomic calculations and those required for the measurement of the absolute surface 
brightness. 

4.4. Discussion of T,.—For NGC 7027, NGC 2440 and IC 418 there can be 
little doubt that the temperatures obtained in Table IX are essentially correct.* 
These temperatures are appreciably greater than the results obtained earlier (5) 
(T,=1 x 10) using the [O 11] ratio alone. So far as [O 11] is concerned, the 
higher temperatures obtained in the present paper result from: 


(i) The use of the improved collision strengths ; 

(ii) The method of correcting for absorption and calibration errors (the 
method of correcting 1(4363) relative to Hy is particularly important for IC 418) ; 
and 

(iii) The use of more recent intensity determinations. 

Further [O 111] intensity measurements have recently been made by Andrillat 
(37> 38), who obtains results which are often markedly different from those 
obtained by all previous workers. In Table X we give the values of T, calculated 
from Andrillat’s [O 111] intensities.t 


TABLE X 
Temperatures obtained from [O 111] intensities of Andrillat 
Nebula oT, Nebula 10“T, 
NGC 7662 4°0 NGC 7009 2'0 
NGC 7027 3°0 NGC 6572 1°9 
IC 5217 2°5 NGC 6543 I'l 
NGC 6803 3°0 Orion rs 


These results can scarcely be accepted without reservation. In the case of 
NGC 7027 it is difficult to reconcile a temperature of 3 x 10* with the other N,, 
T, relations given in Fig. 2. It is also difficult to reconcile the results of 
Andrillat with the results obtained by previous workers. In Table XI we give 
T, for NGC 7027 as obtained from all available measurements of the [O 111] 
intensity ratio. In this table we have made no corrections to 1(4959 +5007) 
relative to H8=100, but have taken J(4363) relative to Hy=50. 


TABLE XI 
T, for NGC 7027 calculated from the [O 111] intensity ratio 
Observer 10“T, 
Andrillat (37) 3°0 
Plaskett (29) 19 
Aller and Minkowski (31) 1°9 
Berman (28, 24) 1°75 
Wyse (23) 16; 
Aller (7) 1°6 


These results suggest that the earlier intensities of Aller (7) give 7, somewhat 
too small and that the intensities of Andrillat (37) give it considerably too large. 


* For IC 418 the temperature of 2-0 x 10‘deg. K is in good agreement with a result obtained 
independently by Dr R. Minkowski (private communication) from a study of the Balmer dis- 
continuity. This work has been published (Pub. Astr. Soc. Pac., 65, 161, 1953) since the present paper 
was written. 

+ Andrillat’s published results give only the temperatures, calculated with the older collision 
strengths. The author is indebted to Mr H. Andrillat for supplying him with the detailed 
intensities. 
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4.5. Further discussion of T,; theoretical considerations.—It is frequently 
assumed that 7’, is constant throughout a given nebula, or at least that the same 
effective value of T,, can be used in the discussion of different processes. ‘This 
assumption requires further examination. 

The factors determining 7, in planetary nebulae have been discussed by 
Menzel and Aller (20) and a similar problem for interstellar clouds has been 
discussed by Spitzer (39, 40) and by Spitzer and Savedoff (41). From these 
discussions it may be shown that, for a given distribution of energy in the ultra- 
violet radiation field of the star, the electron temperature is determined by the 
ratio of the sum of the intensities of the forbidden lines to the intensity of any 
convenient hydrogen line, say Hf. We denote this ratio by /(F))H3. 7, will 
be small when J(F)/Hf is large. 

Variations of 7’, will occur in a given nebula if /(F')/HP varies. ‘T'wo types of 
variation may be considered : 


(i) Large-scale variations, such as could result from changes in ionization 
equilibrium with changes in the distance from the central star, and 

(ii) Local variations due to changes of ionization equilibrium occurring in 
condensations. ‘This has been discussed by Zanstra (34). 


If temperature variations occur it is possible that different effective tempera- 
tures should be used for different processes, since different methods of averaging 
may be involved. Ina number of planetaries by far the greatest contribution to 
I(F) comes from the [O 111] lines. It follows that, if temperature variation occurs 
in these objects, in determining 7’, from the [O 111] ratio we may give the greatest 
weight to those regions where /(/’)/Hf is greatest and where the local tempera- 
ture is least. Such effects could result in the temperatures obtained from the 


[O 111] ratio being less than the effective temperatures which should be used for 
other processes. 

The results obtained in the present paper do not show any convincing evidence 
that such effects are important for the nebulae considered, but the possibility of 
such effects being important in other cases should be borne in mind in future work. 


5. The O/N abundance ratios 
For NGC 7027, NGC 2440 and IC 418 the “O/N” curves of Fig. 2 are 
seen to be consistent with the curves obtained from the other forbidden line 
intensity ratios. ‘This confirms that equation (38) is satisfied to a good approxi- 
mation. It would appear reasonable to assume that the ionic abundance ratios 
in (38) do not differ much from the total chemical abundance ratios N(O)/N(N). 
To calculate N(O)/N(N) for these three objects we proceed as follows. We 
first adopt values of 7, and N, such that (38) is satisfied for the calculated ionic 
abundances; these values of 7’, and N, do not differ significantly from the most 
probable values previously adopted. We then calculate the ionic abundance 
ratios using (36) or (37) and assume that these are equal to N(O)/NM(N). The 
results are given in Table XII. 
Tas_e XII 
Nebula -*T 10*N, N(O)/N(N) 
NGC 7027 1'9 4°9 1°6 
NGC 2440 1°65 4°5 2°43 
IC 418 20 0°75 20 
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The theoretical interpretation of the forbidden line intensities gives satis- 
factorily consistent results for the three nebulae considered and we can therefore 
have some confidence in the reliability of the derived O/N abundance ratios. 
In particular there appears to be little doubt that N(O)/N(N) is a good deal smaller 
for NGC 2440 than for the other two objects, indicating a real difference in 
chemical composition. It should be noted that we have selected objects in which 
the [N1] and [N11] intensities appear to be reliably determined and hence in 
which nitrogen may be particularly abundant. Other planetary nebulae may well 
have larger values for the N(O)/N(N) ratio. 

It is of interest that in many respects NGC 7027 and NGC 2440 are very 
similar: both are bright irregular objects in which no central star has been 
observed, they have similar values for T,, and N, and, apart from the differences 
in the intensities of the oxygen and nitrogen lines, they have very similar spectra. 


TABLE XIII 
N(O)/N(N) 
Mean of eight B stars (42) 4°4 
7 Sco (42) 2°6 
10 Lac (42) 44 
Sun, Menzel (43) 3°0 
Sun, Unséld (43) 13 
Sun, Hunaerts (43) 16 


For comparison we give in Table XIII some values of N(O)/N(N) obtained 
for the Sun and stars. Although there is a considerable scatter among the 
results obtained, it is doubtful if the methods used in these cases are sufficiently 


precise to establish the existence of real differences. 


Department of Physics, 
University College, 
London : 
1953 September 25. 
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METALLIC AMMONIUM 


M. F. M. Bernal and H. S. W. Massey 
(Received 1954 February 3) 


Summary 


It has been pointed out by Ramsey that metallic ammonium (NH,) may be 
important in determining the properties of Uranus and Neptune. The 
pressure—density relation for this metal at pressures of importance in this 
connection is calculated by quantal methods. The transition pressure for 
the phase transition from mixed crystals of ammonia and hydrogen (NH; and 
4 H,) to ammonium metal is also investigated and found to be almost certainly 
less than 250 000 atmospheres, so that the metallic phase will persist through- 
out the bulk of the two planets mentioned. 





Introduction.—The importance of allowing for phase changes which occur at 
high pressures, in determining the constitutions of the planets, has been emphasized 
by Ramsey (1). Unfortunately the investigation of the conditions under which 
these phase changes take place is extremely difficult. The transition pressures are 
usually beyond the reach of available experimental technique, while the processes 


are usually too complex for theoretical treatment. An exception in the latter 
respect is the transition from the molecular crystal to the metallic form of hydrogen. 
This has been treated quite thoroughly theoretically by Wigner and Seitz (2) and 
by Kronig, de Boer and Korringa (3). Their results, which provide, in numerical 
form, an equation of state for hydrogen at pressures of the order of those 
encountered within the major planets, have been applied by Ramsey (4) to the 
elucidation of the constitution of these planets. He finds that whereas Jupiter 
and Saturn are predominantly hydrogen, Uranus and Neptune, which have much 
smaller masses, have lost most of their hydrogen. For these planets the princtpal 
constituents are likely to be water, methane and ammonia. This raises an inter- 
esting opportunity for further theoretical investigation. 

It has long been known that the ammonium ion NH,+ behaves in many ways 
like an alkali metal ion. The existence of the so-called ammonium amalgam, 
produced when sodium amalgam is mixed with a solution of ammonium chloride, 
shows that this similarity persists to some extent into the metallic phase. In other 
words ammonium metal is probably stable above a comparatively low transition 
pressure. Ramsey (4) points out that metallic ammonium is therefore likely to be 
an important constituent of Uranus and Neptune. Owing to the comparatively 
high symmetry of the molecular structure of the NH,* ion it is possible to carry out 
calculations of the properties of metallic ammonium which should be of sufficient 
accuracy for application to planetary constitution. Apart from this application 
such calculations are also of interest as providing an estimate of the order of magni- 
tude of the transition pressure for the phase change from the mixed ammonia and 
hydrogen molecular crystals to the metal. 

In this paper we give the results of a theoretical investigation on these lines. 
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Calculation of the equation of state of metallic ammonium.—The basis of the 
method used is provided by the work of Buckingham, Massey and Tibbs (5) on the 
properties of methane. They obtained a self-consistent field for this molecule by 
first averaging the proton distribution over a sphere. The resulting central 
field problem may then be dealt with by exactly the same methods as for the self- 
consistent fields of atoms. ‘The results obtained were very encouraging and 
suggested that, for many purposes, the approximation is quite a satisfactory one. 
To deal with the ammonium problem the field of the ammonium ion is first obtained 
as for methane by averaging the tetrahedral distribution of protons over a sphere of 
radius ry. The value of rz was chosen from a preliminary investigation (6) carried 
out by one of us (M. J. M. B.) in which the total energy of the tetrahedral ion was 
calculated as a function of ry by a simple variational procedure. The minimum 
value of r, obtained in this way (r,=1-84a,)) was then adopted for the further 
calculations. 

Having chosen rp, a self-consistent field for the ion with spherically symmetrical 
proton distribution was obtained in the usual way. In Table I the self-consistent 
wave functions and the effective nuclear charges (Z,,) for potential are given. 

Once the self-consistent wave functions are obtained the usual procedure of the 
theory of metals (7) may be adopted to calculate the energy of ammonium metal as a 
function of the nitrogen—nitrogen separation. The metal is divided into polyhedra 
centred round each nitrogen nucleus, and the ground-state wave function of the 
metallic electron must contain the correct number of nodes within the polyhedron 
and have vanishing slope at the boundary walls of the polyhedron. ‘To simplify 
the calculation it is usual to represent the polyhedra by spheres of radius a, 2a being 
the separation of the nitrogen nuclei. It has been shown by Wigner (8) that the 
various corrections which should be applied to the results obtained by this 
simplified procedure are not very significant for the alkali metals. Since 
ammonium metal is likely to be essentially similar, we have ignored these 
corrections. 

The equation for the wave function P/r of the metallic electron is, in atomic 
units, 


P"(r)+2[e-—(V+r'Zq)| P(r)-A,, P,,(7) — Ag, Pos(7) 


=-22(al+1)| | Pury) Pni'dn+r' | Para) Prt dr, | 
0 r 


nl 
Pr(r)s (1) 
where V is the potential due to the nitrogen nucleus and proton distribution 


fe - 642, r <1, 
, % 


r2N 


r+ P,, are the normalized NH,* wave functions (the values of P,,, given in Table I 
are not normalized), and Z, is the effective central charge of the NH,* electrons 
Zen=2Z,(1s)+2Z, (2s) +6Z, (2p), 


where Z,,(1s), Z,,(2s) and Z,,(2p) are given in Table I. 





TABLE I 
Self-consistent wave functions and Z, values for NH,* 
P (1s) P (2s) P (2p) Z,(15) Z,(25) Z,(2p) Lot 
©0000 ©0000 ©0000 00000 00000 0°0000 ©0000 
0°1739 0°1737 ©0022 0°1318 00193 0'O159; ©°7957 
0°3028 0°3012 00084 02563 00383 00319, 1°5614 
0°3956 0°3905 0°O177 0°3698 00569 00478, 2°2811 
0°4598 0°4487 0°0295 0°4705 0°0750; 00638, 2°9478 
O°5014 0°4813 00433 05582 00927 ©0797, 375605 
0°5252 ©4930 00587 06334 ©1099; 00956, 4°1219 
0°5353 00-4877 0°'0753 06974 0°1268 o'r116, 4°6361 
0°5348 04686 0°0929 O°7512 0°1433 0°1275) 51083 
05262 0°4384 OvIIII 07962 0°1596 0°14337 5°5436 
05116 0°3993 0'1298 0°8336 0°1756; 0°1592 5°9472 
0°4708 0°3015 0°1679 08899 0°20735  0°1907, 6-6774 
0°4219 0°1869 0:2060 0°9278 0°2387; 0°2219, 7°3296 
0°3709 0:0639 0°2433 0°9530 0°2700, 0°'2528,  7°9263 
0°3214 —0:'0616 0'2791 09695 0°3013;5 0°28334 84837 
02753 —0'1853 0°3130 09804 0°3326 0°31334 90119 
0°2337, —0°3043 = 0344809874 03638; 034277 95 184 
01969 —0-4166 0°3743 0°9920 0°3948;  0°37159 10°0064 
0°1649 —0'5210 0°4014 0'9949 0°4256 0°39974  10°4788 
0°1373.  —0°6168 0°4262 09968 0°4560 0°427Ig 10°9370 
O°1139 +—0°7037 0°4486 0:9980 0°4857,; 0°4538; 11°3809 
00941 —0°7817 0°4688 0°9987 0°5148,  0'47972 11°8111 
0°0638 —o-g119 05031 ©°9995 05706, 0°5290, 12°6296 
0'0429 —I‘OI05 0°5301 0-9998 06226; O'575I1g 13°3914 
00287. —1:0816 0°5508 0'9999 0°6704; 0°6178, 14°0956 
O°OI94 —I°I294 0°5663 10000 0°7139 0°6572, 14°7426 
0°0134 —1°'1580 0°5774 0°7529 0°69351  15°3337 
0°0097  _—I°I709 05849 0°7877 07266, 15°8709 
00078 = —I1'I7II 05893 08184 0°7569,  16°3569 
—1'1611 o'5gII 08454 0°7844 16°7944 
—1°1I79 05879 08893 0°8318 17°5388 
—1°0528 05766 0'9219 0°8700 18-1276 
—0°9725 0°5572 0°9457 09005 = 1185888 
—o'8808 0°5295 09626; 0°9244,  18°9440 
—0°7835 0°4952 ©°97455 09430 19°2141 
—0o'6872 0°4573 09828 0°95723; 19°4180 
—o-5962 04182 09884, 0'9680, 19°5705 
—0°5127 0°3795 09923 09762, 19°6842 
—0°4377 0°3421 09948; 09824 19°7683 
—0'3714 03066 09966, 0°9870, 19°8307 
—0°3135 0°2736 0-9978 0°9904,  19°8765 
—0°2634 02431 0°9985;  0°99299 199101 
—0°2204 O'2152 ©'9990; 0°0048, 19°9346 
—0°1837 0°1894 0°99904, 0°9962, 19°9524 
—0°1§25 0-1671 0'9996 0°99726 19°9655 
—0o-'1261 0°1467 ©°9997; 0°9980, 19°9751 
—o'1038 0°1284 09998; 09985; 19°9820 
—0:0857 O'1120 ©°9999 09989,  19°9870 
—o-0581 = (00848 0°99995  9°99945 199933 
—0°0391 0:0638 I 0000 ©°99972 19°9966 
—0'0262 0°0476 0°9998, 19°9983 
—0°0I74 =: 0°0353 ©°99993 _19°9992 
—O°Orls 00261 ©°9999; 19°9996 
—0°0076 = 00192 ©°9999s + 19°9998 
—0°00§0 = O04 ©9999» 19°9999 
—0'0033 00103 I'0000 20°0000 
—0'0021 0°0075 
N (1s) =12°35788 } " 
N (2s) = 0500120 (N(nl)}* f  {P(nl)}? dr=1 
N (2p)= 1°450171 : 
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The parameters A,, and A,, are adjusted so that r~! P is orthogonal to r~! P,, and 
gt Be 

The terms on the right-hand side represent the exchange interaction. Strictly 
speaking the core functions P,,, should be determined from equations of the same 
form as (1), in which the effect of interaction with the metallic electron is included. 
This effect is likely to be very small and has been ignored throughout. Further- 
more in (1) the contribution from terms involving the Is electron is also very small 
and has been neglected. 

The function P/r is of 3s type. « must therefore be determined so that P has 
two nodes between r = 0 and r=a and d/dr(P/r) vanishes at r=a. 

The numerical calculations were carried out by an iterative procedure which 
converged quite quickly. Numerical solutions were first obtained with the 
exchange terms in (1) ignored. For a particular value of a, the function P, 
obtained in this way was substituted in the right hand-side of (1). The resulting 
equation was solved with A,, and A,, put equal to zero, to obtain a new approximation 
P,, €, to P and e, for this value of a. ‘This approximation does not satisfy the 
orthogonality condition but it is not difficult to guess a modified function P, + 7 P., 
which does satisfy it, at least approximately. This function is then substituted 
for P on the right-hand side of (1) and further approximate solutions for P and « 
obtained. It was found from experience that this approximation is accurate 
enough for the present calculations. 

Fig. 1 illustrates the variation of « with a* obtained in this way. The con- 
siderable importance of the exchange interaction may be seen by comparison of 
curve I, obtained by neglect of exchange, with curve II, in which it is included. 

So far the calculations have referred only to the lowest level in the conduction 
band of the metal. The actual energy per polyhedron relative to the gaseous 
ion is obtained by adding the mean Fermi energy. For alkali-like metals this is 
given quite closely (7) by 


3 for\28 Ilo ‘ lila : , , 
ep= —(— atomic units if a is also measured in atomic units. 


~ 10\4 47?ma a 

This energy is represented as a function of a* by curve III in Fig. 1. The final 
calculated variation of the electronic energy with the volume per nitrogen atom in 
the metal is then represented by curve IV in Fig. 1. This gives the amount per 
nitrogen atom by which the energy of the metal falls below that of the free ion 
NH,*+. The equation of state of the metal at o deg. K may now be obtained from 


the relation 
de 
ae (2) 
p being the pressure when the volume per ion is $7a’*. 

Table II gives the pressure—density relation obtained in this way for a range of 
pressures of interest in connection with the constitution of Uranus and Neptune. 
It remains to consider what can be said about the transition pressure for the phase 
change to the metallic state. 

The transition pressure.—To calculate the transition pressure the energy— 
volume relation for the mixed molecular ammonia—hydrogen crystal and for the 
metal must be available, both of course calculated from the same zero. The slope 
of the common tangent to the two curves then gives the transition pressure a 
o deg. K. 
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To carry out this procedure it is necessary to determine the difference in energy 
between the two phases when each is in its lowest state. We have, per nitrogen 
atom, denoting the energy per molecule of a compound XY by [XY], 

[ NH,, metal] =[NHsg, gas] — Proton affinity of NH; +(e+ eyo, 
[NH,, crystal] + } [H,, crystal] =[NHs, gas] — J_ — 4 D(H.) 
+e(NH,, crystal) + }¢ (Hg, crystal) 
+ €49. (3) 
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‘1G. 1.—Calculated energy—volume relations at o deg. K for metallic ammonium. 
Curve I. Ground-state energy calculated neglecting exchange. 
Curve II. |Ground-state energy calculated with inclusion of exchange. 
Curve III, Mean Fermi energy. 
Curve IV. Calculated total energy. 


The energy zero is that of a free NH4* ion. 

The point O is the estimated equilibrium point for mixed crystals (NH3+4H;,). OP is the 
estimated common tangent to the curves for the metallic and mixed crystal phases. The uncertainty 
in the energy of the mixed crystal relative to the metal due to the uncertainty of the proton affinity 
of ammonia is indicated by the lines AA’, BB’ respectively. 
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I, is the ionization energy of atomic hydrogen and D (H,) the energy of dissociation 
of H,; —«(NHsz, crystal) and —e(H,, crystal) denote the binding energies per 
molecule in the respective crystals. ¢,, is the energy per NH, molecule arising 
from mixing the two molecular crystals. 


TABLE [I 
Pressure—density relation for metallic ammonium at o deg. K 
Pressure (101* dyne/cm?) Density (g cm~*) 
4°0 
3°5 
3 
2°5 
2 
I'5 
axe) 
o'8 
06 
0-4 
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‘To determine the difference between the metal and molecular crystal energies 
we still need to know the proton affinity of NH, and the molecular crystal binding 
energies. ‘The proton affinity of NH, is large and must be known with reasonable 
accuracy for present purposes. It cannot be calculated directly but may be 
determined in principle as follows (g). Let Q be the heat evolved in the reaction 


$ No (gas) + 2H, (gas) + } Cl, (gas) = NH, Cl (solid). (4) 
‘This reaction may be regarded as equivalent to the successive reactions 
} Ny (gas) + $ Hy (gas) = NH, (gas) + Q,, 
} H, (gas) = H (free) + Qy, 
H (free) = H* (free) +e+Q,, 
} Cl, (gas) = Cl (free) + Q,, 
Cl (free) +e= Cl (free) +Q,, 
NH; (gas) + H* (free) = NH,* (free) + P, 


NH, * (free) + Cl- (free) = NH, Cl (solid) + Q,. (5) 
Hence 


QO=Q0,+Q2+03+Q4+ 05+ P+Q,, (6) 


where the Q’s are the heats evolved in each reaction and P is the required proton 
affinity. All of the Q’s in (6) are known from direct measurement except Q,. 
This must be obtained from the heat of sublimation of NH,Cl by the method of 
Born and Mayer (10). 

This procedure can be applied equally well to the other ammonium halides 
and the spread in the resulting values of P gives some indication of the accuracy 
of its determination. In electron-volts we find P= 9°55, 9:11, 8-94 and 8-72 from 
the fluoride, chloride, bromide and iodide respectively. ‘This suggests that 
P=9-06 + 0°5 eV but the uncertainty may be even larger than this. 
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The molecular crystal binding energies are of the order 2 kilocalories per mole 
or o-reV. ‘This is considerably smaller than the uncertainty in the proton 
affinity and can therefore be ignored. We now have 


Proton affinity of NH,= 9-06eV Ionization energy of H =13°54eV 
—(e+ep)y= 5°84eV } Dissociation energy of H,= 2°22 eV 
14:90 eV 15'76eV 


The minimum energy per nitrogen atom of the metal is clearly quite close to 
that of the mixed molecular crystal. ‘The above calculation shows it to be only 
o-86eV larger. On the other hand, the uncertainty in the proton affinity indicated 
above is as much as +0-5eV and other uncertainties arising from the approxi- 
mations made in the calculation are likely to be comparable. ‘This makes a 
reasonably accurate determination of the transition pressure difficult. 

The equilibrium volume in the mixed crystal is not known from direct obser- 
vation, but for present purposes it will suffice to suppose that the NH,—NH, and 
H,~—H, separations remain the same as in the separate molecular crystals. 

The densities of these crystals are 0-82 and 0-087g cm-? respectively. It 
follows that the mean volume per nitrogen atom in a stoichiometric mixture is 
53°9 x 10° cm’. This is indicated in Fig. 1 in relation to the equilibrium volume 
per nitrogen atom in the metal (47 x 10~*4 cm’). 

In this figure the energy level of the equilibrium state of the mixed crystal 
is indicated, together with the uncertainty in this arising from the estimated 
uncertainty in the proton affinity alone. 

To complete the estimation of the transition pressure it would seem necessary 
to know also the shape of the energy—volume curve for the mixed crystal near the 
equilibrium volume. Fortunately this is not really required. Reference to Fig. 1 
shows that the slope of the common tangent to the equation of state curves for the 
two phases is very nearly that of the tangent to the metal curve from the equilibrium 
point of the mixed crystal curve. ‘Taking the energy difference per nitrogen atom 
between the phases as 1-35, 0°85 and 0-35 eV, the range indicated from uncertainty 
in the proton affinity only, the transition pressure is found to be 140000, 100 000 
and 60 000 atmospheres respectively. It must be remembered that the uncertainty 
in this estimate is greater than indicated by this range. A transition pressure 
considerably less than 60 000 atmospheres is certainly not excluded. On the other 
hand, it is extremely probable that it is so low that metallic ammonium is the most 
stable phase throughout the main bulk of the planets such as Uranus and Neptune. 
Thus the calculated cohesive energy of metallic ammonium is not likely to be too 
large by as much as 0-5 eV as judged from comparison of the calculated and observed 
cohesive energy of the metallic sodium. In fact the calculated value is more 
likely to be too small, so decreasing the critical pressure. It is also likely, from 
comparison with sodium, that the equilibrium separation in the metal is over- 
estimated in the calculation. This would also lead to an overestimate of the 
transition pressure. However, even if the energy difference per nitrogen atom 
were increased to 1-85 eV (assuming 0-5 eV overestimate of the cohesive energy and 
the lowest estimate of the proton affinity) the transition pressure would not exceed 
250000 atmospheres. 

Although the foregoing theory cannot yield very accurate results it establishes 
that there is only a comparatively small energy difference per mole between the 
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metallic and mixed crystal phases. It is just possible that the transition pressure is 
low enough to be attained in the laboratory, but the difficulty of carrying out an 
experiment to produce the metal from the mixed crystal is so great that effort in 
this direction hardly seems justified at present. 
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RED-SHIFTS 
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Summary 


A new criterion is obtained, independent of any particular choice of 
world-model, for helping to decide between the recessional and non- 
recessional hypotheses concerning the extragalactic red-shifts. It is found 
that the shift-distance law takes the same form for any non-recessional hypo- 
thesis, provided that there is a definite correlation between shift and distance 
which is independent of time, whereas a wide variety of laws are compatible 
with the recessional hypothesis. When the assumption that the shift-ratio 
of the spectrum of a typical galaxy is independent of wave-length is relaxed, 
more general formulae are obtainable. 





1. Introduction.—The two outstanding observational characteristics of the 
extragalactic nebulae are (a) the general isotropy of their distribution, and (b) the 
red-shifts in their spectra. ‘The former property provides cogent, although not 
conclusive, evidence that the system of nebulae forms the ultimate background of 
the physical universe. The latter is a much more puzzling phenomenon to inter- 
pret, and the question as to whether the red-shifts observed in the light received 
from the nebulae are due either to (i) the Doppler effect, consequent on recessional 
motion in the line of sight, or to (ii) some progressive reddening (actual or apparent) 
of light, on its long journey through intergalactic space, is stillan open one. The 
recessional hypothesis, as usually presented, leads to a decrease in the apparent 
energy of each quantum of radiation, due to the apparent increase in wave-length, 
and also to both a decrease in the number of quanta observed in a unit interval of 
time and an increase in the apparent solid angle in which they are emitted. 
The progressive reddening hypothesis, as usually presented, results merely in a 
decrease of the energy of each quantum of radiation without the other effects. 
From this point of view the non-recessional hypothesis might appear to be 
simpler, but because the recessional hypothesis automatically accounts for the 
existence of red-shifts it has been the more widely adopted, and so far no other 
hypothesis can be regarded as more plausible. 

Moreover, although according to the progressive reddening hypothesis the 
observer attributes an actual progressive loss of energy to the photons or light-waves 
emitted by distant galaxies (unless there is a compensating progressive variation 
in the value of Planck’s constant), on the recessional hypothesis there is no such 
actual loss of energy. A photon or light-wave has no intrinsic wave-length and 
energy but only a wave-length and an energy according to a given observer, and the 
only wave-length and energy which an observer can directly assign to a photon are 
those with which he emits} or absorbs it. If two observers, A and B, are rigidly 
attached to each other, or relatively stationary, then A will assign to any length 
directly determined by B the same value as B. If, however, A and B are in relative 
motion the two values will differ. Thus, on the recessional hypothesis of red-shifts 
the differences in wave-length and energy assigned to light by observers attached 
respectively to source and receiver are due entirely to their relative motion, and 


* Received in original form 1953 November 9. 
Tt Strictly speaking, he must re-absorb it. 
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according to neither is there an actual progressive loss of energy in transit, for A will 
regard the light emitted by B as already reddened at the start of its journey to A. 
Consequently, on the recessional hypothesis there is no actual loss of energy in 
transit from one galaxy to another, whereas on the progressive reddening hypothesis 
(without invoking an appropriate variation of Planck’s constant) there will be an 
actual loss of energy since A will mot regard the light from B as already reddened 
at the epoch of emission. 

The primary object of the present paper is to present a simple theoretical 
criterion which goes some way towards deciding between these two hypotheses. 
Previous discussions have usually been based either on a particular choice of 
world-model or else on a particular hypothesis concerning the origin of the red- 
shifts. Consequently, they are all equally open to objection by an unprejudiced 
investigator who wishes to decide whether or not to interpret the shifts as indicating 
recessional motion before constructing a model. 

2. A criterion for eliminating the hypothesis that red-shifts are both time-indepen- 
dent and non-recessional in origin.—Non-recessional theories of red-shifts are of 
two types: in both, the system of galaxies is supposed to be static, apart from 
individual small proper motions; but in the one, type (a), the shifts are attributed 
to some secular variation in one of the ‘‘ constants” of nature, e.g. the velocity of 
light, the natural wave-length of a standard spectral line, etc., whereas in the other, 
type (b), the shifts are made to depend on some factor influencing light waves and 
photons in their journey from distant galaxies, e.g. the suggestion, due originally to 
Zwicky (1) and afterwards criticized by Tolman, Ehrenfest and Podolsky (2) and 
more recently by Atkinson (3), of a gravitational analogue of the Compton effect, 
or the hypothesis of some new principle of nature which automatically results in 
the “ ageing” of light in transmission over vast distances. By a suitable transform- 
ation of the scale of measurement it may happen that a non-recessional theory of 
type (a) is convertible into a recessional theory, and at the same time a previously 
variable “ constant ” of nature becomes a true constant (and, possibly, another con- 
stant becomes a variable). ‘Thus by a new natural choice of scale of measurement 
we may arrive at a recessional theory. In the case of a theory of type (4), on the 
other hand, this will not in general be possible. Moreover, specific theories of type 
(5) are usually time-independent, in the sense that the spectrum of any one galaxy 
as observed from any other does not change with lapse of time. (Some 
theories of type (a) are also time-independent in this sense.) 

We begin by presupposing that a definite scale of measurement has been 
adopted, so that an unambiguous meaning can be given to the mutual distances of 
the galaxies and to all other measures of length, but we shall make no assumption 
concerning the secular behaviour of any particular ‘“‘ constant” of nature, nor shall 
we invoke any particular hypothesis or theory of either type. We first impose the 
following condition :— 

Postulate 1: In terms of our scale of measurement, the system of galaxies is static; 
also their spectra are time-independent. 

Suppose then that O, A and B be any three points lying in this order on the 
same light-path and that these points (which we can associate with idealized 
galaxies) are all at fixed distances apart. /t will be assumed that these distances obey 
the linear addition law, so that if the distance OA be denoted by x and AB by y then 
OB will be equal to (x+y). Suppose that monochromatic light is emitted by a 
source at O and let its wave-length when emitted be denoted by Ay. If the light is 








182 G. F. Whitrow Vol. 114 


later received by A, let its wave-length on arrival be denoted by A, and on 
subsequent arrival at B by A,. According to classical physics, since O, A and B 
are all postulated to be at relative rest, we should find that A,=A,=A). Relaxing 
this condition, we consider the possibility that these three wave-lengths are not 
all equal but are related in accordance with the following condition :— 

Postulate 2: The ratios ,/Ag and Xz/Ag are independent of Ag. 

This postulate is common to the better-known theories, whether recessional or 
non-recessional. ‘There is a limited amount of direct observational evidence 
bearing on it and that is favourable. The most recent investigation, by O. C. 
Wilson (4) on the spiral NGC 4151 of type S6, shows only a very small apparent 
trend with wave-length, which may well be spurious. Between A 3400 and A 6600 
the graph of shift against wave-length does not indicate any systematic departure 
from the horizontal by more than one per cent. (Dr Hubble kindly informed me 
last year that the spectra of about twenty galaxies have confirmed Postulate 2 to 
within a probable error of about five per cent over a considerable spectral range, 
but the case of NGC 4151 was selected for publication because of its exceptional 
accuracy.) 

Since there is no observational evidence to suggest that red-shifts depend on the 
shapes and sizes of the respective galaxies, we shall assume that A,/A) and A /A, are 
independent of the intrinsic characteristics of the light-sources concerned. More- 
over, we shall assume that the spatial orientation and location of the path OAB are 
also irrelevant and we shall impose :— 

Postulate 3: The shift-ratio (wave-length received/wave-length emitted) is a 
definite function of distance travelled. 

This postulate implies that, neglecting small proper motions, the shift-ratio 
correlating any given pair of galaxies is a function of a single parameter (the distance 
between them) and is the same for all time. Specific non-recessional hypotheses 
usually satisfy this postulate. 

We therefore write A,/Ap =f (x), (1) 


where the functional form f is independent of O, A, A, and any particular epoch, 
i.e. we assume a certain degree of homogeneity and isotropy of natural law in so far 
as we postulate that the functional form of f is universal. Consequently, we shall 


> oho AalMo= f(x +9). (2) 
Furthermore, we shall assume :— 

Postulate 4: Light of wave-length d, emitted at A is received at B with the same 
wave-length d, as that of light emitted at O with wave-length 9. 

Postulate 4, which may be called the “ group property” of red-shifts, is equiva- 
lent to asserting that at a given place a light-wave or photon is completely char- 
acterized by its wave-length (according to a given observer). For, suppose that 
when the light emitted by O with frequency A, is received by A part of it is absorbed. 
Its wave-length is found by A to be A, whichever part of the beam is absorbed. If 
at the instant of arrival at A light of wave-length A, is emitted by A and travels with 
the unabsorbed part of the original beam, it will not be possible for any observer 
receiving the light subsequently to tell which part of the beam travelling on in the 
direction of B was originally emitted by O and which by A, unless some property 
in addition to wave-length characterizes light at a given place. 

Consequently, in accordance with previous postulates, we must have 


Ae/Ay=f(y)s (3) 
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and for equations (1), (2) and (3) to be compatible it follows that f must satisfy the 


functional equation 
f(x+y)=f(*)f(y)- (4) 
It is well known that the only continuous solution of this equation is given by 


f(x)=e™, (5) 
where & is independent of x. Consequently, adopting 

Postulate 5: The function f is a continuous function of a continuous variable, 
we are led to the shift-distance law 
where A= wave-length of light emitted, A + 5A = wave-length of light received and 
r = distance travelled. 

Three different possibilities arise according as (i) k >0, (ii) k=0 and (iii) k<o. 
The second case is the classical one, and the third would imply a progressive blue- 
shift. We shall ignore these possibilities and concentrate on case (i). Writing 
k=1/ro, where r, has the physical dimension of length, we have, expanding in 


series in the usual way, 
r ‘ r\2 
=a + 2\r FH cece (7) 


The first term on the right gives rise to the famous shift-distance law of pro- 
portionality, which is common to all theories. ‘The empirically determined value 
of r, is about 4x 10° light-years, according to Baade’s revised scale of extra- 
galactic distances resulting from his investigation (§) of the classical Cepheids, etc. 


in M31. Recently, Baade’s revised scale has been independently confirmed by 
Thackeray and Wesselink’s investigation (6) of RRLyrae variables in the 
Magellanic Clouds. Moreover, as Dr Baade has kindly informed me, Dr A. 
Blaauw of the Leiden Observatory has recently made a new determination of 
the zero point of the classical Cepheids (type I), using the proper motions of 15 
Cepheids derived by H. R. Morgan from meridian circle observations. He finds 
a correction of the zero point which is in good agreement with that derived by 
Baade and by Thackeray—Wesselink. Mineur (7) also claims to have obtained a 
similar correction. 

From the empirical point of view, the significant feature of the present analysis 
is the coefficient occurring in the second term on the right of (7). We draw the 
following conclusions :— 

(i) if the second-order term in the series expressing the extragalactic red-shift 
as a function of distance is equal to one-half of the square of the first-order term, 
then this shift is compatible (to the second order) with a time-independent non- 
recessional interpretation associated with our general postulates above * ; 

* The referee has suggested that to make this point quite clear it should be stated explicitly 
that, to the same order of approximation, this ratio between the second-order term and the square 
of the first-order term is also compatible with a recessional interpretation when (i) r is identified 
with luminosity—distance L (W. H. McCrea, Zs. f. Astrophysik, 9, 298, 1935), (ii) the metric of 
the model universe is 

ds* = dt®— R*(t){ dr*-+-7" d6?+-7? sin? 6 dg?}/{ c*(1 +kF?/4)*}, 
where k=1, 0, or —1, and (iii) the function R(t) satisfies, at the present moment of observation, 
the condition ts . 
° R/R= —2(R/R)’, 
in accordance with formula 71 a, O. Heckmann, Theorien der Kosmologie, Berlin, p. 64, 1942. 
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(ii) if, however, the second-order term differs significantly from one-half of the 
first-order term, then it will not be possible to adopt a time-independent non- 
recessional interpretation without relaxing one or more of our postulates above. 

We therefore find that, although a recessional interpretation is compatible with 
any particular second-order term, a time-independent non-recessional interpreta- 
tion of either type, based on our general postulates, is impossible if this differs 
from that given by (7). ‘This clear-cut criterion does not appear to have been 
previously isolated. Incidentally, it is readily verified that in the case of a uniformly 
expanding world-model, such as Milne’s original model, the formula expressing 
5A/A as a function of distance r in the world-picture, i.e. distance at epoch of light 
emission, is such that the second-order term is equal to three-halves of the square 
of the first-order term, the complete formula being (A+6A)/A=(1—2r/ry)}?, 
where r,=ct at epoch ¢. 

It thus appears that any time-independent non-recessional hypothesis, pro- 
vided it accords with our general postulates above, leads to a unique formula 
relating distance and spectral shift, the only arbitrary feature being the value of the 
constant 7p. 

3. Non-recessional theories of red-shifts—Having deduced this result from our 
general postulates, it is of great interest to see to what extent formula (6) arises in 
specific non-recessional theories in which red-shifts are independent of time. 

First, we note that it confirms the intuitive suggestion of type (6) advanced in 
1931 by J. Q. Stewart (8) that if, in the absence of recessional motion, light quanta 
travelling over vast distances become automatically “fatigued” then each quantum 
might be expected to lose energy in proportion to the distance of transmission. 
“* A lowering of frequency,” he argued, “ would then result, presumably according 
to the exponential rule v=neoH, 


where v, is the initial and v the observed frequency.” We observe that such a 
formula, x being distance travelled and H a constant parameter, is equivalent to our 
formula (6), if the velocity of light is constant, and also to the energy decay 
law —dE/dx«xE. More generally we could have obtained Stewart’s formula 
directly if, instead of following common practice in taking the shift-ratio to be 
a ratio of wave-lengths, we had regarded it as a frequency-ratio and had replaced 
wave-lengths by frequencies in the enunciation of our postulates. 

Turning to those specific theories of type (a), in which red-shifts are indepen- 
dent of time*, we find that our formula (6) arises both in (i) the theory advanced in 
1938 by Milne and Whitrow(g) according to which red-shifts could be attributed to 
a secular variation in the natural wave-length of any given spectral line, and in (ii) 
the theory advanced in 1941 by Arnot (10) in which red-shifts were attributed to 
the effect of a secular variation in the velocity of light. 

Recently, E. Finlay-Freundlich (11), supported by M. Born (12), has advocated 
a new non-recessional red-shift hypothesis} based on the formula 

dA/A= ATI, 
where A=2x 10°**cm'! deg™*, T is the temperature of the radiation field (about 
1°5 deg. K is suggested for the intergalactic case) and /is the length of path traversed. 
If the formula is time-independent, as appears to be the case, then I would suggest 


that it should be emended and written in the form 
dA/A=e*"'—1 
where k= AT‘. / . 
* These theories form only a subdivision of type (a). 

t Possibly associated with photon—photon collisions without deflections. 
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4. Recessional and time-dependent theories of red-shifts.— Although the empirical 
evidence at present available does not permit us to draw any hard and fast con- 
clusions concerning the second-order term in the shift-distance law, other criteria 
are favourable to the recessional hypothesis. A necessary, but not a sufficient, 
condition for a purely Doppler interpretation of the shifts is the independence of 
shift-ratio and wave-length; and, as we have seen, in so far as this condition has 
been investigated it has been found to be satisfied. Moreover, it has been found 
that the law of red-shifts does not operate within a physically connected system such 
as the local group to which our galaxy belongs, and this result is thought (13) to 
favour a recessional, rather than a non-recessional, interpretation. 

The technique of the present paper can be readily adapted to cover recessional 
theories* and theories in which red-shifts are not necessarily independent of time. 
The simplest recessional theory is the uniformly expanding world-model originally 
discussed by Milne. This theory is noteworthy as being a recessional theory in 
which red-shifts are independent of time. The adaptation of the present method to 
such a model leads to an interesting variation of our previous analysis, because we 
now have to consider a time-independent correlation between shift and velocity, 
and velocities do not obey the linear law of addition. 

We shall postulate that if A recedes from O with uniform velocity u, in terms of 
a scale of measurement in which the velocity of light is a secular invariant c, then 


Ay/Ao =f (u)- (8) 


Similarly, if B recedes from A with uniform velocity 7, 


A2/Ay =f (2). 
As before, if O, A and B are in line, 


A2/ o =f (w), 
where w, the relative velocity of B with respect to O, is given by Einstein’s 
velocity-addition formula 
ut+v 
~ Tue (9) 


Hence sso 
(= 3)-swys. (10) 


I+u/c 
=lo II 
a Nath 1—u/c}’ (r1) 


etc., and defining a new function ¢ such that 


$ (w;) =logf(u), 
p (wy + 2) =$(w,) +9 (wo). 
btihon tive, 


puy={ [ey (13) 


where s is an arbitrary constant, which is positive in the case of red-shifts, 


Writing 


we find that 
Hence 


and consequently 


*In the case of recession red-shift means shift for each particular galaxy and not for each 
particular distance. 


13 
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In Special Relativity, and in Milne’s original model, we have s=1, corre- 
sponding to the postulate that atomic vibrations and photon standard frequencies 
are uniform in the scale of time employed.* The more general formula here 
deduced is, however, explicable on the theory of clock-regraduation. Suppose 
that O and A part company at time zero and are equipped with congruent clocks. 
If t’ denotes the local time of reception by A of a light signal emitted by O at 
local time t), we have 

ty’ = ato, (14) 


where «= J (= “) . If O, A regraduate their clocks by the rule T=t’, 





where & and s are positive constants, 
T,' =«°Tp, (15) 
where 7,’ is the new time of reception by A and 7, the new time of emission 


from O. In the new scale of measurement, O and A are still in uniform relative 
motion but their relative speed of recession is now U, where 


1+U/c_ (1+u/c\* 
(FS) , (16) 


If atomic clocks keep 7-time then the Doppler-effect formula for light travelling 
between O and A is given by 


meen?) e en 


where u is the relative velocity of O and A according to the original time-scale ¢. 
This is the formula (13) obtained above. 

We deduce that the postulated uniform expansion of the world-model does 
not of itself suffice to isolate a unique scale of time. Conversely, however, we 
can give a definite meaning to saying that the model is in a state of uniform 
expansion if our clock is already given, e.g. as the uniform time-scale of atomic 
vibrations and photon frequencies. 

Turning now to the world-models of General Relativity, we find that they 
are characterized by non-uniform expansion and by red-shifts which are not 
independent of time, whereas both in the non-recessional theories discussed 
above and in Milne’s uniformly expanding model the shift-ratio of light from any 
one galaxy as observed from any other was independent of epoch. We now relax 
this latter restriction (Postulate 1) and replace Postulate 3 by :— 

Postulate 6: The shift-ratio is a definite function of epoch of emission and epoch 
of reception. 

Introducing a definite set of proper epochs, 79, 7,, T2, to denote the respective 
instants (local times) at which a photon passes O, A and B, we now have, in place 
of (1), if we still invoke Postulate 2, 

A1/Ao =f(7 1 To)s (18) 
and similar formulae for A,/A, and A,/A,. 

We also impose Postulate 4 (the “group property”) in the following more 
stringent form, viz. 

Postulate 4 (revised): Light of wave-length d, emitted at A at 7, is received 
at B at tr, with the same wave-length r, as that of light emitted at O at tr. with 
wave-length xy. 








* Formula (13) with s=1 leads to the formula correlating shift and apparent distance mentioned 
at the end of Section 2 above. 





No. 2, 1954 On the interpretation of the extragalactic red-shifts 187 


We now find that our previous functional equation (4) is replaced by 


S(7 25 70) =A» To) f(a 74)- (19) 
To solve this equation write 


SE, 1) = exp o(€, 0), 
whence (72, To) = $(71, To) + P(72 71); 


the general solution of which is 


H(E, 0) = WE) — Wn), 


where denotes an arbitrary function of a single variable. Hence, we have 


Ay/Ao = exp {(71) — 179) };, (20) 
and it is clear that all shifts will be to the red provided that (7) is a monotonic 
increasing function of its argument. Introducing a new monotonic function 


R(r) = exp (7), (21) 
we can write (20) in the form 


Ay/Ag = R(71)/R(70).- (22) 


If Ay, A, are measured in a scale of length associated with a constant value c 
of the velocity of light, two different kinds of model can give rise to formula (22) :— 
(i) O,A... are all relatively stationary, i.e. the model is static. In this case 


To=71—7/6, (23) 
where + is the constant distance between O and A. Hence 
tages (24) 
Ao R(t; — ric} ‘ 
In general, this shift-ratio observed at A will not be independent of epoch 7, 
at A, unless 
R(r) ce”, (25) 
where a is a (positive) constant, in which case we have the situation already 
considered in Section 2. 
(ii) O, A,... are all in relative recessional motion, i.e. the model is expanding. 


With each monotonic increasing function R(r) we can, if we wish, associate a 
homogeneous expanding relativistic model of metric 
dy” + 7? dé? +77 sin? 6 dd?) 

(1 +h7*/4)? J’ 
where R= +1, 0, —1. It is well known that the red-shifts arising in this model 
are governed by formula (22). 

5. General analysis of the law of red-shifts.—Although, as already remarked, 
present observations appear to confirm the view that red-shifts are independent of 
wave-length, as soon as we consider the possibility of a non-recessional inter- 
pretation of the red-shifts the theoretical necessity for a constant shift-ratio 
throughout the entire spectrum becomes less evident. Although we have at 
present no reason for rejecting the recessional hypothesis, it is of theoretical 
interest to discover what limitation is placed on the possible variation of wave- 
lengths in the spectra of different nebulae, assuming that they are relatively 
stationary and that their spectra are time-independent but relaxing Postulate 2. 

13* 





ds? = c? dr? — {R(r) 4 
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First, we shall consider the consequences of adopting Postulate 1 and the 
following :— 

Postulate 7: The wave-length of monochromatic light received at any point A 
is a definite function of the wave-length of light emitted at point O and the distance OA. 

Thus, with our former notation, we write 


Ay =f(%, Ao), (26) 
where f is a function which may involve A, in a non-linear manner. Similarly, 
Ag=f(x+y, Ao), (27) 

and, appealing to Postulate 4, we also have 
A2=f(y,A1)- : (28) 


Replacing A, by variable z we deduce that the function f must satisfy the functional 
equation 
AIA% B)}=H{x + ys 3}. (29) 
To determine the form of f satisfying (29), it is convenient to write f(x, 3) 
as W,(z). We then find that 


ty {b(2)} 7 $2+y(2) i $alp(2)}, (30) 
i.e. Y, and ys, are commutative functional operators. Assuming, as a generaliz- 
ation of Postulate 5, that these operators form a continuous group, for all x, y>o, 
it follows (14) that there exists a functional operator 0, independent of x and y, 


such that 
$2) =F %a(x)A(z), (2) =F a(y)0(2), (31) 
where «(x) satisfies the functional equation 
x(x)a(y) = a(x + y). (32) 
Hence 
a(x) =e, 


where & is a constant parameter, and we deduce that 
Ay =O {e6(Ap)}. 
Consequently, there exists a functional form @, independent of O, A and Aj, 





such that 
8(A,)/8(Ag) = e*. (33) 
Hence, in place of (6), we have the shift-distance law 
O(A+8A) 
my ~*" (34) 


where A + 6A = wave-length of light received, and 7 = distance travelled. Expanding 
to the first order, we have, in place of (6), 


+ =k at (35) 


To obtain a red-shift increasing with distance it is only necessary to postulate that, 
assuming 6(A)>0, its derivative @’(A)>0, i.e. (A) is a monotonic increasing 
function of A. Provided that @(A) is proportional to a positive power of A we 
obtain formula (6), giving 5A/A the same throughout a given spectrum. In the 
general case, however, relaxing Postulate 2 but retaining Postulate 1, we no 
longer have this condition holding but instead there is still some function 6 such 
that 6(A + 5A)/6(A) has this property. 





No. 2, 1954 On the interpretation of the extragalactic red-shifts 189 


Finally, we consider an even more general case in which we impose only 


Postulate 4 (revised) and 
Postulate 8: The wave-length of monochromatic light received at any point A 
is a definite function of the wave-length of light emitted at any point O and of the 


epochs of emission and reception. 
Adapting our former notation, we now have 


Ar=S(71 700)» Ag=Sf(T2 71» Aa), (36) 
and similarly, appealing to Postulate 4 (revised), we have 
Az =S(T 2» To» Ao): (37) 
Consequently, f must satisfy the functional equation 
fi t2 Ta» A719 Tos Ao) } =f 729 Tos Ao} (35) 
To solve this equation, we temporarily replace 79, 71,7, by symbols x, y and z, 
not to be confused with our previous use of these symbols, and write A for A). 
Writing x, y, z as suffixes, we have 


i i a a(A)} =f, 2(A), 


i.e. omitting commas, the functional operators f,,, f,, and f,, are related by the 


law 
Sadun =fax- (39) 
To find the form of f,,, it is convenient to introduce the operator f,,. The 
equation on the left of (36) gives the wave-length of reception as a function of the 
epoch of reception, the epoch of emission and the wave-length of emission, in 
that order. Imagine now a reversal of the flow of time, so that all processes are 
seen in reverse. ‘Then we should have 
Ao =F (To 71 Ar)» 
since Ay would now be the wave-length of reception, A, the wave-length of 
emission, 7, the epoch of reception and 7, the epoch of emission. Thus, if we 
write A, =fy2(Ap), we can also write Ay=/,,(A,), i.e. the operators f,, and f,, are 
functionally inverse. Hence 
Sox =f ey a, Sua =fyx a, (40) 
Consequently, (39) can be written as 
; a . =f... (41) 
Since y does not appear on the right-hand side of (41), it must play the role of a 
dummy suffix on the left and we can write on this side 


| t=, fey= 2 , (42) 


See sek pe . (43) 
Replacing x by 79, y by 7, and z by 7, we find that 

(725 Ag) = (79, Ag), (44) 
and we deduce that a function ¢(7, A) exists so that the wave-lengths A assigned 


to a given photon at different points on its path are related to the respective 
proper (local) epochs r at which it passes these points by the law 


(7, A) = constant. (45) 


so that 
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Our previous formulae are particular cases of this general law, which holds 
even when O, A, etc. are in relative motion, provided that 7 is always interpreted 
as proper time. ‘Thus, if the variables 7, A are separable, so that 


UA, 7) =4(7)8(0), 
we find in particular that 
(i) ¢(7)=e° “" corresponds to (34), provided r=c(r,—7»), 

(ii) 6(A)=A corresponds* to (22), 

(iii) 0(A) =A, (7) =e-**” corresponds to (6), provided r= c(t, —7»). 
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A STUDY OF THE EQUILIBRIUM OF GLOBULAR CLUSTERS 
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Summary 


This paper examines the theory of the distribution of density in globular 
clusters by a method which gives prominence to the velocity distributions of 
the stars. The technique reproduces standard generalizations by Jeans and 
Eddington and leads to the well-known “ isothermal’”’ distribution when 
the velocity distribution is Maxwellian and the stars all have the same mass. 
The general equation for the density and potential is given for any isotropic 
velocity distribution, where the stellar masses are not all the same. This 
equation is solved numerically (1) for stars of equal mass, having a “‘ trun- 
cated Maxwellian’ velocity distribution; (2) for a full Maxwellian velocity 
distribution, but with the number of stars of mass m proportional to mPe—m, 
The solution is then connected with Emden’s equation of index —(p+2) ; 
and (3) for a truncated Maxwellian velocity distribution with the number 
of stars of mass m proportional to mPe~*™. The case p=1 is examined in 
detail and the distributions of the stars of various masses discussed. 





1. The idea that the stars in the galaxy, more particularly those in relatively 
dense aggregations like globular clusters, influence each other like molecules in 
a gas, is frequently to be met with in astronomical literature, and indeed coloured 
all early attempts to account for observed distributions of stars in clusters. The 
analogy with a gas, if it can be sustained, provides a powerful short cut to the 
solution of the equations of motion of a large number of gravitating bodies: 
and the analogy appears likely to be valid when one considers that the actual 
behaviour of gases can be well represented by considering the molecules as point 
centres of force. But the analogy is useless unless collisions among the stars 
are sufficient to set up equipartition of energy, and not very valuable unless the 
mean free path for stellar encounters is small compared with the dimensions of 
the stellar aggregation considered. Calculation shows that the mean free path 
is probably large compared with the effective diameter of a cluster, so that in 
clusters the analogy with a gas is an idea that cannot be pushed very far. 

It is the author’s point of view that a satisfactory theory of the velocity distri- 
bution can only be built up from an elaborate theory of relaxation: and that 
until this is done progress can only be made by constructing models for comparison 
with observation, that is to say, introducing arbitrary mathematical restrictions 
into general solutions of the problem, such as, of course, has often been done 
before. It is here supposed that it is best to introduce these restrictions in the 
velocity distribution, and to abandon both gas analogies and “‘ models” of the 
potential. 

2. According to Plummer and von Zeipel the distribution of star counts in 
globular clusters is closely represented by the fifth polytrope, and in an important 
paper Eddington, who supposes that the density follows the same law, examines 
a number of hypotheses without arriving at a justification for this particular 
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polytropic distribution. Attention has now rather shifted to globular galaxies, 
which are supposed to follow a different law of distribution. For example, 
Hubble showed that in these the projected density resembles that of a projected 
isothermal gas sphere—but with a difference. An isothermal gas sphere has an 
infinite mass, and to obtain a model which was capable of comparison with a 
globular galaxy Hubble simply cut off the isothermal gas sphere at a suitable 
distance from its centre. Recently de Vaucouleurs has introduced the very 
important modification into Hubble’s model of populating it with a variety of 
stars with very different masses and luminosities. 

The mathematical result that an isothermal gas sphere has infinive mass may 
be compared to the physical idea that it is impossible to set up a finite isothermal 
body unless one provides it with an ideally non-conducting boundary—a perfectly 
lagged jacket. Ina real sphere of gas, held together by mutual attraction of the 
particles and extending into space with, say, an exponential decay of density, 
there must be a temperature gradient: and so if the cluster stars behaved like 
molecules in a gas we should expect to find a high kinetic energy per unit mass 
in the centre, and a transfer of this energy by collisions from inside to outside, 
where a lower kinetic energy per unit mass would be found. Conduction only 
takes place from a source or sources to a sink, and in the case of a cooling sphere 
this sink is outer space: so we should expect the sink in the case of the “‘ gaseous”’ 
cluster to be the escape of high velocity stars from the gravitational control of the 
cluster. This “cools” the outer shell of the cluster by systematically lowering 
the average kinetic energy. Instead of an isothermal sphere, we are led to the 
model of a cooling sphere. 

Doubtless this model would serve a useful purpose if the mean free path of 
the stars was small compared with the dimensions of the cluster, but this does not 
appear to be the case. Instead of the stars exchanging energy within a small 
volume in the cluster and setting up an identifiable ‘‘temperature’”’ in this 
volume, and other ‘‘ temperatures” in other volumes, the stars seem able to pass 
many times through the cluster without sensibly exchanging their energy with 
that of other stars, in effect keeping the “‘temperature”’ of the cluster well mixed. 
It seems safer to assume that transfer of stars from one spherical shell to the 
next takes place with no collisions at all, and can be worked out from their gravi- 
tational orbits. In this event the only effect of collisions is to impose some 
restriction on the velocity distribution in the cluster as a whole. When this is 
done, however, some remnant of the analogy with a cooling sphere can yet be 
found, for while it may be shown that if the velocity distribution is anywhere a 
Maxwellian distribution, it must be a Maxwellian distribution everywhere and 
the cluster must be isothermal (i.e. has everywhere the same parameter in the 
Maxwellian velocity distribution): yet if account be taken of the absence of stars 
with high velocity, the kinetic energy per unit mass decreases outwards from the 
centre of the cluster. 

3. Consider a spherically symmetrical cluster of stars and consider two shells 
in the cluster so close together that the motion of a star from one shell to the other 
can be described with sufficient accuracy by a Keplerian orbit. It is assumed 
that the probability of a collision in passing from one shell to the other is so small 
that the effect of collisions can be neglected. 

Now consider the class of stars whose velocities lie within a given narrow 
range in one shell (and a corresponding range in the other). Then if the velocities 
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are distributed isotropically in both shells, the ratio of the densities of the stars of the 
given class in the two shells is equal to the ratio ef the velocities themselves. 
The equation of motion of any star is 


d*u m 
qe t+ Y= pp (3-1) 


where , is the constant of gravitation multiplied by the mass inside a sphere of 
radius Iju. Ina restricted range of values of u (from u, to wu.) we may suppose 
to be given with sufficient accuracy by 


B= Hy +(u—m)e’, 

and then the solution of (3.1) is 
LL, 

u= [3 (1 +ecosn8) (3.2) 

if nm? =1—p'/h? and py =(", — up’) /n?. 
Now consider a star at 7, moving with a velocity 7, in a direction making an 

angle « with the outward radius vector. Then h=r,v, sina and, if the motion 
is given by (3.2), the constant e is found from 


ht 2h? —sh®v,* cos* x 


2 

e-1=— - 
2 2 2, 2 

Mol) Bolt no 


The outward component of velocity at r, is (dr/dt),=v,cos«. At a point rz 
close to r, we have 


(dt/dr),,/(dt/dr),, = {a® — b® — sin? «}!?/x cosa, (3.3) 


where x=/7y/7,, @?=x?+2x(1 —x)u,(uy—Uyp’)/0,?, 6? =u,2u'(x? —1)/v,?. When 
a® —b? >1 all v, stars at r, reach r,, whatever the direction «. If 1 >a?—b?>0 only 
those proceeding in directions sin~! ,/a?—6® >|«| >o reach (or come from) 7, 
and if 6? >a*, nov, stars at 7, reach 7. 

Now the ratio of the number of stars with velocities from v, to v,+dv, at r, 
to r,+dr to the number with velocities between v, to v,+ dv, in the shell from 
r, to r.+dr (subject of course to v,? — vg? = 244(u, — Ug), Vy dv, = Vy dvq) is propor- 
tional to the times spent in the shells, averaged over all values of «. If the 
velocities at 7, are distributed isotropically the fraction of stars with directions 
between |x| and |a+dz| is sin|a|da. Let m(v,,7,) dr dv, be the number 
of stars between 7, and 7,+dr and with velocities between v, and v,+dv,, and 
let n'(vg, 72) dv, dr be the number of stars between r, and r,+dr with velocities 
between v, and v,+dv, and which pass through r,. Then from (3.3) 








N'(Vg, 1%) dv, _ [ x cos «sin x dx 


n(v,,7;)dv, J 9 {a®@—b?—sin?«}!?’ (3-4) 


the upper limit / being given as follows: 
If >a l=o 
I>a*-—b?>0 l= sin~! (a? — 5?)!? 


a®*—b? >1 l=n/2. 
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Accordingly 
n'(W ya) de _ if B>a? (5-42) 
ny, 7,) dv, 
= x(a? — b?)!2 if 1>a*-—b?>0 (3-42) 
= x{(a?— 5?)'2—(a?—b?-—1)'?} if a?—b? >I. (3-43) 


Again, remembering that 


Ug? — Uy? = 2g — Uy) = 2py(Ug— Wy) + L'(Ug— uy)’, 
we have, after a little reduction, 
xv3?/0," =a? — b?. (3-44) 
Equation (3.44) also shows that if b?>a? the stars have imaginary velocities 
at 7, i.e. do not get there. If a?=b5? the stars are at apocentre at 7r,. If 
o<a?—6? <1 all the v, stars at r, pass through 7,, but if a2—5?>1 some of the 
V, Stars at r, (i.e. those with nearly circular orbits) miss 7,. We enquire what 
value of ~ at 7, is such that the path just touches 7. 
We have h=r,v, sina, =rgv2,sina,: if the path touches r, then sina,=1 and 





. rv, j%&Iv 
TeV 4 V2 
and using (3.44), 
sin? %» = (a? — b?)-}. (3-5) 


If the velocities of the v, stars at r, are distributed isotropically the fraction 
making directions between « and « + d« with the radius vector is sina dx: hence 
the fraction of v, stars at 7, which pass through r, (in the case a?—5?>1) is 

a/2 
| sind which is 
(t—cosa), or 1—{(a?—b*—1)/(a?—b%)}12, (3.61) 
Now let (v2, 72) dr dv, be the number of stars between 7, and r, + dr with velocities 
between v, and v,+dv, whether the paths pass through r, or not. ‘Then 
N(Vg,%2)=N'(V, T2) if a®*—b?<1, (3.62) 
and from (3.61) 
N( Vg 72) = N' (Vg; To) X (a? — b®)*?/ {(a? — b?)42 — (a?@—b?—-1)'*} if a*@—b?>1. 
(3-63) 
Combining (3.42) with (3.62) and (3.43) with (3.63) we get 
Ma» F2)dP2 _ qe payne 
ne, r,)dv, = x(a —5 ) (3-7) 
in all cases (except b? >a? when n(v,,7,)=0). Using (3.44) we find 
N(Va,1%)dV, 4g Va | T,*v, 
n(v,,7;) dv, oie % 7°, (3-8) 
The number densities v(v, r) are related to the total numbers in the shells n(v, r) 
by n(v, r)=47r*v(v, r) and therefore 
V(Ug,%q) dU, Vg 
V(4,7;)dv, -%” (3-9) 
the result sought. It has been shown to be true if the paths between 7, and 7, 


are governed by (3.2) and if the distribution of velocities is isotropic at both 
r, and 75. 
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4. Let v,,(v,7)dmdv be the number of stars per unit volume at r with 
velocities between v and v+ dv and with masses between m and m+dm. If stars 
in this group move from 7, to 7, the velocities change from the range v, to v, + dv, 
to the range v, to v,+dv2, where 

027-0? =2(d2—91), U2 dvg=0, 42, 


and ¢ is the gravitational potential at r. 
Let 
Vin(U, 7) = Ofn(P, &*). 


Vn( P71) = Pf Py @1”)- 


But from (4.2) and (3.9) with v,dv,g=v, dv, 
Vin( Vay ¥2) = Ve"fin( Pr) V1”): 
But Vin( Vy V2) = Uo"fn( Po, U2); 


so that Fin $1» 01) =Sin( $2» V2”). 


In bo» v2") = falas v,") + (2-44) a 


and since v4” — v,? = 2(¢. — ¢,) equation (4.3) gives 


Of n(h, v? Of, (d, v? 
Lalo), 4 Yah 2) (4.4) 


This is equivalent to Jeans’ Theorem*, as (4.4) is satisfied by /,,(¢, v*) = F,,(v? — 2¢). 
Notice, however, that Jeans’ Theorem does not hold for an ellipsoidal velocity 
distribution. In particular, it is not valid for the system of radial oscillations 
proposed by Belzer, Gamow and Keller.+ 

We can now show that if the distribution of velocities is Maxwellian at any 
point it is Maxwellian everywhere in the cluster. If the distribution is Maxwellian 
at some radius 7, 


Then in particular 


Now of, 
+ (v_”— v,?) 522 ’ 





H(i» 71) = Noa Pa)(Bm)®2 S02 exp {— Bmo,*} do, 
where £ is independent of m and v,._ Then 
ful ?12)=No(t(Bm)?®5- exp {—Bme,%}, 


and (4.4) is satisfied if N,,(r)=C,, exp {2Bm(¢ — 9)}, where C,,, depends only on m. 
The distribution anywhere is therefore given by 


Hn (s 7) =Cip(Bm)®? S- 08 exp {Bm(2d — by —0*)}, 


showing that the velocity distribution is everywhere Maxwellian and has every- 
where the same parameter f, i.e. everywhere the same “‘temperature”’. 

Now C,,, is independent of ¢@ and N,,=C,, exp {2mB(¢—5)}, so that Cc. s 
proportional to the mass function at the centre of the cluster, or C,,=A fol), 


* J. H. Jeans, M.N., 76, 78, 1915. 
+ J. Belzer, G. Gamow and G. Keller, Astrophys. 7., 113, 166, 1951. 
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where N, is the total number of stars per unit volume at the centre and £,(m) the 
central mass function. ‘The central density py is given by 


Po= No| m£ (m) dm 


and p, the density at any point, by 


Po= |" m£,(m) exp {2Bm(¢ — bo)} dm |" mé,(m) dm. 


Only if the stars all have the same mass do we recover the ‘‘ equation of the isothermal 
gas sphere”, namely p=p,exp(const x 6—¢ 9). The star cluster on which the 
stars move in and out from the centre with small chance of a collision, leaving 
equipartition to be set up as the result of many journeys through the centre, 
differs from the isothermal gas sphere setting up a distribution for a mean molecular 
weight. 

5. Eddington’s relation.—Let E be the kinetic energy per unit volume, and let 


E(m) dm be that part of it possessed by stars with masses between m and m+ dm. 
Then 





2E(m) = j. mv*f,,(p, v?) dv. (5.1) 


In this integral the upper limit L may be infinite, but we shall want later to study 
cases where L is finite, and L?—2¢ is constant. Differentiating (5.1) with 
respect to r 











ide m) “L Of, ab = * 
ee ee 
But by (4.4) we have 
ph Of n + ee = 4 Am 
| Pond a6 d du= 2m "de! dv 
d is 
aa ait = L?) a | ‘ v*F(P, v*) ac}. (5-3) 
dL di ; 
But L — = -¢, so that (5.2) and (5.3) give 
dE de 
2) 5 an mo%f,,(¢, 08) de. (5-4) 


Now »,,(v,7r)=v*f,,(¢, v?), so that mi v*f,,(¢, 0?) dv=p(m) the component of 


density due to stars of mass m at r. ‘gree 


dE(m 
sad ae 2 o(m m) 





and integrating with respect to m 


pS. (5-5) 


. SS 
wie 
en 


This equation was given by ito | 


* A. S. Eddington, M.N., 76, 581, 1916. 
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Eddington’s equation may be compared with the equation of hydrostatic 
equilibrium, because E, the kinetic energy per unit volume, has the same dimen- 
sions as pressure : and in the case of a gas dominated by collisions E= #p. Hence 
if we define a virtual pressure P by P= §E equation (5.5) gives 

dP=p dd, (5.6) 
the familiar equation of hydrostatic equilibrium. Since gravitating stellar 
systems must satisfy (5.6) and Poisson’s equation, their equilibrium is controlled 
by two equations common with those of a single star. The difficulties in the 
cluster case lie in the concept of temperature. 

6. Before considering more general cases we examine in this section and the 
next the special case in which all the stars have the same mass m. First, if the 
velocity distribution is Maxwellian, 


P = po Exp {2mB ($ — $o)} = po Exp 2j” ($ — $5). (6.1) 
It is convenient to set p/py»=7 and 2j/7(¢—¢))= —¥%. Then n=e-’. 
Poisson’s equation is 
1 d/,d 
Rar (» +) a 
where I is the constant of gravitation. If we put z=r x (87I"poj*)"” this becomes 
S (= +) =72*= 2%", (6.2) 
the equation of the isothermal gas sphere with constant mean molecular weight. 
Solutions of this equation are tabulated in Emden’s Gaskugeln, Eddington’s 
Internal Constitution of the Stars (p. 90) and in Chandrasekhar’s Principles of 
Stellar Dynamics (p. 234). The projection has been given by Hubble.* 

It is well known that it can be objected against the “isothermal” solution that 
it has an infinite mass. We can certainly suppose that an initial condition for a 
globular cluster is, that it has initially a finite mass. Let us suppose further that 
another initial condition—before any time has elapsed for relaxation (equipartition 
of energy)—s that it has a finite actual plus potential volume. Indeed this cannot 
be escaped for a finite number of stars. If any of them have the velocity of escape, 
they move off and do not remain in the cluster: they need not be counted as cluster 
members at all. Of those that remain, one has the greatest energy per unit mass, 
and may be supposed to move ultimately to its apocentre. Until collisions 
modify the velocities, all the cluster stars remain inside the sphere whose radius is 
this apocentric distance, say R. There is of course nothing to prevent this being 
very great, but on the other hand, if we consider a cluster as nothing more than an 
elaboration of a double star, triple system, etc., there is no difficulty in supposing 
that in some cases at least systems with a large number of stars start out with 
initial orbits all contained within a finite volume—with a radius small compared 
with the distance from one cluster to the next. 

If this is the case the velocities at the centre start out with an upper limit—the 
velocity of fallin from radius R. As time goes on collisions, taking place for the 
most part near the centre of the cluster, will modify the velocity distribution and 
will tend to set up a Maxwellian distribution in which all velocities are represented : 
but this may be a very slow process as far as the high velocities are concerned be- 
cause a high-velocity star is more likely to lose energy in a collision than gain it and 

* E, Hubble, Astrophys. F., 71, 231, 1930. 
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especially because once a star has got velocity anywhere near the velocity of escape it 
spends almost all its time well away from the centre of the cluster and experiences very 
fewcollisions relative to the number experienced by a low-energy star which spends 
its time near the centre. In other words the time of relaxation for high velocities 
is much greater than the time for average velocities. Thus even if the time of 
relaxation calculated for the centre of a cluster is of order 10° years, a cluster may 
be supposed to be 10” years old or more, and yet to have had very little change in 
the velocity distribution of its most energetic stars : so that if it started out with 
a limiting radius, it still has a limiting radius, not greatly changed from the original 
one. We may therefore suppose that there is some interest in the case 


S(¢, 0?) = F(v,2-—v*) ~~ for v?¥<v,? 
=0 for v?>v,?, 


where v,?=2¢+¢, so that v,* vanishes at some r= R, say. Thus Eddington has 
discussed the case 

F(v,?—0*)=C (2? - 0)’ (6.3) 
and also the “truncated Maxwellian” distribution 

F(v,? — v*) = Cexpj?(v,?—v*). (6.4) 
As Eddington showed, (6.3) leads to a polytropic distribution, since in this case 

v(v,r) = Cv® f (¢, v) = Cv? (v,? — v2)’, 
so that 
Uy : K. ‘ = , 
o(v,r)=C | v%(v,2 — v*) dv=}3Cv,"+8 | (1 —t)'t!? dt =4Cv, ***BA+1, 3), 
/0 0 
i te (6.5) 

where B is the Beta function. Then from (6.5) p=const x (¢’)****, where ¢’ isthe 
gravitational potential adjusted to vanish at r= R, the outside of the cluster. This, 
with Poisson’s equation, gives the polytrope of index A+ $. However, no physical 
argument has been advanced for the form (6.3). 

7. The truncated Maxwellian distribution.—Only a study of relaxation can deter- 
mine whether or no the truncated Maxwellian velocity distribution (6.4) is a 
good approximation to the actual velocities in the state of affairs discussed in 
Section 6. In this section we adopt this velocity distribution as a model, that is, as 
an arbitrary or plausible simplification introduced into the mathematical descrip- 


tion of a quantity, the exact expression being unknown. We consider in this section 
the case where all the stars in the cluster have the same mass. Now from (6.4) 


v(r)=const x exp { j?v,?} | g? exp { —j?v?} dv. (7.1) 
0 


“zr 
As the integral | x* exp { — x} dx will occur often it is convenient to use a symbol 
0 


for it, so we set 
= 


er 
| x exp { —x?} dx=erf, x. 
Valo 


As x-> ©, erf,.x-1. In this notation (7.1) becomes 


v(r)=const x exp {j?v,?} erf, (jz). (7.2) 
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Here z, is the limiting velocity atr. It is related to vp, the limiting velocity at the 
centre, by 
v,* = 2¢ = Uo" 7 2p. 
so that, since 2j7(d—do) =, J? (Up? —v,") =¥. 
Write j2v,2=k : then j?v,2=k—y and (7.2) becomes 
v(r) =const x exp {—#} erf, 1/(k—Y). 


Accordingly » = p/p, = exp { —} erf, \/(k —)/erf, \/k and Poisson’s equation (6.2) 
gives 


d/dz (2%) = 2* exp {—p}erf, /(k—y)/erf, Vk (7-3) 
with the boundary conditions ys =0 and % =0 at z=0. 
The kinetic energy © per unit mass is given by 


= ; o v‘ exp { ~j%3} do/ (" vexp {—j*v"}dv= =aerty V/(k—-d)/erf, V/(k-¥) 
Jo Jo 2 


3 , 
-P {1— $(k—p)>? exp {—(k—-)}/-V/7/4 er, V(k—-Y)}, 

so that © decreases from approximately 3/4)? at the centre of the cluster (ys =0) to 
zero as ys —k at the outside of the cluster. If © is considered as the counterpart of 
temperature, the truncated Maxwellian velocity distribution in a cluster presents 
some resemblance to a cooling rather than an isothermal sphere. 


Since, by (6.2), p/p» = = < (z*s) the mass inside radius 7, m(r), is given by 


m(r)= | 4npr®dr=4npp (8x0 pgj?)-8? 2%. 
“0 
Hence if (zy), is the value of 2 at the outside of the cluster, the total mass M of 


the cluster is given by 
M = 47 (2*p).0/(S71j?)>* po. (7-4) 


Now the velocity of escape at the outside of the cluster is u,. = \/(2l' M/R) and R, 
the limiting radius is z,, x (87Ipgj?)-¥7._ Then by (7.4) 

jru*.. = (27) 0/%ac» 
the subscript oo denoting in every case the final value at the end of the solution, 
p=0. The velocity of escape at the centre of the cluster is vu), where 

Jug? =j7u?,, +k. 
But v, the limiting velocity (for truncation of the velocity distribution) at the centre 
of the cluster is given by j?7)?=k. Hence 

te 

Up? (22s) + REx” (7-5) 
which shows how the degree of truncation below the velocity of escape can be 
found for any value of the parameter k when the end values z,, and (2*#),, of the 
solution have been calculated. 

Notice that this is the degree of truncation relative to the velocity of escape and 


not relative to the parameter of the Maxwellian distribution. In fact, since 
v? = 2/37? for a full Maxwellian distribution, 


u2/0* = $k. (7.6) 
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8. Returning to the general case where the masses of the stars are not all the 
same, we may adopt Jeans’ solution of (4.4) and write 
Vm (0, 7) =v? F,, (v? — 24), (8.1) 
where F,,, (v? — 2¢) is related to the central mass function £,(m) by 
£,(m)= | “v® Fy, (e?— 249) do / | [°c Fy (0-269) do dm (8.2) 
0 0-0 
and the limit v) may be finite or infinite. The density p is given by 


°C rT 
p= | |” me? F, (08-26) do dm (8.3) 


and as 7, is a function of ¢ the density is known as a function of ¢ when the function 
F.,, (v® — 2¢) is settled—by adopting a model of the velocity distribution, in default 
of a real determination. Equation (8.3), together with Poisson’s equation and 
appropriate boundary conditions (6 = ¢y, 6 =0 at r= 0), gives a general solution of 
the distribution of ¢ and p with 7, when the function F,,, (v? — 2¢) is settled. 

g. To show how equation (8.3) can be applied in practice we consider the 
model 


Fa (v? wie 2¢) = Cy exp { om (v? aa 2¢ as do)} for v? <U,? —2 (¢ a do) 1) 
=0 ot ey? 2(b— dy). - 
This form includes the Maxwellian and truncated Maxwellian distributions, 


obtained by setting «,,= $m, the full Maxwellian distribution being of course 
Vg= 00. From (8.2) and (9.1) we have 


Eq(m) 7" Cn Lm - erf, (U Ven) | ‘ Cr Sm - erf, (u% V %m) dm 


and in practical cases we may take vy to be so large that erf,(vp1/«,,)=1 without 
serious error, and then 


£ (m) = Cu om ' 2} f Co —" dm. (9.2) 


In the first case, where «,,, is independent of m, =«, say, from (8.1) and (9.1), 


Vn (Vv, 7) =C,, exp { —2%(d—go)}. v%exp{—av?} for v<y, 
so that 


wi 
Vm (r) - | P Vm (z, r) dv=a*? Cn exp { — 2a (¢ = o)} ” + erf, (a? . 2%); 


leading to 
p/po= exp { — 22( — do)} . erfy (x? v,)/erf, (x? v9), (9.3) 
which is the same as (7.1) if 2a(6—¢))=—%. Tf, then, x is independent of m, it 
makes no difference whether the stars have the same mass or not. 
In general, where ~,,, is a function of m, 


m 


r2 rrr 4 00 
p= ‘ f Vm (v, 7) dv.mdm= Teh, MC yp%p,-8? erfy (a222,,04) EXP {2%m (b—by)} dm. 


(9-4) 
If there is equipartition of energy among the stars with small velocities «,, = Bm, 
and, writing Bv,?=w?, 


p=p-**. =. | m-*? erf,(w/m) C,, exp {28m (¢— d9)} dm. (9. 
7™Jo 
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Now let n= | mC, dm | | m-"2C, dm, 
0 Jo 
2pm (> — do) = —¥, 


m/m=p, C,,dm=C,,du. 
Then p= -ma(fimy |" erhyfwV/ (nim) ™*Cyexp{—nH}du. (9.6) 


So long as w is large this reduces to the isothermal for stars of varying mass, 
namely 


p=constx |e 4*C,exp {— ph} dp. (9.7) 


From (9.7) we see that when ¢ is small 


ploo=1— HJ w*Cudu/ |” ™® Cyd 
= I-y, 


by virtue of the definition of m. With this choice of m the solution starts out at the 
origin in the same way as the isothermal for stars all of the same mass, in which 
pP=poeXP{—Hj- 

10. At this point it is necessary to specify a particular case of the mass 
function. (For anything that we know, any distribution of masses could perfectly 
well occur inacluster.) A case very well suited for general examination is 

C,, = const x m* exp { — bm}. (10.1) 
(As we see from (9.2), the central mass function £)(m) is then approximately 
proportional to m*-9 e~™, so that if a>} we shall get a likely representation of the 
mass function, as it is known or supposed.) From (10.1), 

‘itp m= (a+4)/6, 
and reintroducing ».=m/m, 
C,,= const x “exp {—(a+ })u}. 

Combining (9.7) and (10.2) gives 


7 =p/po=const x (a+ $+y)-@r2® | s-US e* dx 
0 


and since » = 1 when ¢ — 0, at the centre of the cluster, 


a+} a+1/2. 
7 iets) 


Let y=1+¢4/(a+4). Thenyn=y-“*4®, Now Poisson’s equation is 
1 d/,d¢ 
ae Pe 
r2 dr (+ =) = —4n0p, 

so that if we write z=7 x (87Ip, Bm)*?, 


@ 2 = a8 yet) 
dz dz ? 
Hence 
d dy ; 
aie 2 ale = 2 —(a+1/2) 
(a+) 3(2 2) aty-@+l/2), 


or, if x=2z(a+})?2, 


d 
y (x? dy/dx) = x2y-(@+1® 
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with the boundary conditions y=1 and y=o at x=0. Since a>} the index 
—(a+4)is negative. Equation (10.4) differs from Emden’s equation by having a 
negative index, and by having a positive quantity on the right-hand side of the 
equation. 

We select for special investigation the case a= § or C,=const x 5? exp { — 3y}. 
The central mass function is then &()=const exp {—3y}. At the centre the 
maximum number of stars per unit volume occurs at p= 4, i.e. m= 4m, and the 
maximum contribution to the light (assuming luminosity « (mass)**) is made by 
stars of mass m= $m. If then m is roughly equal to the solar mass the maximum 
contribution to the light is made by stars of mass 1-50 and absolute magnitude 
about + 3, which would give a spectral type of about F for the centre of the cluster. 
(Since the outer parts are richer in less massive stars the mean spectral type is 
later than this, and more detailed calculation might show that a cluster which has a 
colour corresponding to F or G must have a central mass function with a mean 
mass a little greater than the solar mass.) 

11. With C,=consty®*exp{—3u} and £(u)=const xpexp{—3yu} the 
differential equation to be solved is 





=(* =) =x*¥y8 (11.1) 
with y=1andy=oatx=o. The solution in series (useful when x <1) is 
Y eee 29 4 _1717 190138 
yurt ge 40° + 5040 1 088 640 + 399 168 000~ , 


Values of y from x = 0 to x= 1 were found from this series, and values of y for larger 
values of x found by quadrature. Results are described in Section 13. 

It is a great convenience of the form (10.1) that once y has been found a number 
of physically important quantities can be found without further numerical 
integration. We have, for a=$, ¢=3(y—1) and p=y-%, and in addition : 
relative number of stars per unit volume of mass m =m, 

¥4(*)/v,(0)=exp {mp} = exp {— 34 (y—1)}. 
Light emitted per unit volume (assuming luminosity oc (mass)**) proportional to 


| pw? * v(x) dp, ie. to y-™*, Number of stars per unit volume brighter than a 
“0 
given magnitude, i.e. for which m>pym, 


iw I+3umyy : 
| mai vAx)du oF (I+ a exp { — 3441 (y—1)}. 


12. As in the isothermal gas sphere, the mass of the cluster given by solution 
(11.1) is infinite. (The asymptotic solution of (11.1) is ¥/($).x", so that the 


8) 2) 
Mass 479 | r°y-3 dr is proportional to \, x2 dx, infinite.) Truncation of the 
“0 


Maxwellian velocity distribution, after the manner of Section 7, can be effected by 
using equation (9.6) instead of (9.5). We have, in fact, from (9.6) when the mass 
function is that of Section 11, or a=, 


p=const x |" uexp{—(3+0)u}.erfe (wy/(uin)} du 


cone PCS 9) ASH OATES H+ Hs + wa Gn, 
12.1 
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For each mass v,?— 2¢ is constant, so that mBv,2— 2m8¢ is constant, that is to 
say mw* + is a constant of the model, =, say. 
Accordingly 


p=const x (k—p)9?(3 +) 8k-7? (3 (3+)? +$(3+¥)(k—-Y)+(k—-Y)?}. (12.2) 
This reduces to p=const x (3+%)-* when ko, i.e. when the Maxwellian 
distribution extends to infinite velocity, as in the last section. Writing as before 
y=1+ 4u and setting c=1 + 4k, 

p=const x (¢— y)?® y~9 (359? + 28y (c—y) +8 (c—y)*} 
and since y=I, where p=po, 
~ y\3/2 24 rocy + 1sy2 
prmeeeve (=) m 8c? + I2cy + 15y 





c—I 8c? + 12¢4+15 


The differential equation to be solved, instead of (11.1), is therefore 


d dy c—y\32 8c? + 12cy + 15y” 
aw | BZ.) ee gty-8 2 
= (* ) ihe (=) 8c? + 12¢+ 15 nat 


with the boundary condition y=1 and y=o at x=o0. The solution of course 
terminates at y = c because here p = 0 by (12.3). 

Numerical quadrature of equation (12.4), with a particular value of c, presents no 
special difficulties, and once y has been found as a function of x the various 
quantities of physical significance can be found without further numerical inte- 
gration. 

We have at once the density p from (12.3) and the potential ¢ from ¥ = 3(y—1). 
The relative number of stars of a particular mass, per unit volume, is found from 





Ym (r) = const x 4 exp { <1 (3 + p)e} ° erf, \ ‘{u(k ss p)}, 
of which the practical form is 
Vm (7)/%m (0) = exp {— pp} erfy V/{u(k—) }/erfy V {uh}. (12.5) 
The number of stars per unit volume down to a given magnitude is found from 


(my?) = | _ Yq(7)dm=const x | wexp{—(3+)u}erfs v/{u(k—¥)} dy 


‘ium 
const 


— [ (2+ 3my)exp {—3nry} orf visu (e—y)}+8 (1-2) 


x (1442) fer Vise} 


ley 3/2 
+ 3 (1-2) "2 Guortexp (- ais) | 


The last two terms are not important in practice, and neglecting them we have the 
simple formula 


I+3¢1.y 
(I+3-1)¥” 


x erf, 1/ {3441 (¢ —y)}/erfev/ (341 (¢ — 1)}- 


n(my,1)/n(m,,0)= exp {— 341 (y—1)} 
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Lastly, if we take the luminosity to be proportional to the (mass)**, we can calculate 
the relative amount of light emitted per unit volume, which is proportional to 


| wm (7) des 
that is, to 


|, u®2exp {—(3+)u} erfy V/{u(k—-¥)} . dp. (12.7) 


This may be found from 


+. | wexp{—x} erfy {ay/x} dx=P (p+ 1) {1—Ta say (P +1 9)} 
ViJo 
4 V(p+8) a 


Va 2p+3 (1+a*)?*9?” 





B, (p, 9) 
B (,q) 


| x”-1(1—x)*1 dx, which is tabulated by Pearson.* 
0 


where /,(p, q)= and B,(p,q) is the incomplete Beta function 


ol 
Also B(p, q) is the complete Beta function | x?-1(1—x)*1! dx. Alternatively 
0 














-2 v7 2 (p+) tan—'a aii a 
|, xe“ erf, {ay/x} dx= 2p +3 {{. cos?” +1 6d6 — aa. 
In (12.7) we have p=$ and a?=(c—y)/y. Since 
net Ad =n 2-0 99. 9 erat 9272 conto +. 92.7:5 9-7-5-3 
| cos 6d6 = 5 8in0.cos 0 {os 0+ 3 008 6+ 3-6 0° 6+ .6.4°"* 8.6.4.2 
9-7-5:3 9 
= 8.6.4.2 10’ 
and since 
sin(tan-ta)=a/1/(1+a*), cos(tan-4a)=1/1/(1+4*)=y/c 
the expansion in y/c for the light / is 
_ yn), Ctam V(cly—1)+ Viv(e—y)j- P(V/e) 
bey * Tan! (e-1) + V(e=1) Pe)” ay 
where P(x)=1+ Gxt %.$x2+ §.§. S08 + 3.5.8.8 (x4— 1025). 
The total mass M of the cluster is given by 
M=4n7 (= 2). / (= r) 9/3 518, (12.9) 


Let L, be the emission of light per unit mass at the centre of the cluster. If the 
unit in both cases is the quantity for a star of mass m 


Ly= |” u8#exp {— 34} erfs V(hu) du ||” wtexp {—34)} erfe V(hu) du 
=3°521(4¥)/T (3), or 1-679, for large k. 


‘The total emission from the cluster is 


ro 3/2 
L=L,x 42 | ist dx | (= v7) Po". (12.10) 
0 





* Tables of the Incomplete Beta Function, ed. Pearson, London, 1934. 
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Accordingly the total mass divided by the total luminosity is 


I dy - 2 
[-L* (#2) / Ix? dx, (12.11) 


0 


the factor by which 1/L, is multiplied being a constant of the model, i.e. indepen- 
dent of py andj*. The degree of truncation relative to the velocity of escape at the 
centre of the cluster, in the notation of equation (7.5), is given by 


e2 kx, 


2 
“o 3 (= 2). +kx, 


but the degree of truncation relative to the Maxwell distribution depends on the 
mass, since v*,, = 2/(3j), so that 





’ 


(12.12) 


Uo?/v?,, = § pk. (12.13) 


13. Some calculations have been made with the formulae given in previous 
sections. These have not been made in any attempt to produce immediately results 
which can be compared with observation but in order to study from a mathe- 
matical point of view the effects of mass function and velocity function on those 
numbers which must ultimately be compared with observation. The reasons 
for taking up this rather detached position are that, in the author’s opinion, not 
enough is known about the mass function or the mass-luminosity relation in the 
objects concerned ; that we are quite at a loss to decide where to put the cut-off ; 
and that there are some rather disquieting indications that equipartition is not 
complete. The last remark is occasioned by the presence of cluster variables— 
presumably above average in mass—in the outer parts of the cluster w Centauri. 
For all these reasons the author feels that the comparison between theory and 
observation (at all events, of globular clusters) is not a matter to be made 
summarily, and it will not be attempted in this paper. 

Nevertheless, the calculations that were made do show how the truncation of 
velocity brings about a finite solution ; how the introduction of a mass function 
alters the density distribution quite significantly from that of an “isothermal gas 
sphere”; and how when there is a mass—and therefore a luminosity—function 
the various observable quantities are distributed differently with radius. On this 
last account it is quite wrong to compare an expected density distribution with an 
observed light distribution. The actual calculations carried out, or transcribed 
from the works of others in one case, are as follows : 


Designation Description Equation 


A Isothermal gas sphere for stars of only one mass. (6.2) 


B A truncated Maxwellian velocity distribution for stars of 
one mass only: k=7-82. (7-3) 


A full Maxwellian velocity distribution but with a central 
mass function proportional to » exp { — 3}. (10.4) 


Asin C, but with the velocities truncated: c= 10andk=27. (12.4) 


As D, but c=5 and k=12. (12.4) 
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The potential / is shown as a function of z (not x, which is z\/3) in Tables 1a 
and 1b. ‘The potentials do not differ, to the figures shown, between z=0 and 
z=1°5, in the five solutions shown together in Table 1a. Neither table is 
illustrated. 


TABLE I a TABLE I b 
Potential ys as a Potential ¢ as a function of radius = 
function of radius z, for the five solutions 
all five solutions 
z ib z A B Cc D E 
0°00 0000 2 0560 0'560 0566 0564 0°564 
0°25 o-olo 3 1°064 1:064 1090 1-090 1°085 
o'50 0'O41 4 1°571 1°57! 1650 1°646 1°635 
0°75 o-OgI 5 2°04 2°04 2°20 2°19 2:18 
1°00 o'158 6 2°47 2°47 2°73 2:72 2°68 
1°25 0°242 8 3°18 3°16 3°70 3°68 3°62 
1°50 0°338 10 3°74 3°72 4°58 4°54 4°43 
12 4°19 4°17 5°36 5°32 5°15 
16 4°89 4°84 6°83 6°7 6°37 
20 (541 5°34 7°96 7°87 7°36 
30 6-29 6°17 10°49 10°29 9:18 


4° 6°87 6°70 12°59 12°23 10°04 
50 7°31 7°05 14°43 13°86 11°28 
60 7°66 7°29 16:08 15°28 11°92 
70 =—-7'95 7°51 17°60 16°54 ae 
80 8-20 7°66 19°00 17°66 
90 8-41 7°78 20°3 18-68 
100 8-61 ey 21°6 19°59 


The reduced density p/p) is shown in Fig. 1 for all five models, the logarithm 
of the density being plotted against +/z. 





log” v ' T T T T T 
Pe 














iL. A. i nailte 1 i A. V 


i 2 3 4 5 6 7 
Fic. 1.—Log p/p» plotted against »/z for five cases. 
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To show how a number of important quantities are very differently distributed, : 
Fig. 2 has been prepared from the calculations for case D. The logarithms of 
the space densities (not projected !) of the numbers of stars of different masses, of 
the number of stars of greater than a given mass, and of the light emitted (assuming 
luminosity proportional to (mass)**), are shown plotted against 1/z. 

The distributions are quite different, so that it should be wrong, for example, to 
count all the stars in unit area on a plate, regardless of magnitude, and expect this 
count to behave like the density. 





." t 1 | 








1 1 1 l I 


oO | 2 3 4 5 9 


Fic. 2.—Distribution of various quantities in case D (unprojected). The logarithms of the 
relative numbers of stars t= 4, p=2, and u=1, of the relative number of stars brighter than those 
for which p,=2, and of the emission of light, all per unit volume, are plotted against ~/ z. 





The mass function for the object as a whole is not of course the same as the mass 
function for the centre. It has been estimated, for model D, by numerical 
integration. The result is shown in Table II. 


TABLE II 


Mass function for model D 
Logarithm of the relative numbers of stars of mass m=yum 


be 4 _2 I } } 4 
log number 5°37 3°38 0'00 1°51 3°03 3°55 


14. The overall properties of the finite cases B, D and E are as follows, in 
dimensionless units : 


Quantity , B D E 
Mass (2745) 86°6 1528 
Extreme radius 2, 94 256 


@ 
Total light J 27] dz 33'S 


0 
Central truncation v9?/u 9" ° o'82 


Total mass 
Total luminosity 





27 
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We can now show how the properties of case D, for example, in real units, 
depend on various parameters. We have 


2 = (Ar) = (821 poj?)*?r. 
If L, is the emission of light per unit volume at the centre, the total light emitted is 


L=L,x 4nx* | lz?dz. If the unit of light is that emitted by a star of mass m, 


Ly po= | w?texp{—3u}erty v(uk) dh / |” wexp {— 31} erfs (uk) 
for the central mass function chosen. This ratio tends to the limit 
3°? 1 (F)/T (3) =1-679 


Then L=4npyX*x 1-68 [ Is*ds. 
Jo 
Lastly the total mass is given by (7.4) or 


M = 4m (2%) po. 
Accordingly the total mass divided by the total luminosity is independent of the 
parameters A and pg, i.e. j? and pp, and is, for case D, equal to 27, the units being 
the mass and luminosity of a star of mass m. 


Half of the projected light in case D is emitted inside z= 8-4, so the effective 
radius of the object is 


r,=8 4x}, 
For m= ©, absolute magnitude supposed +5: for j?=(10km/sec)-*, and for 
Po = 1000 parsec”’, the dimensions of model D are as follows : 
Effective radius 25 parsec (proportional to py '?7~") 
Total mass 5 x 1070 (proportional to py -¥*7~*) 
Absolute magnitude —10™7 (emission proportional to 
Po “*7-*) 


For comparison, Shapley* gives a ‘‘ diameter”’ of 64 parsecs and an absolute 
magnitude of —11 for w Centauri. The fit means nothing, since two arbitrary 
constants are disposable. 

15. Despite the disclaimer, at the beginning of paragraph 13, of any intention 
of making a comparison between any of the models and any observations, we may 
allow ourselves to remark that the projected light in case D does fit de Vaucouleurs’ 
mean curve for elliptical nebulae remarkably well. This is shown in Fig. 3. 
No adjustment has been made, the quantity « in both cases being the projected 
distance from the centre divided by the projected effective radius, that is the radius 
within which half the light is emitted. This fit probably does not mean that the 
mass function and the cut-off are right (the cut-off is not, if the observations are to 
be relied on where the emission is very faint): but it may mean that the truncated 
velocity distribution is not a bad approximation and it certainly seems to mean that 
introducing a mass function is much better than supposing that all the stars have 
the same mass. This de Vaucouleurs has already suggested. 


* H, Shapley, Star Clusters, Appendix A. 
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Isothermal 








Fic. 3.—de Vaucouleurs’ empirical law log B= — A(«''*—1) plotted against reduced radius 
a, compared with points on the projection of case D (open circles). The projection of the isothermal 
and of Hubble’s model are also shown, these two with an arbitrary radius scale. 
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Summary 
When a flat photographic plate is bent on a spherical former, as in the 
Schmidt camera, strains are produced in the plate which must be allowed 
for in making measurements after exposure. In a previous paper one of 
the present writers obtained the necessary corrections appropriate to 4 
circular plate. In this paper the corrections for a square plate are obtained. 





Introduction.—In a previous paper* one of the present writers investigated, 
at the suggestion of Dr E. H. Linfoot, the deformation of a photographic plate 
in a Schmidt camera. In that paper only circular plates were considered. In 
this paper square plates are considered. When dealing with the circular plate 
certain assumptions were made, and similar assumptions are made where appli- 
cable in this paper. The modifications appropriate for a square plate are men- 
tioned in the section on boundary conditions. The solution for the square plate 
is, of course, more complicated and much of the simplification due to symmetry 
is lost. The problem is, however, shown to be equivalent in the mathematical 
sense to the problem of the bending of a clamped square plate by uniform 
transverse pressure. 






Y 


«Unstrained plate 


Fic. 1.—Diagrammatic representation of the strained and unstrained states of the plate. 


Displacements and strains.—The plate which is flat in the unstrained state is 
bent by means of a frame so that it is in contact at all points of its edges, with a 
spherical former. ‘The strained and unstrained positions are shown diagram- 
matically in Fig. 1. Rectangular axes Nx, Ny, Nz are taken such that Nx, Ny 

* W. M. Shepherd, M.N., 113, 450, 1953¢ 
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are in the plane of the unstrained plate and Nz is the inward normal at the point 
of contact N of the unstrained plate and the former. The displacements of the 
point (x, y, z) in the middle surface of the plate are u, v, w parallel to the axes 
Nx, Ny, Nz and the corresponding strains in the plane of the plate are s,, 5, @z,- 

The plate subtends only a small angle at the centre of curvature of the former 
so that the displacements are small and will be related to the strains by 


equations* (1). 
Ou 1 (dw\2 ) 
se" Ox 7 2 (5) 


_ dv, 1 (du)? 
sy = Oy + 2 dy ’ 
_ Ou, 2, dw du 
Pou By * Ox + Ox dy” | 
If small quantities of order higher than (x/R)* are neglected in the strains, 
the displacement w of a point S,(x, y,0) is given by the ordinate z of the point 


(x, y, 2) on the former. 
To the required degree of accuracy w is given by 





w = (x? + y)/2R, 


where R is the radius of curvature of the former. 
Equations (1) now take the form 


a he lh 
$e= 9 T2\R) ° 


ov I/y\* 
BAR) f 
Ou dv xy 
C gy = dy + Pa + R® J 
On eliminating the two displacements u and v from equations (3) the strains 
are found to satisfy the compatibility condition 





(4) 


The differential equation for the stress function.—Since the displacements are 
small the equations of equilibrium of an element of the plate in its own plane 
may be written in terms of the original position of the element. It follows that 
the mean stresses xx, yy, xy in the plane of the plate can be derived in terms of an 
Airy stress function x(x, y) from the body-stress equations, following the usual 
assumptions that x3, yz vanish on the faces of the plate, and the mean 23 across 
the plate vanishes. For reasons given in the previous paper 22 can here be 
neglected in comparison with the stresses in the plane of the plate, and unless the 
coefficient of friction is remarkably large, the surface shear stresses x2, yz will 
be likewise negligible. Here we may write 


ie. ae co 
m= 54, =a, w-- Zk. (5) 





* Cf. S. Timoshenko, Plates and Shells, pp. 304-5, 1940. 
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The strain components in (3) are independent of z so that the means of the 
stress-strain relations for the stresses in the plane of the plate are 


re ee a°x 
Es, =38— oF) = 54s — 05° 
ae ade e 
Es, = yy — 08% = 5% —o 4, > (6) 
a2 
Eo,,=21 + 0)8) = —2(1 +0) 5 











where E, o are Young’s modulus and Poisson’s ratio respectively. 


Equations (4) and (6) give immediately the differential equation for the 
stress function 


an we £ 

4 Br oa 

Wx= 7a +? Saas + FA =~ RE (7) 
The boundary conditions.—If it is assumed, as in the previous paper, that the 

edges of the plate are free from load acting in the plane of the plate, then the 

stress function y must, at all points of the boundary, satisfy the relations 

oe ey 

ast ~ Bnds ~°* (8) 


where ds and dn are respectively elements of the tangent and normal to the 
boundary. ‘These conditions are in this case equivalent to 


x= 2 =0 (9) 


on the boundary. 

The solution of the problem.—In the mathematical sense the problem is 
completely defined by equations (7) and (9). These equations for x are identical 
with the equations determining the transverse displacement in a plane plate 
clamped at its edges and subjected to a uniform transverse load qg per unit area 
where 

q E 
D =— R - (10) 


and D is the flexural rigidity of the plate. 

The solution of this problem has been obtained in a variety of ways. The 
most convenient for the present purpose is the numerical solution given by 
Timoshenko.* 

The displacements in terms of the stress function.—When the stress function x 
is known the strains are given by equations (6) and the displacements u and v 
are obtained by integrating equations re 


u= Js. dx — a +fi(y), ie 
v= | s,dy— Sa thle, 


Wy 


where f,(y) and /,(x) are functions of y and x only respectively. 


* Ibid., p. 222 ff. 
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Since u=o when x=0 for all values of y and v=o when y=0 for all values of x, 
it follows that equations (II) may be written 


: x 
4 = | sede 6R?’ 


y y® 
v= | svay— 6R?’ 


I * Ory 
u= =F), Dye dx 
I (¥d*x (12) 
v= E I. ax? rau 
Radial and cross-radial displacements. Correction of angles.—The polar 
coordinates (r,¢) in the plane of the unstrained plate of a point S, are related to 
the coordinates x and y by the relations x=rcos¢, y=rsing. In the strained 


y 


























Fic, 2. 


state the point S, moves to S and the displacements in the directions r and ¢ 
increasing are respectively u, and u, as shown in Fig. 2. These are related to 
the Cartesian displacements by the equations 


u,=ucos¢+vsing=(ux+vy)/r, | 
u, =v cosd—usind =(vx—uy)/r. 
The correction to the position angle ¢ is then given by 
5¢=u,/r. 
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The total correction to the off-axis angle A is 


ig 
ZS,0S =5A= > +R: , (15) 


Without radial stretching of the material of the plate the correction would be 
4,8, so that the correction due to stretching only is 
ae 4, 
x + R: 
The actual elongation of the length NS, of the mid-surface of the plate (the 
displacement u of the previous paper) is 

= Re=} RA? +u,. (16) 


The stress function and the displacements 
(a) The stress function.—F or a square plate of side a the stress function can be 


taken in the form 


Ee= 


X=X1t X2t Xs (17) 
where 
(—1)"-D2 max 
ee Cc —_—— 
m 
(m7 /2) tanh (mz/2)+2 my (mry/a).. 
osh —— — —-—~=——_ s 
2 cosh (mz/2) a 2 cosh (m7/2) 





mry | 
h af’ 


( i rj 1/2 mx E,, 
sabes als co 
m a cosh(m7/2) 


sh m7 _ TY ml 
a a a 
( am r)m—D2 mry En 
x 8S 7 cosh (mr/2) 


mir 


mir e 
x <— tanh — cosh —— — —— sinh — 
2 2 a a 


max mnx et 
, 
a 


where the numerical values of the first four coefficients E,, are 
E, =0°3722; E,= —0°0380; E;= —0-0178; E,= —0-0085. 


(b) The displacements.—As the displacements are symmetrical with respect to 
the line d = }1, so that v(x, y)=u(y, x), we need only evaluate the displacement u 
for a quadrant of the plate. 

On introducing the stress function x given by equation (17) into the first of 
equations (12) it is found, after some manipulation, that 


Co 


cosh (mzry/a) y sinh (mry/a) 
amt + Pm cosh (m7/2) — Qu Cosh a 
R sinh (m7x/a) x cosh (m7x/a) x 
{ ™ cosh(mn/2)  ~"a Shi ~ 6R?’ 


where 
K =4a3/R*z'*, A,,=(—1)"""2, 
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and 


mr o E.. 
Pn = UI +o) tanh — [m+ E,,] + a om 


7mm 
QO, =(1+0)[m*+E,,], 

mx I-cE£,, 
R,,=}(1 +0)E,, tanh — + mr pom 
Su=H1+0)Eq. 


By arguments similar to those used by Timoshenko it may be concluded that the 
errors due to terminating these summations after the first four terms will be very 
small. 

Numerical results for the glass plate (o=0-25).—The total correction of the 
off-axis angle A (i.e. the quantity to be added to the angle calculated from the 
measurements of the plate after release in order to obtain the angle A at the time of 
exposure) is given by A= K, y*, where 2Ry =a, and 6A is measured in radians. 
The coefficient K, is given in Table I for the fourteen points shown in Fig. 3. 














Fic. 3. 


@ Points for which the displacements are given in Tables IV and V. 
© Corresponding points in the strained state (displacements much magnified). 


The correction of the off-axis angle due only to stretching of the plate is given 
by «= K,y’, where K, is given by Table II. 

The correction of the position angle is given similarly by 54 = Ky? and K; is 
given by Table III. K,is zero when ¢=0° or 45°. 

In each of Table I, II, III the maximum possible error is + 0-0005. 
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The corrections in seconds of arc appropriate to a square plate which circum- 
scribes a circle which forms a spherical cap of angular diameter 6° are given in 
Tables IV and V. 


TABLE I 
Values of K, 





j 





0°1260 0°3714 
0°1293 0°3742 
0°1354 0°3855 
01389 0°3948 


























Tasce II 
Values of Ky 








0°0381 
0°0409 
00522 
00615 


























TABLE III 
Values of Ks 





























TABLE IV 


Values of the correction (€) to the off-axis angle A due to stretching only 


eg 1° > 3 28’ 


°° 0”-49 

15 o”'50 

3° 0°53 

450 0""54 
Circular plate 0” 46 





























The last line in the table gives the values obtained in the previous paper for a 
circular plate of angular diameter 6°. 
In this table the maximum possible error is + 0"-02. 
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TABLE V 


Values of the correction 5 to the position angle 4 


= 


15 
30 


























In this table the maximum possible error is + 0"*3. 
If values of the corrections are required for intermediate values of A and ¢ they 


can be found with sufficient accuracy by graphical interpolation using the values 


given in the tables. 
Fig. 3 is not accurate but only serves to indicate the nature of the distortion 


represented by the values given in Tables IV and V. 


Bristol : 
1954 February 19. 





MAGNITUDES AND COLOURS OF 123 STARS IN THE 
NEIGHBOURHOOD OF THE SUN 


G. G. Yates 


(Communicated by the Director of the Observatories, Cambridge) 


(Received 1954 March 24) 


Summary 
Magnitudes and colours of 123 nearby stars have been determined photo- 
electrically. ‘The material tends to support the view that the main sequence 
in our neighbourhood is narrow. It appears that errors in trigonometrical 
parallax are too great to permit the detection of fine structure in the diagram. 
There is some evidence of the existence of a sub-dwarf sequence lying about 
1™-5 below the main sequence. 





The magnitudes and colours of 123 nearby stars have been obtained, using the 
pulse-counting photoelectric photometer attached to the 15-inch Huggins 
Refractor at the Cambridge Observatories. The observing list was compiled 
from the Yale Parallax Catalogue (1936) following these rules : ‘The stars were to 
be (a) north of declination + 15°, (6) within 20 parsecs of the Sun, and (c) with 


visual magnitude between 5™-0 and I0™-5, 

Criterion (a) ensured that all the stars could be compared directly with N.P.S. 
stars at equal altitudes, (6) is necessary if even moderately well-determined 
absolute magnitudes are to be secured, in view of the relatively large errors of 
trigonometrical parallaxes. 

The magnitude range was chosen to suit the photometer. For stars brighter 
than 7™-0 a seven-inch diaphragm was placed in front of the O.G. At the other 
end of the range a red star of visual magnitude 10™-5 has a photographic magnitude 
near 12™-0, when the star-light intensity becomes about equal to the background 
intensity (sky plus noise), "This was regarded as the useful working limit for the 
instrument. 

Before the completion of the programme the ‘‘ General Catalogue of Trigono- 
metric Stellar Parallaxes’’, Yale, 1952, came to hand. The parallaxes quoted in 
‘Table II and used for the calculation of My are taken from this catalogue. It will 
be noticed that these more recent determinations of parallax have placed a few of 
the stars beyond the 20 parsec limit. 

‘The measurements were made in two colours, using the same technique as for 
the photoelectric part of an earlier programme (1). Blue magnitudes, not far 
from the international pg system, were obtained with a Wratten 39 filter; yellow 
magnitudes, near the international pv system, with a Wratten 5 filter. Fig. 1 was 
produced by combining the published spectral response curve of the 931-A 
photomultiplier with the transmission curves for the two filters, and shows the 
two spectral regions covered by the measures. ‘The colour reductions do not, 
however, depend on these curves, but were taken from the same two non-linear 
correction curves as were used in the work just mentioned. ‘These were obtained 
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directly from measurement of N.P.S. stars, and reduce the observations to the 
international pg and pv systems. No attempt was made to vary the corrections 
from night to night. All observations were made at about the same altitude, and 
in that case our experience is that night-to-night changes of the colour corrections 
can be neglected in comparison with other sky irregularities. ‘There is a fuller 
discussion in the paper quoted. 

Three comparison stars were used: (i) BD 88° r12 (N.P.S. 2r), which was 
used as a standard for the brighter stars (m<7™-o), (ii) BD 88° r14 (N.P.S. 4r) 
and (iii) BD 88° 9 (N.P.S. 8). A bigger selection of comparison stars would have 
been desirable, but would have consumed more observing time. Experience 
has shown that weather conditions in Cambridge make progress in such a photo- 
metric programme as this very slow. ‘The comparison stars had been well 
observed during the previous programme and are believed to be constant in bright- 
ness. ‘This earlier work also produced a good estimate of both blue and yellow 
accidental errors. ‘The standard error for a single observation was there found to 
be 0™-06 in the blue and 0™-04 in the yellow. These values have been confirmed 
in the present work where the errors are 0™- 054 | for blue, 0™-041 for yellow and 
o™-046 in colour index. 


lee | 


Relative 100 T 
Sensitivity 











0 
3000 4900 5900 6000 7000 A 


Fic. 1.—Approximate spectral regions covered by the blue and yellow measures. 





One complete observation of a star involved a sequence of measurements in 
each colour as follows: sky near comparison star, comparison star, sky near field 
star, field star, sky near comparison star, comparison star, twelve readings in all, 
each lasting for one minute. The time taken for one cycle of this kind was usually 
20-25 minutes; it was never less than 18 minutes and seldom exceeded 30 minutes. 
Since the transparency of the sky can change considerably during such an interval 
it seems likely that most of the error comes from this source. 

In view of these errors one would have preferred to observe in each colour 
about ten times, but this was impracticable, and in fact the total number of 
observations made on the 123 stars in each colour was about 500. ‘The standard 
errors of the mean blue and yellow magnitudes and of the colour index are then 
respectively 0™-027, 0™-o21 and 0o™-023. ‘The standard errors of the means for 
the separate stars are given in Table II. These errors in apparent magnitude, 


15* 
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although considerable, are of secondary importance in the determination of the 
absolute magnitude, the errors in the latter being controlled almost entirely by 
the errors in parallax. ‘The standard error in the distance modulus M — m for all 
the stars is o™-38, so that the error in absolute magnitude M contributed by the 
uncertainty in parallax is more than ten times the error due to uncertainty in 
apparent magnitude. In only a few stars with very well-determined parallaxes is 
the parallax error small enough for the total error in absolute magnitude to be 
influenced appreciably by the error in apparent magnitude. ‘The standard error 
of M—~m was derived as follows: 

Suppose the parallax is p and its catalogue probable error Ap. ‘The standard 
error in p will be 1-483 Ap. Also 


=5+2-17Inp, 


whence 5(M—m)=2172. 


0-10 


0-05 ~~ a e 
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Fic. 2.—Error cf one observation. 





(a) against photographic magnitude. 
(6) against colour index. 


(c) against declination. 


Writing 5p = 1-483Ap the standard error in M—#m due to a probable error Ap in 
parallax is 
oy a 

It might, perhaps, be anticipated that the errors in the measures of individual 
stars would show a dependence on magnitude or colour or declination (the latter 
measuring the separation between the star and the polar comparison star), but any 
correlation of this kind appears to be slight. Fig. 2 shows the dependence of errors 
of one observation on these three parameters. Each point in the diagrams is the 
mean from ten stars. ‘The diagrams indicate somewhat greater errors amongst 
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the fainter stars and possibly also amongst the redder stars. ‘These two ettects are 
not unrelated, since amongst a group of stars within a definite distance the fainter 
stars will on the whole be the redder. ‘There appears to be no correlation between 
error and declination. 

As it is not practicable in Cambridge to restrict observing to first-class nights a 
system of weighting based on the variation of the comparison star during the time- 
interval of each observation was adopted. If the variation was less than 2 per cent a 
weight 3 was assigned ; if the variation lay between 2 per cent and 6 per cent the 
weight was 2; if the variation lay between 6 per cent and 10 per cent the weight 
wasI. Inthe few cases where the variation exceeded 10 per cent the observation 
was rejected, as were the few observations made on thick nights when the intensity 
of the comparison star fell to less than half its normal value. It is not claimed 
that this system of weighting is ideal, as it samples the sky variations only at the 
pole, but it has the merit of being independent of the subjective judgment of the 
observer, who might otherwise be tempted to underweight any observation show- 
ing a large deviation from the mean. 


TABLE I 





Deviation 6 from mean line | Observed number} Expected number 





° 6<1:0S.D. in M 
o S.D.<8<1°5 S.D. 

+5 S.D.<ib< ° S.D. 

0 S.D.<—6< 3:0 S.D. 

°o S.D.<56 


own ze 

















A small correction was applied for differential extinction between the pole and 
the field, with a mean extinction coefficient 0™-5 in the blue and o™-3 in the 
yellow. Since the stars were always observed at or near the zenith distance of the 
pole the correction rarely exceeded 0™-02 in either direction in the blue, or 
o™orinthe yellow. ‘The maximum correction applied (once only) in the blue was 
o™-08 and the mean correction for all stars taken without regard to sign was 0™-012; 
the corresponding values for the yellow are 0™-05 and 0™-007. ‘These figures 
show that quite a large error in the assumed extinction coefficient would not be 
likely to change the final mean magnitudes significantly. Owing to the variability 
of the Cambridge sky it is impracticable to measure the actual zenith extinction 
whenever observations are made. 

Discussion of results.—Fig. 3 is a plot of absolute visual magnitudes against 
colour index. From the diagram it appears that the majority of stars of colour 
index less than about 1-22 fall near a straight line. As the colour index approaches 
1-4 the line bends over and finally runs almost parallel to the luminosity axis. 

Most of the scatter in the diagram is due to uncertainty in parallax, i.e. 
uncertainty in M. ‘This may also be demonstrated, at any rate for the linear 
portion of the sequence, by comparing the deviation of the individual stars from the 
mean line with those predicted from the parallax errors. Omitting the known 
binaries, there are 86 stars with colour index less than 1-22, giving a distribution of 
deviations in absolute magnitude as in ‘Table I. 

‘rom this analysis it appears that apart from seven stars whose actual deviation 
exceeds 3-0 times the S.D. in M the scatter in the diagram may be attributed 
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wholly to uncertainty in parallax. ‘The material tends to support the view now 
commonly held, that the main sequence for stars in our neighbourhood is quite 
narrow; it also suggests that the probable errors assigned to trigonometrical 
parallaxes are fairly accurate. Clearly, however, there is at present no hope 
whatever of determining a fine structure of the H—R diagram through trigono- 
metrical parallaxes. Whether such a fine structure really exists is still quite 
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Fic. 3.—Colour-luminosity diagram for 123 stars near the Sun. 


| e| Binary stars. @ Possible sub-dwarfs. 





uncertain, for redeterminations of the H—R diagrams of the Pleiades by Johnson 
and Morgan (2) and of the Coma Cluster by Thiessen (3) have failed to confirm the 
results obtained by Eggen. 

‘The seven stars which appear to stand off significantly from the mean line are 
numbers 7, 9, 18, 26, 62, 71, 102 of Table II. Of these seven stars, 7, 9, 62, 71 and 
102 lie below the main sequence. ‘The evidence here is inadequate, but possibly 
these stars form a part of a sub-dwarf sequence, lying about 1™-5 below the main 
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sequence. Star 62 (H.D. 103095) is described in the 1952 Parallax Catalogue as a 
sub-dwarf. Numbers g and 62 are also included in Eggen’s observations and are 
described by him as sub-dwarfs(4). Star 95 which is situated well below the main 
sequence (= 2-3 S.D. in M) and the binary 96 would also appear to belong to 
this sub-dwarf sequence. The two remaining stars of the seven singled out 
above, numbers 18 and 26, stand above the main sequence. The position of 26 
rests on one observation only, so that both colour and absolute magnitude may be 
seriously in error. Star 18 seems to be of higher than average luminosity; its 
position is relatively well determined, the S.D. in M being less than average, 
whilst the colour rests on four observations, the S.D. of the mean being 0™ 012. 

Ten of the stars included in the list are binaries (or more complex systems). 
Eight of these—15, 17, 47, 58, 68, 73, 83 and g7—stand above the main sequence as 
might be expected. ‘The absolute magnitudes of the components are greater than 
the absolute magnitude of the combined system, and hence the individual com- 
ponents will deviate less from the main sequence than the system asa whole. This 
receives rough confirmation when estimates of the absolute magnitudes of the 
components, based on their relative magnitudes as given in the 1952 Parallax 
Catalogue, are made. ‘The two remaining binaries are numbers 50 and 96. We 
have already suggested tentatively that 96 should be classed with the sub-dwarfs. 
The apparent magnitudes of the components are 9-7 and 10-1, and the spectral 
type is K4, according to the Parallax Catalogue, the absolute visual magnitudes are 
thus 8-2 and 8-6. 

Number 50 has been included in the list of binaries although it is not certain 
that it is a double star. This star is BD + 20° 2465 (=Ci 1244) and is included 
amongst Eggen’s (4) stars (his number 33). ‘The star has been discussed by 
Reuyl (5) who, from an analysis of the variations in proper motion over the period 
1915~1942, deduced an orbit and estimated the mass of the dark companion to be 
about 0-1 of the mass of the primary. More recent proper motion measures by 
van de Kamp and Miss Lippincott (6) fail to confirm Reuyl’s elements but still 
indicate the possibility of orbital motion. On the basis of this evidence it would 
appear that the companion, if any, would be too faint to affect the photometric 
measures. Our photometric measures would place the star on the main sequence, 
in agreement with Eggen, but contrary to Struve’s statement (7): ‘‘ A star which 
appears to be noticeably above the main sequence is BD +20° 2465, whose 
spectral type is M5 and whose absolute magnitude is 10-9.” 

The relation between colour index and spectral type is shown in Fig. 4, whilst 
‘Table III lists all the types included in the programme with their mean colour 
index and also the number of stars of each type. The binaries are not included in 
this table. 

Comparison with other work.—Recent photoelectric measures of some of the 
stars in our list have been published by Johnson and Morgan (8), and by Eggen (4). 
We have measured 15 stars in common with Johnson and Morgan and 21 in 
common with Eggen; 13 starsare common tothe three lists. ‘The blue and yellow 
magnitudes given by the various observers are collected in Table IV. Johnson 
and Morgan’s magnitudes have been reduced to our system by means of the 
equations : 


By — By =0-212 Cy —0°372, 


Vy — Y; =0-086 Cy — 0-074, 
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TABLE II 
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Serial number 
1952 Parallax catalogue 
number 


Spectral type 
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| Total weight 
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TABLE II (cont.) 
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The errors in parallax are probable errors : all other errors are standard errors. 
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where By, Yy, Cy are our blue and yellow magnitudes and colour index; B;, Y; 
are Johnson and Morgan’s B and V respectively. Eggen’s magnitudes were 
reduced using the equations : 

By — Bg =0°252 Cy —0°205, 

Yy — Yy=0-032 Cy +0018, 
where By is Eggen’s Pg,, and Yy is equal to (Pg, — C,) of Eggen’s work. 

Over the colour range concerned in this comparison (from about C.I. + 0°30 to 

+ 1:50) the relation between the three magnitude systems is adequately represented 
by linear functions of colour index, which have been obtained by least squares. 
The existence of these colour equations calls for comment. Both Eggen and 
ourselves had supposed that our magnitudes were on the N.P.S. system, so that the 
systematic differences represented by the equations ought not to have existed. 
Similarly Johnson and Morgan’s V magnitudes (yellow magnitudes) are claimed 
to be very close to the I,, system; but here also there are systematic differences 
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Spectral Type 


Fic. 4.—The larger circles are mean points derived from three or more 
stars of the same spectral type. 


both between Johnson and Morgan’s yellow magnitudes and ours; and between 
Johnson and Morgan’s yellow magnitudes and Eggen’s. It appears that the 
problem of defining a unique magnitude system is not yet satisfactorily solved. 

‘Table 1V shows that when the three sets of measures have been reduced to a 
common system the agreement between the magnitudes, in both colours, assigned 
to the individual stars by the different observers is satisfactory. ‘The deviations of 
the different measures from the means of all are generally small. 

Acknowledgments.—It is a pleasure to acknowledge indebtedness to Professor 
R. O. Redman who suggested this programme and assisted with the observational 
routine. ‘Thanks are also due to other members of the Observatory Staff who 
assisted with the observational work, in particular, Dr Dewhirst, Mr Ovenden, 
Mr Hignell and Mr Wallman. 


The Observatories, 
Cambridge : 
1954 March 22. 





Magnitudes and colours of 123 stars near the Sun Vol. 114. 


References 


(1) A. Beer, R. O. Redman and G. G. Yates, Pg and Pv Magnitudes of 7-10™ stars in 
the +-15° Selected Areas. Mem. R.A.S., 67, Pt. 1, 1954. 

(2) H. L. Johnson and W. W. Morgan, Ap. F., 114, 522, 1951. 

(3) G. Thiessen, Z. f. Ap., 32, 59, 1953. 

(4) O. J. Eggen, Ap. F., 112, 141, 1950. The two stars are numbers 7 and 41 in Eggen’s list. 

(5) D. Reuyl, Ap. F., 97, 186, 1943. 

(6) P. van de Kamp and S. L. Lippincott, Ap. 7., 55, 16, 1949. 

(7) O. Struve, Stellar Evolution, p. 36 ff. 

(8) H. L. Johnson and W. W. Morgan, Ap. 7., 117, 313, 1953. 





METEOR STREAM-WIDTHS AND RADIANT DEVIATIONS* 


Fred L. Whipple and Frances W. Wright 


(Communicated by the Director of the Harvard College Observatory) 
(Received 1954 April 1) 


Summary 
Among eight showers of photographic meteors, the average deviation 
from the moving mean radiant increases nearly linearly with the spatial 
stream-width. 





We have assembled the measured average deviations of individual photographic 
meteor trains from the moving mean radiants for seven major meteor showers 
photographed on the Harvard Meteor Programme and for one photographed by 
P. M. Millman. 

To derive the width of a stream in space at the chord of the Earth’s crossing, we 
have reduced the observed photographic duration of the shower by the sine of the 
angle (y) between the vectors of the Earth and the stream in space. Let Vy be the 
heliocentric velocity of a meteor stream, V’, its geocentric velocity, neglecting the 
Earth’s attraction, and A the elongation of the corrected radiant from the apex of 
the Earth’s motion. ‘Then the angle y is given by 

siny =V,/Vysina. (1) 

Fig. 1 shows the results graphically for the seven photographic showers: 
Draconids (1), Leonids (2), Geminids (3), Orionids (4), Southern and Northern 
Taurids (5), Perseids (6) and for the 6 Aquarids(7). The stream-width in days is 
the abscissa and the average deviation in degrees (the |O-C| residual in the 
Geminid, Leonid, Perseid and Taurid papers) is the ordinate. A least-squares 
solution gives the straight line. ‘The 6-Aquarid point has been added later and is 
not included in the solution. 

Table I lists pertinent average data for the eight photographic showers. ‘The 
quantities A, V,, Vy and sin y have been defined above. Dur. represents the 
observed photographic duration in days, Str.-Wd. the spatial stream-width in 
days and Av. |O-C | the average deviation, in degrees, of an extended trail from 
the moving mean radiant. The last column lists the number of trial residuals 
upon which the deviation depends. We have listed the showers in Table I in 
order of increasing stream-width. 

The greatest uncertainties in the data of Table I and Fig. 1 lie in the durations 
of the showers. ‘The adopted values are based upon the earliest and latest dates 
of identifiable photographic members and may differ considerably from durations 
obtained by visual or radio techniques. The duration of the 6 Aquarids is quite 
uncertain, as indicated by the horizontal dotted line in Fig. 1. 

‘lhe observed correlation between stream-width and radiant scatter proves that 
the spatial motions of individual meteoroids differ increasingly from the average 
motions of other meteoroids near them in the stream as the width of the stream 


* This research has been carried out under a contract with the U.S. Office of Naval Research, 
Contract No. N5ori-07647. 
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Fic. 1.—Scatter vs. stream-width. 


increases. Hence, the stream must derive much of its cross-sectional area from 
some random type of disturbance to the motions of individual meteoroids. 
A detailed explanation of the correlation will involve a rather complex study of 
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possible forces which affect the motions of meteoroids and will not be attempted 
here. ‘lhe fact that the T'aurid and Geminid streams follow the scatter stream- 
width relationship shown by such recognized cometary streams as the Draconids, 
Orionids and Perseids suggests a common type of disturbance which affects them 
all. This similarity perhaps casts doubt on Hoffmeister’s proposal (8) that the 
Taurids and Geminids are “ ecliptical currents ”, different in character from the 
cometary streams. 

A comparison of the present measures of radiant scatter with those from 
fainter meteors, either visual, radio or photographic, by the Super-Schmidt 
cameras, should provide valuable information concerning the development of 
meteor streams. An increased scatter for fainter, and presumably smaller, 
meteoroids would prove that physical forces are significant in spreading the streams. 
These forces might act during the separation from the parent body or after the 
separation. A scatter essentially independent of brightness, on the other hand, 
would suggest that the perturbative action of the planets is the major factor in 
dispersing the streams. 


Harvard College Observatory, 
Cambridge, Mass., U.S.A. : 
1954 March 26. 
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ON THE OCCURRENCE OF COMMENSURABLE MEAN MOTIONS 
IN THE SOLAR SYSTEM 


A. E. Roy and M. W. Ovenden 
(Communicated by the Director, University Observatory, Glasgow) 
(Received 1954 April 6) 


Summary 


It is established that, in the solar system, among the planetary and satellite 
systems, the occurrence of commensurability between pairs of mean motions 
is more frequent than in a chance distribution. It is suggested that the 
observed distribution may be due either to a property of the mechanism of 
formation of the solar system, or to an inherent stability of commensurable 


configurations. 





1. Introduction.—Of recent years, a number of theories of the origin of the 
solar system have been developed, notably those of O. J. Schmidt (1, 2), G. P. 
Kuiper (3) and C. F. von Weizsiacker (4, 5), which suggest that the distribution of 
planetary and satellite orbits is not due to chance. ‘These theories, in turn, have 
re-directed attention to the many empirical relationships which have been 
discovered between the distances of the planets and satellites from their primaries, 
of which the best-known but by no means the most accurate is the Bode—Titius 
law relating the semi-major axes of the planetary orbits. It is less well known that 
a law similar in form to the Bode—Titius law is applicable to Saturn’s system of 
satellites (6). 

These and other similar relationships have often been considered to be 
entirely fortuitous. However, in 1913, Miss M. A. Blagg (7) proposed a more 
complicated relationship which is applicable to the planetary orbits and to the 
satellite systems of Jupiter, Saturn and Uranus. One of the present authors (8) 
has recently drawn attention to the fact that all the bodies discovered in these 
systems since the publication of Miss Blagg’s paper (viz. Pluto; Jupiter IX, X, XI 
and XII; Uranus V) fit this relationship to a high degree of accuracy. (‘There 
are insufficient satellites at present known in the Neptunian system to apply Miss 
Blagg’s formula.) 

The above considerations have led the present authors to a re-examination of 
the problem of the distribution of planetary and satellite orbits from a rather 
different viewpoint. It has long beeg noted that certain pairs of bodies have 
mean motions which are very nearly commensurable. Outstanding examples are 
Jupiter and Saturn, whose mean motions are approximately in the ratio of 5 to 2, 
‘Tethys and Mimas and Dione and Enceladus where the ratio is I to 2, and Hyperion 
and Titan where the ratio is 3 to 4. In the Jovian system, not only are the mean 
motions of Europa and Io and Ganymede and Europa very nearly in the ratio I to 2, 
but they are also connected by the well-known Laplace relation 


Ny — 3g + 2N3 = 0, 


where n,, 7, and n, are the mean motions of Io, Europaand Ganymede respectively. 
‘This case will be discussed later. 
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The present paper gives an account of an analysis of the mean motions of the 
planets and satellites made with a view to establishing whether or not the number 
of such near-commensurabilities could be due simply to chance. 

2. The occurrence of commensurabilities in a chance distribution of numbers.— 
A pair of numbers n, and n, can be said to be approximately commensurable if 


tp Mp 
™, = A ’ (2.01) 
where n, <n, and A, and A, are integers with A, < A,. It is clearly necessary in 
practice to set an upper limit to the sizes of 4, and Ag, since for any pair of numbers 
N/n,, a pair of integers A,/A, can always be found to satisfy (2.01) to any required 
degree of accuracy if there is no such limit. It is also necessary to specify the 
accuracy with which (2.01) must hold to be considered significant. If € is a small 
positive quantity defined by 
Ne A 

e= = - z : (2.02) 
then ¢ is a measure of the degree of commensurability of n, and 7, (m,, mg, A, and A, 
being positive). 

Quite arbitrarily, the upper limit to 4, and A, was taken to be 7. Within this 
limit, with A,<A,, there are 17 fractions between } and 1 (excluding unity), 
symmetrically disposed about 4. The smallest difference between adjacent 
fractions is 4=0-:02381. If values of « <¢g<0-01Ig0 in a chance distribution 
of numbers are considered, then there is no ambiguity in the choice of integers A, 
and A,. The total extent of the regions in which a ratio m,/n, must lie to be 
commensurable within ¢, in the range } to 1 is 34€). Thus the probability of any 
one ratio ”,/n, in a chance distribution lying within ¢, of one fraction A,/A, is 


17 X 2€) x Z = 39°67 €,. 


Limiting values of n./n, to the range + to 1, the number of ratios m,/n, lying within 
the prescribed limit €, of a fraction A,/A, in a chance distribution of v pairs is 7N,, 
where 

?N,, = 39°67 egv. (2.03) 

3. Method of analysis.—Because of the adopted condition that A, < A, <7, 
only pairs of bodies (planets or satellites) whose mean motions do not differ by 
more than a factor of 7 are considered. The mean motions have been taken from 
Connaissance des Temps for 1949, with the exception of the mean motions of 
Pluto and Miranda (Uranus V) derived from data in The Handbook of the British 
Astronomical Association for 1953. For every pair of bodies in the solar system 
falling within the above criterion, the ratio of the mean motions was found, and the 
deviation « from the nearest fraction A,/A, computed. The results of these 
computations are given in Tables I to IV for the total of 46 pairs. The numbers of 
pairs for which e <«, were then found for various values of ¢, < 0-011 go, and com- 
pared with the numbers predicted by (2.03) for a chance distribution. 

A similar analysis was then performed for the semi-major axes of the orbits of 
the planets and satellites. ‘This analysis was limited to the same sample of 46 
pairs of bodies. ‘The range of ratios of semi-major axes was thus (7)~?* to 1, or 
02733 tor. Within this range, there are 13 fractions (excluding unity) with both 

16 
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Taste [I 


Saturnian System 
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208 616 +0:008 616 
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‘212 322 +0°o12 322 
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Jovian and Martian Systems 
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numerator and denominator <7. Hence for the semi-major axes the number 
7, to be expected in a chance distribution of v numbers is given by 
7N,, =13 X 2€9 x (0°7267) 1 
= 35°79 €ov- (3-01) 

As a check on the validity of formula (2.03), an analysis, identical in method. 
with that used for the mean motions in the solar system, was made of a control 
distribution of numbers taken at random from a table of seven-figure logarithms 
and arranged in ascending order to form various systems of mean motions of bodies 
in a “fictitious solar system ’”’. 

In Table V, the numbers of pairs of bodies (in a sample of 46 pairs) with « <«, 
for various values of €, in the actual solar system and in the control distribution are 
compared with the predictions of (2.03), and the number in the distribution of 
semi-major axes of the solar system is compared with the predictions of (3.01). 
This table shows that, while the numbers for each value of «, for the control 
distribution and for the distribution of semi-major axes agree well with the 
predictions from a chance distribution, the distribution of mean motions in the 
actual solar system shows a marked excess of observed over predicted numbers. 


TABLE V 


€o 00089 0°0030 
'N,, : 16:2 5°3 
Observed no. control 
distribution 13 3 
Observed no. mean 
motions 26 12 
o. * 5 14°6 
Observed no. semi- 
major axes 15 


4. The assessment of multiple counts.—In one respect, the numbers given in 
Table V might appear misleading. If « represents the number of a planet or 
satellite, considering the body nearest its primary to be «=1, then, if « is com- 
mensurable with («+ 1) to within ¢, and (« + I) is commensurable with (« + 2) to 
within €9, the commensurability of « with (« + 2) to within €) may or may not be an 
independent commensurability. 

Formally, if 

N(x,41) ~e 
“rales 


€1 (4.01) 


a 
and 
Ny +2) ) 


+€ -02 
My +1) ’ (4 ) 
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|<2 | <€ 9, &9 being a small positive quantity), then 


Mare _ [ & g 
“a2 ~(F+5) (+4) 


eg g e 
~ th + (fat ‘er) > 





238 A. E. Roy and M. W. Ovenden, On the occurrence Vol. 114 


, ‘ a eg. ; ; 
It is therefore necessary to calculate the probability that B isa vulgar fraction with 


eg <fh<7, while oc + “6 <e,. Considering first the chance distribution, 
equation (2.03) shows that the probability p that a pair of numbers m, , ,, 7, give a 
ratio satisfying (4.01) is p=39°67«€). Hence the probability that two pairs of 
numbers m,,. 3), #,, and m,,,.)) M%q.4.1) Will occur satisfying both (4.01) and (4.02) is p?. 
The number of ways in which the 17 vulgar fractions lying between } and 1 
(excluding unity) can be multiplied together in pairs to give vulgar fractions lying 
within the same range is 99, but of these, only 31 have both numerator and 
denominator not exceeding 7. Therefore, ina chance distribution, the probability 
of there being three mean motions related by (4.01), (4.02) and (4.03) with 3 
being a vulgar fraction and eg < fh <7 becomes # p?. 
This probability must now be multiplied by the probability that in these cases 


é 
Feat t*? Ep. (4.04) 


A precise evaluation of this probability in an actual case would be tedious. How- 
ever, it is clear that the probability cannot exceed unity. Furthermore, in the 
simplified case of the relation 
| my, + mM, | < €o 
a eT &g 

where |m,|<«, and |m,|<e,, the probability is ?. Since 7 <1 and i <1, the 
probability in the actual case must lie between } and 1. A value of # has therefore 
been adopted. 

The final probability, P, that in a chance distribution there will be three 
numbers satisfying conditions (4.01) to (4.04) is thus 


—ty 31 2 
P==x%xp 
4 
0 


= X 3p X (39°67) 9? 
= 394°2 €9”. 
P=0-056. 


In the control distribution there are 15 pairs of type [~, («+ 2)], and the expected 
number of commensurabilities arising according to the above argument is 
15p=0-84. ‘The observed absence of any such commensurabilities in the control 
distribution is therefore in good agreement with the prediction. 

As can be seen from Tables I to IV, the total of 46 pairs of bodies in the solar 
system i$ made up as follows :— 


For €g9 = O'OIIO, 


24 pairs of type a, (+1) 

17 pairs of type a, (~+2) 
4 pairs of type a, (x +3) 
I pair of type a, (x+4). 


For the distribution of semi-major axes, P=? x 2 x (35°79)" 9? =0°043 for 
€9=0°011g. Hence the expected number of “spurious” counts is 0-73. ‘The 
observed number is 1; again, the agreement is good, as is to be expected in view of 
the agreement (in Table V) of the control distribution and the distribution of 
semi-major axes with the chance criteria (2.03) and (3.01). 
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With the distribution of mean motions, a different situation arises, since the 
proportion of commensurable pairs exceeds the chance value. It is found above 
that, of the 24 pairs of type [«, («+1)], 16 fall within the criterion « <¢y=0°011g. 
Hence the final probability P becomes in this case 


Px? x vy (=) P 
5 99 \24 
Orr. 
Thus, of the 17 pairs of the type [«, (« + 2)] in the actual solar system, 1-9 should be 
commensurable in this way simply as a result of the commensurability of the type 
[x,(a+1)]. In fact, there are three cases of commensurability of the type [«, (« + 2)] 
which can be attributed to this cause. The effect of including them in Table V is 
small, since by (2.03) these would by chance give 1-4 cases of commensurability. 

The above argument shows that the probability of finding a commensurable 
pair of type [«, («+ 3)] due to commensurabilities already present for [«, («+ 1); 
(x+1), (+2); and («+2), («+ 3)] is negligible. 

5. Significance of the results.—If p is the probability that, in a chance distri- 
bution, one pair of numbers will be commensurable within the limit «), and 
q=1-—p is the probability that that pair of numbers will be outside the same 
criterion, then the probability that, in a sample of v trials, there will be 7 pairs 
within the criterion is the (7 + I) th term in the expansion (q+ p)’, i.e. 


yO es (5.01) 
Table VI gives the probabilities, given by (5.01), of obtaining by chance the figures 
of Table V. 
TABLE VI 


Probabilities of obtaining the results of Table V by chance 


€o O'O11g 00089 00059 00030 O°0015 
Control distribution 0046 0°077 o'lig O'1l3 0°247 
Semi-major axes 0089 O°124 o'141 o°150 0°263 
Mean motions 0°00 42 00016 0'0014 0°0039 0°034 


In considering the significance of the figures of Table VI, itis instructive to quote 
Yule and Kendall (g) :—‘‘ Frequencies differing from the expected frequency by 
more than three times the standard error are almost certainly not due to fluctuations 
of sampling. They point to some departure of the sampling from simplicity, 
which may in turn point either to some flaw in the sampling technique or to causal 
effects in the universe itself.” In the present investigation the sample of 46 pairs 
is rather small for the application of the usual standard error formula (which is 
itself derived from (5-01) on the assumption of a very large sample). However, a 
frequency of three times the standard error is equivalent to a probability of 0-0044, 
and this value may be taken as a criterion of significance. 

Table VI shows that, whereas for the distribution of semi-major axes and for 
the control distribution, the probabilities are well above this critical value, for the 
distribution of mean motions, for all values of «, except the smallest, the prob- 
abilities are below the critical value. The agreements of the figures for the control 
distribution and for the distribution of semi-major axes with chance predictions 
prevent the results for the mean motions being attributed to flaws in the sampling 
technique. It may be concluded that the distribution of mean motions departs 
from simplicity, 
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This departure from simplicity may not be due entirely to a preference for 
commensurability. It is found that no ratio of mean motions exceeds 0-75, 
and since the mean density of fractions in the occupied range is somewhat higher 
than in the complete range of 4 to 1, this will affect the chance criterion. How- 
ever, with the limited sample of mean motions, it is not clear physically where 
precisely to limit the range. A chance criterion which takes this distribution 
into account is , 

ae (7N,)'= —Seongt = 42°8eqv. (5.02) 

4 7 
Using this criterion, the probabilities that the results for the actual mean motions 
are due to chance are still, in general, below the critical value. 

It should be noted that (a) the control distribution had the same limitation 
of range as the actual mean motions, and (4) the distributions are actually derived 
from a combination of rectangular ones, and the chance of an excess deviation is 
rather less than on the normal law*, thus re-inforcing the conclusion that the 
bodies in the solar system are distributed in such a way that pairs of orbits show a 
preference for commensurability of mean motion. 

6. Commensurability relationships with more than two bodies.—Some evidence 
supporting the belief that commensurability in mean motion is preferred is given 
by particular cases of the occurrence of commensurability relationships involving 
more than two bodies. The Laplace relation, noted before, is 

Ny — 3M, + 2ng=0. (6.01) 
This relation holds to at least 12 significant figures, viz. :— 
Ny — 2N, = 0-739 469 O9I, 
Ny — 2Ns = 0°739 469 0092. 
Indeed, according to Laplace’s analysis, this relation is exact. 

A similar relationship holds, although to a lesser degree of accuracy, for the 
four major Uranian satellites. If m,, m,, m3; and m, are the mean motions of 
Ariel, Umbriel, Titania and Oberon respectively, then, with the values of the 
mean motions adopted in this analysis, 
so that the relation ealliaionne-.itenlaiitaatai 

Ny — Ny — 2n3+n,=0 (6.02) 
holds to an accuracy of about 0-002 per cent of Ariel’s mean motion. As far as 
the present authors are aware, this relation has not been noted previously. 

By a rough extension of the criterion (2.03) to higher dimensions, it can be 
shown that the probabilities of (6.01) and (6.02) arising by chance are very small, 
of the orders of 10-* and 10-° respectively. An analysis of the mean motions in 
the solar system in groups of three and four bodies may be fruitful; it is hoped 
that such an analysis will form part of a future paper on this subject by the present 
authors. 

7. Discussion of result.—It is sometimes stated (10) that where two bodies have 
mean motions which are in a simple ratio to each other, the phenomenon of 
resonance appears. The perturbations are said to re-inforce each other in 
successive revolutions of the bodies, so that commensurability is rapidly departed 
from. The asteroid gaps are held to be supporting evidence of this effect, although 
it is noticeable that the few asteroids with periods exceeding seven years cluster 
about distances corresponding to mean motions of } and $ that of Jupiter (11). 


* The authors are indebted to the referee for drawing their attention to this point. 
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As is well known, in the development of the disturbing function when 7, 
and m, are almost commensurable, the long-period terms have large coefficients, 
since the factor ' : 

(41%, + tang) * 

enters into the expression for the amplitude (7, and 7, being small positive inte- 
gers) (12). The problem of resonance has been discussed by Brown and Shook 
(13). They find that their analysis contains nothing that precludes the existence 
of asteroids with osculating mean motions exactly twice that of Jupiter, although 
it shows that the elements of the orbits of such asteroids would suffer large 
oscillations. In passing, it may be remarked that the Trojan groups of asteroids 
provide an example of a stable orbit commensurable with that of Jupiter. 

The present analysis of the actual mean motions in the solar system indicates 
that, far from orbits with commensurable mean motions being unstable, there is a 
distinct preference for them. Even where the long-period inequalities are very 
large, as with the four Galilean satellites of Jupiter, the pairs of satellites Tethys 
and Mimas, Hyperion and Titan and Dione and Enceladus, and the planets 
Saturn and Jupiter, the bodies do not depart rapidly from commensurability. 

Two possible interpretations suggest themselves. The observed distribution 
of orbits may be explained by supposing either that the mechanism of formation 
of the planets and satellites was such as to favour orbits with commensurable 
mean motions, or that such orbits are relatively more stable over long periods of 
time than are neighbouring orbits, the planets and satellites thus tending towards 
commensurable configurations. 

In this connection, it should be noted that Laplace found a physical origin for 
the relationship (6.01) among the terms involving the squares of the disturbing 
forces. He showed that, provided that the relationship exists within a certain 
degree of accuracy, the mutual attractions of the satellites would subsequently 
render it rigorous. 
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FURTHER MEASUREMENTS TO DETECT A GENERAL 
MAGNETIC FIELD OF THE SUN 


H. von Kliiber 


(Received 1954 April 22) 


Summary 


The attempts made in 1949 to detect a general magnetic field of the Sun 
by an interferometric method were continued at the Cambridge Observatories 
during 1950, again supposing that the field would be of dipole form, with an 
axis approximately coinciding with the Sun’s rotational axis. The results of 
the measurements of about 700 interference knots are shown in Fig. 1. The 
m.s.e. of a single point is about +3 gauss. We conclude from these measure- 
ments that at the time of the observations, in 1950 July/August, such a genera! 
field must have been smaller than about 2 gauss and was not detectable by 
the present method. This is in agreement with the 1949 result. 





Attempts to detect a possible general magnetic field of the Sun by observations 
of the Zeeman effect of selected spectral lines (1) were continued at the Cambridge 
Observatories during the summer of 1950. The same excellent Lummer plate was 
used as in the previous year, again combined with the McClean solar spectrograph. 
The plane grating in this auxiliary spectrograph had an area of g x 12 cm and 600 
lines per millimetre. It was used in auto-collimation with a focal length of 450 cm, 
in the second order, and gave there a dispersion of about 2 A/mm and a resolution 
of about 100000. The quartz Lummer plate has a length of 200 mm, a width of 
30 mm and a thickness of 3-42mm. The range at (5250 is about 0-33 A. 

In 1949, the Lummer plate was used before the spectrograph, but in 1950 was 
placed behind the image plane of the spectrum. This arrangement had been 
tried experimentally in Potsdam in 1943 (2), on that occasion with a Fabry—Perot 
interferometer. Since about 1948, Treanor (3), Adam (4) and Priester (5) have 
used Fabry—Perot or Lummer plate interferometers very successfully in this way. 
The polarimetric optics were the same as described earlier, Wood’s method being 
followed, with a thick quartz plate (1, 6). Again the two Fraunhofer lines 
A 5250'22 Fe of Zeeman-type 0 (3)/1 and A 5263-32 Fe of type 0(5)/2 were used for 
the reduction. 

As far as possible the whole optical arrangement was protected from tempera- 
ture effects, due to the incoming light, by use of heat filters. ‘The Lummer plate 
was particularly carefully protected against temperature changes and care was 
taken as well that the whole quartz optics remained in perfect temperature 
equilibrium throughout the exposures, which lasted about Io to 20 minutes. This 
point has been checked repeatedly by examination of the fine structure of the green 
mercury line A5461. The general optical adjustment was verified both with this 
line and with terrestrial absorption lines. ‘The actual resolving power of the 
Lummer plate, tested in this way, was about 600000. Because of unfavourable 


weather conditions only a small number of really good plates could be obtained in 
1950 July and August. 





No. 2, 1954 Detecting a general magnetic field of the Sun 243 


As on former occasions, measures were made only at the two points of helio- 
graphic latitude + 45° on the central meridian of the Sun. In every case two 
photographic plates were combined in the manner suggested by Wood. For each 
of the two Fraunhofer lines five interference knots on each side of the plane of 
optical symmetry of the Lummer plate were measured and reduced in the 
customary way. ‘To obtain a clear picture of the results, these were immediately 
converted into magnetic field strengths and plotted in Fig. 1. Each of the points 
of this figure depends on the measurements of 40 interference knots and represents 
the mean from 160 single measurements. 
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Despite the fact that the interference knots seemed to be slightly better defined 
than in the previous year’s work, thestandard error (m.s.e.) of a single measurement 
was not improved, and remained between +0-002 and +0:005mm. The 
interference knots were imaged on the photographic emulsion by an excellent 
Tessar lens of 25cm focal length, placed behind the Lummer plate. The 
separation of two corresponding interference knots varied between 18 and 23 mm. 
An increase of this separation by 0-00I mm corresponded on the average to an 
increase of the corresponding wave-length by about 0-0003A. Furthermore, 
taking into account the Landé factor of the two Fraunhofer lines which were used, 
and of the wayin which Wood’s method works, we found thata residual of 0-001 mm 
between a corresponding pair of plates was equivalent to about 2-5 gauss. ‘The 
m.s.e. for each pair of plates, as far as internal consistency alone is concerned, is 
about +3gauss. ‘The results in Fig. 1 summarize altogether more than 3000 
single measurements. 

Looking at Fig. 1 there can be no doubt that the single points scatter around the 
zero line only within the permissible errors, i.e. single values show no significant 
deviation from zero. If we form the mean of all measurements for each hemi- 
sphere, we find for the investigated point on the central meridian of the northern 
hemisphere the value — 0-4 gauss, and for the corresponding point on the southern 
hemisphere +0-15 gauss. The positive sign denotes a north-seeking magnetic 
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TABLE I 


Sunspot groups from Fig. 2 for which magnetic field strengths were measured at Mt. Wilsor 
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pole. Both values, however, can be considered only as a purely formal result of 
calculation, without any real significance. The only conclusion we can draw is 
that during the observing period 1950 July to August no permanent magnetic 
field of the Sun larger than about 2 gauss can have existed. This assumes that it 
would be possible to establish the existence of such a general field by the study of 
longitudinal Zeeman effects and that the Sun possesses a dipole field orientated 
approximately with the rotational axis. Alfvén (7) has pointed out possible 
limitations on results obtained in this way. 

As already mentioned in (1), it is possibly important to consider the general 
activity of the Sun at the time of observations. In the year 1950 the Sun was still 
rather active: the Ziirich relative numbers in 1950 July/August reached values up 
to 130. In Fig. 2 we reproduce a synopticchart of the sunspot activity at the time 
of the observations, using data from (8). The crosses indicate the dates and the 
heliographic positions of our measurements. It can be seen at once that the 
northern hemisphere of the Sun was considerably more active at this time than the 
southern hemisphere. Furthermore, the present measurements show more 
scatter on the northern half of the Sun than on the southern half, but the number of 
observations is so small that this may not be significant. Photographic measure- 
ments of the magnetic field in sunspots for the period of the Cambridge observa- 
tions are this time, unfortunately, not available at the Potsdam Observatory (g), but 
in order to give some assessment of the importance of the main centres of 
disturbances shown in Fig. 2, we have combined in Table I the maximum 
magnetic field strength measurements at the Mount Wilson Observatory (10), 
together with the heliocentric longitude and latitude of the spot area, the Ziirich 
classification (8) at the time of the greatest activity of the group, and the Mount 
Wilson number of the source. 


Again we are deeply obliged to Professor S. Tolansky, F.R.S., of Royal 
Holloway College (London), for the kind loan of the Lummer plate used for these 
observations. 


The Observatories, 
Cambridge : 
1954 April 22. 
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THE SPECTROSCOPIC DOUBLE-LINED BINARY HD 77464 
M. W. Feast 


(Communicated by the Radcliffe Observer) 
(Received 1954 May 3) 


Summary 


HD 77464 is shown to be a spectroscopic double-lined binary consisting 
of two almost equal B2 V stars. The binary has a period of 64-892 and 
nearly circular orbits. ‘The minimum masses are 5°61 and 5:41 times the 
solar mass and the minimum semi-major axes 11°53 < 10° and 11°94 X 108 km 
respectively. 





Introduction and observations.—The first part of the programme on the deter- 
mination of the radial velocities of southern O and B type stars at the Radcliffe 
Observatory has recently been completed (1). As was to be expected, a consider- 
able number of variable velocity stars were discovered, a few of which show double 
lines. In view of the programmes already in hand it is not intended, for the time 
being at least, to study the single-lined binaries in detail. However, because of 
the extra information that can be obtained in the few cases where the spectra of 
both components are visible, it has been decided to obtain sufficient observations 
for the determination of the orbital elements of these systems. 

The first double-lined binary discovered at the Radcliffe Observatory was the 
star HD 77464, R.A. 8"59™-1, Dec. — 51°22’ (1950), 6™-6pg(HD). The HD 
type is B3 whilst the revised Radcliffe type (1) is B2 V for both components. 
Between 1952 March and 1954 February, 34 spectrograms of the star were obtained, 
27 of which showed doubling of the lines. The observations are listed in Table I. 

The columns of this table are as follows: (1) A serial number. (2) The date 
(U.T.). (3) The Julian date. (4) The phase calculated from the finally adopted 
Ty, the time at which V, equals y and dV,/dT is positive. The phase is given as a 
fraction of the adopted period and not in days as is sometimes done. (5) The 
velocity, V,, of the primary. (6) The number of lines used in forming V,. (7) 
The velocity, V, of the secondary. (8) The number of lines used in forming V4. 
In the cases where the lines were measured single the velocity is entered in brackets 
in column (5). The single-lined phases are helpful in fixing the period. They 
have not been used in computing the other elements. (9) The weight assigned to 
each plate, both components having equal weight. (10) The Ca* K line velocity. 
Both stars show a stellar K line and an interstellar K line is also present. When 
the three lines were resolved the stellar K velocities were used in forming the 
stellar velocities of columns (5) and (7), and the interstellar velocity is given in 
column(10). When the K line was measured single its velocity is entered in column 
10 in brackets. In some cases the line was found to be double, consisting of a 
resolved stellar component (used in forming the stellar velocity) and a blend of the 
other stellar component and the interstellar component. The velocity of this 
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blend is also entered in column 10 in brackets. The bracketed velocities were not, 
of course, used in any of the calculations. Frequently one or two stellar lines were 
unresolved at double-lined phases; the velocity derived from these lines is omitted 
from the table. (11) The camera used. Details of the various cameras of the 
Cassegrain spectrograph of the 74-inch Radcliffe reflector with which all the 
observations were made are givenin(1). Inthe present case either the c(49 A/mm 
at Hy) or b(29 A/mm at Hy) cameras were employed. Also in column (11) are 


TABLE [ 


(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 
No. Date (U.T.) J.D. Phase V; N, V, N. Wt. Ca*+K Camera 
2 434 000+ km/sec km/sec km/sec 
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5227190 0°047 + 72°4 
522°326 0066 + 793 
756°476 o-o41 (+ 27°0) 
769°456 0924 — 26:0 (+11°7) 
778°451 o'229 +1482 10 : +25°8 
781-381 0944 (+ 27°38) 9 (+ 6-1) 
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letters denoting cases where observers other than the writer obtained the spectro- 
grams. They are :—T (Dr A. D. Thackeray), W(Dr A. J. Wesselink) and M 
(Mr A. Morrisby, Cape Staff). All the measurements were carried out by the 
writer. A discussion of the methods of measurements of the spectrograms, the 
system of wave-lengths employed and the types of photographic plates used will 
also be found in reference (1). 
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The two components of the binary are of almost identical spectral type and 
brightness. However, the observations are such that the approximate period 
can be fixed without ambiguity. The designation of one star as the primary and 
the other as the secondary is merely a matter of convenience and is not to be taken 
to imply anything about the relative brightness of the two stars. 

The derivation of the elements.—Most of the observations were obtained in the 
early months of 1953 and clearly the values of all the elements except the period 
depend mainly on these observations. The 1953 observations yield an approxi- 
mate period but the accuracy of this is greatly improved by the few observations in 
1952 and 1954. ‘Thus in a solution for the orbits by the method of least squares 
the error in the period is, in the present case, practically independent of the errors 
in the other elements, so that it may safely be determined alone first and left fixed 
in the rest of the work. This greatly facilitates the work, especially in that it 
allows the formation of normal places. The period was therefore determined 
graphically by trial and error and was fixed at 64-892. The mean error cannot, of 
course, be determined by this method but periods of 64-890 and 64-894 are both 
definitely inferior to the adopted one. When the observations were assembled on 
this period the preliminary elements 7), e and K,, given in Table III, were 
deduced graphically, the symbols having their usual significance. Since the 
velocity curves clearly showed that the orbits were practically circular, Sterne’s 
second method (2), in which differential corrections to preliminary circular orbits 
are obtained, wasemployed. However, use was made of a modification suggested 
by Dr A. J. Wesselink, which considerably reduced the amount of computation 
without impairing the accuracy. It is easily shown that for any binary the relation 


V,=y(1+ K,/K,)— K,/K,. V2 


must hold, so that from it y and K,/K, may be obtained. To solve this equation 
the method of least squares could be employed (though the variables would need 
to be changed, probably to (V,+ V,) and (V,— V,), in order that the observational 
error should have a noticeable effect on one variable only). However, the method 
of solution by grouping discussed by Jeffreys (3) is very suitable in the present case. 
If the observed values of V, and V, were uniformly distributed between their 
extreme values the efficiency of the solution would be = that of the complete 
least-squares solution, as was shown by Eddington (see 3). Since in fact values of 
V, and V, near their extremes predominate, the efficiency in this case is even higher 
and the method is obviously quite satisfactory. The values of K,/K, and y thus 
found are given amongst the final elements in Table III, together with their mean 
errors (computed for a least-squares solution). ‘These values are then used to 
reduce the V, values to the V, velocity curve and the solution reduces to the treat- 
ment used in the case of a single-lined binary except that y is not varied. Sterne’s 
second form of the equations of condition (2, equation (6)) becomes in this case: 


8V =cos L.8K,+cos2L.Ky.e.cosw,+sin2L.Ky.e.sinw,+sin L.u.K,.5T, 


As is usual, normal places were formed to facilitate the solution of these equations. 
The 27 observations at double-lined phases were reduced to the 16 normal places 
givenin Table II. 

The columns of Table II are as follows: (1) The serial numbers from Table I 
showing from which observations each of the normal places was obtained. 
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(2) The weights. (3) and (4) The weighted mean velocities. (5) The mean, V, of 
V, and V, reduced to the V, curve as discussed above. (6) The weighted mean 
phase of each normal place from the preliminary 7). (7) The residuals of V from 
the final orbit (in the sense observed minus calculated). The 16 equations of 
condition were solved by least squares, yielding four normal equations for the 


TaBLe II 

(1) (2) (3) (4) (5) (6) (7) 

. Final 
“ae Vv V V residuals 
from Wt. - . ; Phase 

= km/sec km/sec km/sec (O—C) 
lable I 

km/sec 
+ 72°4 — yg 10°7 0°038 —t-9 
- 79°3 24°4 24°9 0-060 o"7 
+ 106°9 60'9 57°4 O'1I2 +10 
+128°8 68-6 72°6 O'144 tr 
+ 144°5 86°6 89°8 o°185 
+1504 97°7 98°4 0°229 
+143°6 94°9 93°5 0-287 
+1348 ah wal 79°8 0°325 
+116°8 — 87:8 76:0 0°337 

49°7 * 932 + 100°8 0585 

75°6 +1342 +134°7 0°639 

71°0 + 143°6 +137°0 0°661 

97°8 ~157°3 +1$7°7 0°737 

92°9 +149°8 +I51'4 0-783 

55°3 +112°8 +113°5 0°874 

26°0 + 81°5 + 82:6 O'9I5 


29 
39, 13 
18 
20 
9 
5 
22 
23 
10 


I 
4 
3 
4 
4 
4 
3 
3 
3 
I 
2 
3 
4 
4 
2 
2 


unknowns and other equations for their meanerrors. The final elements and their 
mean errors are entered in Table III. The eccentricity is hardly significantly 
different from zero. The other elements do not appear to call for any particular 
comment. It was found that the mean error of an observation of average weight 
was +3°3km/sec for each component. The observed velocities are shown 
plotted in the figure, in which the curves were computed from the final elements. 


TABLE III 


Element Preliminary value Final value and mean error 
64-892 
JD 2 434 094628 JD 2 434 0947567 + 0-008 
143° + 22° 
323° + 22° 
° 0:016+0°006 
K, 122'24+ 0°63 km/sec 
K, 127°1 km/sec 12663 + 0°63 km/sec 
K,/K,=m,/mz 1°036+0-011 
y -27-9+1°'1 km/sec 
Interstellar Ca II -13°9 km/sec 
a, sini II 530 000 km=0°'077 1 A.U. 
a, sint II 940 000 km=0°079 9 A.U. 
m, sin® 7 5°610 
my sin® i 5410 
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Fic. 1.—The velocity curves of HD 77464. 


@ = Velocities at double-lined phases. 
x =Velocities at single-lined phases. 
O=Cat* K interstellar velocity when the line is completely resolved from stellar K. 


The curves are drawn according to the final elements. 
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MOTIONS IN THE SUN AT THE PHOTOSPHERIC LEVEL 
V. VELOCITIES OF GRANULES AND OF OTHER LOCALIZED REGIONS 
H. H. Plaskett 
(Received 1954 May 5) 


Summary 


Velocities and surface brightness in three Oxford spectra, showing granu- 
lation, have been measured. The results are in good agreement with those 
earlier found by Richardson and Schwarzschild in a Mount Wilson spectrum. 
In two of the three Oxford spectra a significant correlation coefficient of 
r= —o-'28 has been found between the surface brightness and the sight-line 
component of the velocity, showing that the bright granular matter is moving 
vertically upwards. Assuming that the top of the granulation zone occurs at 
an optical depth of unity, 0-75 km sec~! is a lower limit to the velocity of 
ascent of the bright granular matter. The smallness, or in one spectrum 
the absence, of any correlation between surface brightness and velocity 
suggests that the granular velocities are being partially or completely masked 
by some other perturbation of velocity. Removing from the observed 
velocities that part due to granulation, the perturbing velocities are found to 
originate in an oscillatory velocity field with an average period of 5” of arc 
and an amplitude of the order of 0-5 km sec~!. 





Introduction.—From a beautiful spectrum, showing well-marked spectra of 
granules, taken with the 150-foot Tower Telescope at Mount Wilson in 1949, 
, Richardson and Schwarzschild (15) have measured a random, turbulent velocity 
of 0-37kmsec"!. Though there can be no question of the accuracy of their 
measurements, the velocities which they find are so small and apparently so random 
that confirmatory measures and a new analysis may not be out of place. Between 
1947 and 1950, 224 spectra of granulation were obtained with the Oxford solar 
equipment (1x, 12). The present paper contains in succeeding sections an 
account of the observation and the measurement of three of these spectra, 
secondly a determination of the velocity of solar granules, and finally evidence 
for the existence of other localized velocity-fields in the photosphere. 


1. Observation and measurement 


1.1. Observation.—Because of the difficulty of estimating ‘“‘seeing” at the 
instant of exposure, early work on the spectra of granules (10) proved very 
wasteful. In the present investigation this difficulty was avoided by the use of 
a guiding plate (13) with which the solar image on and near the slit could be 
viewed before and during the exposure with an angular magnification of 390 
times. The seeing was estimated from the appearance of a spot seen in the field 
of view of the guiding plate; the spot was held so that no part of it was within 
50” of arc from the slit. 

Of the 224 spectra thus obtained under conditions of moderate to good seeing,. 
three were selected as showing the best contrast and the finest detail in the streaks. 

17* 
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due to granulation. Observational details are given in the self-explanatory 
Table I. The heliographic coordinates of the centre of the slit, and the direction 
cosines of the sight-line to the corresponding point on the Sun are found by the 
methods described elsewhere (13). The fourth spectrum listed in the table is a 
standard spectrum of the centre of the disk made with the solar image focused 
34cm in front of the slit. It is against this standard centre spectrum that the 
three preceding granulation spectra are measured. 


TABLE I 


Observational details 









































Exp. Slit Centre slit 
Date U.T. spectrum time 
Width Height L—Ly, B 
sec mm mm 
1949 Sept. 11-4306 3 0030 12°7 —~2°52 4°74 
1949 Sept. 11°4319 3 0030 12°7 —2°§2 4°74 
1950 May 11°7277 I 0-040 12°7 —6:25 +-1°65 
1950 April 9°3392 I 0°040 12°7 0°00 6-02 
Date U.T. spectrum cos y, COS Y2 COS Y3 
1949 Sept. 11°4306 —0°0439 —0'0440 —o-'9981 
1949 Sept. 11°4319 —0°0439 —0°0440 —o-'9981 
1950 May 11°7277 —0'1092 +-0°0816 —0°9907 
1950 April 9°3392 ©0000 00000 — 1°0000 




















Spectra centred at 5900 A on plate. 

Emulsion : Ilford Rapid Process Panchromatic. 

Dispersion : = 0-68 mm A~’. 

I mm on spectrum normal to dispersion= 10"-41 of arc on Sun. 
I mm on slit =10"+34 of arc on Sun. 


1.2. Measurement of velocities. —Prior to measurement each of the four spectra 
was marked by two sets of 12 fine scratches exactly I mm apart in a direction 
perpendicular to that of the dispersion. One set of these scratches was at 
5884 A and the other at 5959 A. The scratches were made with a special dividing 
machine fitted to a Cambridge micrometer. 

The spectrum to be measured is mounted on a Hilger Photo-measuring 
Micrometer, and carried on a subsidiary stage by which the spectrum can be moved 
perpendicular to the direction of dispersion by a micrometer screw. A mask 
in the plane of the cross-wire of the micrometer limits the visible height of the 
spectrum to 0-145mm. ‘The scratches on the spectrum facilitate its orientation 
so that the direction of motion of the microscope is exactly parallel to the direction 
of dispersion, and also serve in calibrating the vertical micrometer. 

Measurement is made of the seven atmospheric and six solar lines listed in 
Table I]. The wave-lengths with asterisks are those recommended by the 
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Commission on Wave Length Standards (2); those without asterisks are 
taken from St. John’s Rowland Revision, corrected as recommended by the 
Commission (3). 

The standard centre spectrum was measured, direct and reverse, only at the 
heights ofthe 12 scratches. The positions in millimetres of the 13 lines in Table II, 
referred to a common origin and averaged over the 12 heights, gave the standard 
settings of these lines for a point distant 0-205 mm west of the point corresponding 
to the centre of the slit. Assuming the photospheric matter over the height of 
the slit has on the average no vertical component of velocity, the sight-line 
component of the velocity of this mean point in the standard spectrum is given by 


Cy is the correction for the component of the observer’s velocity in the sight-line, 
U, is the solar rotation and y,(0) is the appropriate direction cosine—all with 


reference to the mean point, 0-205 mm to the west of the centre of the slit. C, and 
yi(0) are computed by the methods described elsewhere (137). 


TABLE II 





Atmospheric water 


: Solar lines 
vapour lines 





Wave-length Rowland Wave-length Rowland aa Velocity 
1.A.U. no. L.A.U. no. _— factor 





km sec™!mm~! 
5885 977 
5891-661 
5901 468 
5919°054* 
5919°644* 
5924°272 
5932 °092* 


5892-883 * : Ni 73°67 
5905 ‘680* 4 Fe 74°29 
5916°257* Fe 74°80 
5927°797* Fe 75°36 
5929682 Fe 75°46 
5930°191* Fe 75°48 


nau Oat w 


ue 
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Each of the three granulation spectra was then measured, direct and reverse, 
at each tenth of a millimetre over the clear height of the mask (0-145 mm) from 
0-5 mm to I2 mm inclusive, that is at 116 heights. At any height the difference 
between the measured position of an atmospheric line and its standard position in 
the centre spectrum, both referred to a common origin, gives the correction appro- 
priate to this wave-length to reduce the dispersion of the granulation spectrum to 
that of the standard centre spectrum. A correction curve (in practice a straight 
line) may then be constructed for the atmospheric lines, and from it the appro- 
priate correction for any solar line at this height may be read off. The difference 
between the so corrected position of a solar line at a given height and the standard 
position of the same line in the centre spectrum then gives the displacement of 
the solar line (in millimetres of the standard centre dispersion) due to the 
sight-line component of the velocity relative to that of the standard centre spec- 
trum. ‘This displacement in millimetres can be converted into km sec~' by 
multiplying by the appropriate velocity factor given in Table II. The velocity 


+ On p. 423 of this paper it should be noted that sin 59’ is given, not by Z, but by 
Z/[X2+ Y?+Z?}!/2, 
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factor is found from the dispersion of the standard centre spectrum, given by the 
Hartmann constants computed from the standard settings of the atmospheric 
lines of this spectrum and the Doppler factor appropriate to the wave-length. 
Denote the velocity component so found from the mean of the six solar lines at 
the given height by V*. Then the sight-line component, V, of the solar velocity 
of the photospheric matter at the point corresponding to this height in the spec- 
trum is given by 

V=V*+C,—U, cosy, + Vo, (2) 
where C, is the correction for the sight-line component of the observer’s velocity, 
U, the solar rotation and cos y, the appropriate direction cosine, all for the given 
height in the granulation spectrum. FH 

In the second column of Table III are given values of V, the average over the 

116 heights of the so found values of V. Since it was anticipated that V would 
be zero, it was first suspected that an error had been made in determining the 
corrections, C, and C,, for the observer’s velocity. However, a re-calculation of 
these corrections by independent methods revealed no errors. The consistently 
positive sign of V for all three granulation spectra then suggested the hypothesis 
that the photospheric matter had an upward motion when the standard centre 
spectrum was exposed. In this event equation (I) should read 

Vo=V_ —Cy+ U,cosy,(0), (3) 
where V,’ is the sight-line component of the photospheric velocity. From the 
standard positions of the lines in the centre spectrum and the wave-lengths in 




















TasBce III 
Mean velocities of spectra 

Date U.T. spectrum V V; V; 
km sec™?! km sec™! km sec™! 
1949 Sept. 11°4306 +0o°18(7) —o'15(7) +0:02(8) | 
1949 Sept. 11°4319 +0°15(0) 0°19(4) 0°00(9) | 
1950 May 11°7277 +0°29(9) ~0°04(5) +o'14(0) | 
1950 April 9°3392 —0°34(4) —0°15(9) | 
| 
“Vieans —o'18(5) 0'00 | 
\ 








Table Il we may find the velocity, V4, in the usual way from the Hartmann 
dispersion formula. It turns out to be —0-08(4)kmsec™, or from equation (3) 
V,' = —0°34(4)kmsec™*. If the apparatus function of the Oxford spectroscope 
were as asymmetrical at 5900 A as it is at 4300A, the velocity V, found in this 
way would probably be systematically in error by an amount depending on the 
difference of the mean equivalent widths of the solar and the atmospheric lines 
{12). However, a recent investigation by R. J. Bray shows that in the region from 
5461 A to 5945 A the violet asymmetry of the apparatus function is much reduced. 
The values of V, and V,’ thus found are probably substantially correct. 

In the third column of Table III are given the averages, V;, for each granulation 
spectrum of the values of V when the new value, V,= —0-08(4)kmsec~, is used 
in equation (2). The fourth entry in this column is the corresponding value, 
V,', for the standard centre spectrum. If we assume, as seems reasonable, that 
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the velocity at the centre of the disk averaged over a number of days is zero, then 
the average for the entries in the third column, namely —0-18(5)kmsec"!= 
— 0°18(5)/50-72 A = — 0-003(6) A, must be due to the solar wave-lengths in Table II 
being too large relative to the wave-lengths of the atmospheric lines by 0-003(6) A. 
This is in good agreement with Miss Adam’s (1) interferometric determination 
at 6100 A where, with atmospheric lines as standards, she finds the I.A.U. wave- 
lengths of solar lines too large by 0-0033 A. Applying a correction of —0-003(6)A 
to the wave-lengths of the solar lines given in Table II, we finally obtain the 
averages, V,, given in the last column of Table III. From these values it is 
evident that we may on occasion have substantial velocities of ascent or descent 
over regions subtending more than two minutes of arc at the centre of the disk. 

In the present investigation we are interested in the velocities of small regions 
along the slit in excess or defect of the average over the whole length of the slit. 
These required velocities are given by V—V, where V is the velocity directly 
measured at a given height and V is the average of all such velocities over the 
whole height of the spectrum, that is the quantity given in the second column of 


Frequency 
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Fic. 1.—Frequency histogram of the distribution of 348 sample r.m.s. errors. 
The smooth curve is the theoretical sampling distribution. 


Table III for each granulation spectrum. These required values, V—V, which 
are given in the second column of Table IV for each of the 116 heights in each of 
the three granulation spectra are of course independent of the corrections for the 
observer’s motion, of the adopted wave-lengths of the lines, or of any motion 
shown by the standard centre spectrum. 

1.21. Errors of measurement.—The r.m.s. error of a single solar line may be 
found in the usual way from the deviations of the line-velocities at a single height 
from the mean of the line-velocities at that height. The distribution of the 348 
sample r.m.s. errors from 116 heights in each of the three granulation spectra is 
shown in Fig. 1. The smooth curve is the theoretical sampling distribution on 
the assumption that the hypothetical, infinite population of line-velocities has a 
Gaussian distribution. The observed and theoretical distributions are in good 
agreement. From the observed sampling distribution we find the average 
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r.m.s. error to be s=0°354kmsec™!, or the r.m.s. error o of the hypothetical, 
infinite population of line-velocities, given by 


whe 2 I(n/2) 
= A n—-1 Tin—12\” 


will be o=0-373kmsec™, since in our casen=6. The r.m.s. error of a itself 
is approximately o/\/2m=o-107kmsec"!. The anticipated r.m.s. error of 
V —V, the mean of the six line-velocities at a single height, will be 
0°372 
v6 





=0°152km sec}. 


The dispersion in the values of V — V, given by the standard deviation of the 
values in Table IV, is 0-20(3), 0-26(4) and 0-29(5)kmsec™! for the spectra of 
1949 September 11-4306, 114319 and 1950 May 11-7277 respectively, or 
0°25(6)km sec”! for all three spectra together. Assuming for the moment that 
the velocities, V—V, along the height of the spectra are constant, this standard 
deviation of 0-25(6) implies an r.m.s. error of the hypothetical, infinite population 
of line-velocities of 0:25(6)\/6=062(7). This differs from the r.m.s. error of a 
line-velocity, deduced in the previous paragraph, by 0-62(7) — 0-37(2) =0-25(5), 
or by approximately two and one half times the r.m.s. error of itself. This is 
well above the 5 per cent level of significance; we may therefore reject the 
hypothesis that the values of V — V are constant, and conclude that the dispersion 
in V — DV is primarily due to differences in photospheric velocity along the height 
of the slit. 

The measures of V — V by Richardson and Schwarzschild (15) have a standard 
deviation of 0-26(2), agreeing very closely with the value of 0-25(6) found in the 
present investigation. They give the probable error of measurement of V —V 
as 0°035 kmsec"!, that is, an r.m.s. error of 0-052. This is roughly one-third of 
the corresponding quantity, 0-15(2), estimated for the present measures—a 
ratio which, since it corresponds fairly closely with that of the two spectroscopic 
dispersions employed, might have been anticipated. 

1.3. Measurement of surface brightness.—Of the three granulation spectra 
listed in Table I, the first two are on a plate carrying a wedge (W. iv) and a blank 
spectrum, each made with the same exposure time as that of the granulation 
spectra. The calibration curve at 5885 A for these two spectra, corrected for 
non-uniformity of illumination over the height of the slit, is determined in the 
usual way (10). The third granulation spectrum is on a plate carrying no 
calibration. There are available, however, two plates from the same box, both 
exposed on the same day and for the same time as the granulation spectrum and 
developed under the same conditions with Kodak Press Contrast developer. 
Each of these plates carries a wedge and a blank spectrum, and the calibration 
curve at 5885A from these two plates may probably be safely applied to the 
granulation spectrum. 

Between 5884-4A and 5885-7A, with a length of approximately 0-9mm on 
the Oxford spectra, is a stretch of effectively pure continuous spectrum containing 
only two lines of Rowland Number —1 and four lines of —2. Three separate 
and non-overlapping microphotometer (9g) tracings perpendicular to the disper- 
sion were made of this region in each granulation spectrum with a slit width of 
0-057 mm (perpendicular to the dispersion) and a slit height of 0-28 mm (parallel 
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to the dispersion). The deflections of each tracing were measured at each of the 
116 heights at which velocities were measured, and the deflections converted 
into incident light fluxes by the calibration curve. The mean of the three 
determinations of incident flux for each height should give a value in which 
accidental error, due to fluctuations in the number of silver grains, is reduced. 
The incident flux measured at a height mm on the spectrum is proportional 
to the illumination, J(€), of the slit by the image of the Sun times A(€), the trans- 
mission of the slit. For an ideal slit B(€) would be a constant independent of 
height, but in any practical case, particularly where a narrow slit width is used 
and slight irregularities in width are magnified by diffractional spreading, A(£) 
must be determined. The following method, a modification of that earlier 
employed (10), has been used. On each of the two days on which the three 
granulation spectra listed in Table I were exposed, there were available four 
other calibrated spectra taken with the same slit width and height at the same 
part of the solar image, but showing less marked granulation. Single tracings 
perpendicular to the dispersion were made of each of these four spectra at 5885 A, 
and the incident flux at each of the 116 heights was found from the measured 
deflection and the appropriate characteristic curve. Let J;(&) B(€) be the incident 
flux measured at height € on the ith of these four spectra, expressed in units such 


that 12 
= JAB) 


16.5, 

is unity. In this expression we may put 
J{E)=JSolEMt + €i), 
where J,(€) is the mean illumination and ¢; is a small number, absolute value not 
greater than 0-05 and as frequently positive as negative, which arises from the 
presence of granulation, itself variable with time. Thus the average of the 
incident fluxes at height € on the four spectra will be 
4 


£ = SAE)BE) =JolE)BEMT + 2%€1)] = FolE)ALE).- 


Hence if we divide the measured J(£)A(€) of one of the granulation spectra of 
Table I by the just determined J,(£)8(€) for that day, we shall obtain kJ(€), where 
k is a constant and where variable slit transmission with height and any constant 
trend in the surface brightness of the Sun have been removed. The values of 
kJ() thus found are listed in the third column of Table IV. 

The standard deviation of the surface brightnesses iisted in Table IV is 
0-023, 0-026 and 0-021 for the granulation spectra of 1949 September 11-4306, 
11-4319 and 1950 May 11-7277 respectively. From the frequency of maxima 
and minima, not counting some of the less marked turning points as being possibly 
due to silver-grain fluctuations, we find the mean diameter of the bright maxima 
(‘‘ granules”) to be 4”-2, 3"-7 and 3”-6 of arc respectively for the three spectra. 
These results may be compared with those found by Richardson and Schwarz- 
schild (15). The standard deviation of their measured surface brightness is 0-22 
magnitudes, or 0-021 in surface brightness units—a dispersion comparable with 
that found for the Oxford spectra. On the other hand, the mean diameter of 
their bright maxima is 1"-8 of arc, roughly one-half the size of those measured 
on the Oxford spectra. This superior solar resolution is due in part to a solar 
image twice the size of that given by the Oxford telescope, but probably even more 
to better seeing when their spectrum was exposed. 
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TABLE IV 
Values of V-V and kJ 














é V-—V kJ é V—V kJ é V—V RJ 
mm km sec™! mm _ | km sec”! mm _ | km sec™! 
1949| Sept. 11°14306 50 -O'17 1°007 9°7 +017 0°967 
I - 7 1-009 ‘8 + ‘17 0°933 
o's FOrrs 0994 2 04 1028 ‘9 + -41 1 "007 
‘6 r *§8 | 0-999 3 etl 1054 
"7 "42 0°993 "4 — ‘27 1052 10°0 02 1°O17 
8 + "4 1°009 5 — EE 1°043 “s + 1-009 
‘9 "24 I‘O17 6 — -06 1°O41 2 + -o1 0'996 
se + "49 1°027 _ + ‘20 0°986 
1°o 13 1°026 8 Il 1‘OI7 4 + Os 0°973 
‘I “24 1°035 ‘9 ay 1013 5 ‘47 0°957 
2 ‘22 1°044 6 ‘25 0°948 
3 7 1°044 6-0 13 1‘O12 7 "30 0°954 
"4 05 1°033 ‘I 22 1°029 8 "14 0°963_ 
“5 — +30 1'038 3 “_— 1°042 ‘9 — 04 0-967 
6 — °27 1:028 ae "2 1034 
“" — °33 I'O14 “4 "O02 1°029 I1°o + *82 0978 
8 — 18 1:006 “5 — 18 1°O10 “I + -o08 °°990 
‘9 + +03 1'O13 6 — 03 0°984 2 00 0996 
se — ‘02 0-982 4 + -28 1-002 
20 — ‘07 1'016 8 ‘II 0'980 "4 + ‘25 0987 
18 I°OI! ‘9 - +19 °'990 < + +17 0-980 
s “45 1°007 6 + °13 °°992 
3 ‘28 | 0984 7° ‘39 | 0999 7 + 06 | o-g91 
"4 + -12 0'987 ‘I — *30 0-981 8 + ‘OI 0987 
5 13 0'987 2 + 07 o'981 ‘9 — ‘14 0°99! 
6 ‘Ol 0°997 3 ‘14 | 0°997 12°0 ae 0-992 
oo — +34 I'O14 "4 — ‘28 I‘O1l 
8 03 1009 5 — ‘or 1009 1949] Sept. 11°14319 
9 “9 1'O13 6 + -08 1°000 
a "29 0-992 o's +0'27 1004 
3°0 — 03 I‘O19 ‘8 25 0985 6 + ‘15 1°O12 
“< + -08 1-006 ‘9 — ‘I§ 0'977 oe + ‘12 1-009 
2 + *14 0993 8 + 1 ‘004 
‘3 + -19 0°993 8-0 08 0-976 ‘9 + -08 1-014 
‘4 08 0°995 I + ‘09 0'977 
5 18 1:008 2 — +16 0956 1°0 — ‘04 1°029 
6 ‘OI 1:006 3 07 o'951 | + +27 1°022 
“y '§3 1°00! "4 "3 0°957 ‘2 + -38 I°O15 
‘8 08 0°997 5 ‘2 0968 3 + -46 1‘O17 
‘9 “35 1006 ‘6 — -28 0983 ‘4 + ‘29 1°030 
"7 + °17 0985 *s + ‘oI I°O17 
4°0 ‘20 1°007 ‘8 ‘19 0'969 6 + °39 1'O12 
I 03 1-018 ‘9 07 0960 5 + 22 0-988 
“2 08 1035 ‘8 + 33 1005 
3 24 I'O19 90 + -18 0'978 ‘9 + ‘21 1016 
‘4 ‘25 1004 a 09 1002 
5 ‘O7 °0°995 2 "19 o'992 2°0 + ‘03 1-013 
6 ‘10 ©-999 3 “14 0°978 ¥ + ‘02 1006 
"7 08 1°005 "4 — -18 0'977 2 — ‘12 1°005 
8 "05 1-006 *s + °22 o'981 “7? + °24 0-989 
‘9 + 0°05 1°009 ‘6 +0-02 0982 ‘4 —0'03 0-989 
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TABLE IV (cont.) 
Values of V— V and kJ 
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TABLE IV (cont.) 


Values of V— V and kJ 











é V—V kJ é V-V kJ é V-V kJ 
mm km sec™? mm _ | km sec™! mm_ | km sec! 
4°6 —0'13 1°029 71 —0°36 1°005 9°7 0'00 o'971 
sy — 32 1°036 ‘2 + ‘41 I°O10 8 “$9 0981 
8 — "94 1°023 3 ‘OI 0°997 ‘9 — 2 o'961 
‘9 — 32 1003 ‘4 + +12 0989 
5 + *54 I‘OoI 10°0 — °°953 
5°0 + 37 1°005 6 + +34 1025 ‘I — "44 0958 
oS — ‘21 0995 | + +18 I‘O10 2 + ‘or o°951 
oF + 39 1002 8 "04 1‘OI! 3 — 18 0-958 
~ + 26 *1-O12 ‘9 66 I°O15 4 — ‘09 0965 
“4 , 1005 5 — 06 | 0-987 
5 + ‘22 0°998 8-0 ‘12 I‘OII 6 — ‘21 I°O15 
“6 + +18 1°008 “! + +42 1038 "7 — ‘02 0:997 
Se + 95 1°026 2 — +38 1°044 8 ‘OI 0980 
8 + 43 1'016 3 — +25 1°024 9 "12 0°973 
‘9 45 1‘O17 ‘4 07 I°O14 
5 34 I‘Oll Ilo 04 0965 
6-0 + ‘60 0996 6 — "O4 1016 “I — +29 0-964 
1 + -29 | o-981 7 — 44 1009 2 — +30 | 0976 
2 ‘29 | 0°979 8 — +31 | 07984 3 — 66 | 0°993 
3 ‘17, | 0°977 ‘9 43 | 07988 “4 — +28 | 0983 
"4 “31 0968 5 — ‘22 0'971 
5 + °20 0984 90 + +42 0982 6 — ‘10 0'976 
6 r 1007 ‘I — "54 | 0°967 7 — *36 | 0972 
"7 ‘61 1016 2 — +28 0°975 8 _— 0°977 
‘8 “49 1000 3 — -25 | 0-962 ‘9 — +33 | 0-988 
‘9 ‘34 | 0°993 ‘4 — 18 | 0°965 
5 — °23 0'972 12°0 — +18 0994 
70 -+ 0°28 IOI! 6 +004 0°979 
































2. Velocities of granules 


2.1. Correlation between velocity and surface brightness.—An inspection of 
Fig. 2, a graphical representation of the results in Table IV, shows some evidence 
for the 1949 spectra of a correlation between velocity and surface brightness. 
Bivariate frequency tables between velocity and surface brightness were set up, 
and correlation coefficients were computed in the usual manner. The correlation 
coefficient for the spectrum of 1949 September 11-4306 is — 0-262, and that for 
September 11-4319 is —0-294. The scatter diagram and the computed regression 
lines for the latter spectrum are shown in Fig. 3. Using Fisher’s approximately 
normal z—¢ distribution for the correlation coefficient (7), we find the deviation 
of = from (=o (a completely uncorrelated parent population) is respectively 
2-86 and 3-22 times the standard deviation of z— ¢ for the first and second spectra. 
In other words the probability of getting r= — 0-262, or algebraically smaller, in 
the first spectrum from a completely uncorrelated parent population is less than 
I in 460, and the similar probability of getting r = — 0-294 in the second spectrum 
is less thantin1520. It is therefore unlikely that these are chance correlations. 
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On the other hand, the correlation coefficient for the spectrum of 1950 is 
+0:076. In this spectrum z—{, where {=0, is less than the standard deviation 
of the distribution, or the correlation is meaningless. Applying the same analysis 
to the 340 pairs of measures of Richardson and Schwarzschild (15) we find 
r= — 0-304, once again a significant correlation since z—, for =o, is more than 
five times the standard deviation of the distribution. Finally the 348 pairs of 
measures from all three Oxford spectra give r= —0-153, where z—{ for [=o 
is 2°87 times the standard deviation of the distribution, or for all the Oxford 
measures there is a significant negative correlation. 





v—V T T T T T l ! 
+05 
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+0°3 +- 


+01 }— 


—O1 


~03 


— Fs 





—0°7 





! | | ! L le | 








092 094 096 098 1:00 1:02 1:04 1:06 1:08 
kJ 
Fic. 3.— Scatter diagram of V—V on kJ for the spectrum of 
1949 September 11°4319. The two straight lines are the 
regression lines computed by least squares. 


Comparing the correlation coefficient, r = — 0-153, for the 348 pairs of Oxford 
measures with that of r= — 0-304 for the 340 pairs of Mount Wilson measures, 
the smaller numerical value of that for the Oxford measures is undoubtedly a 
consequence of the three times larger error of velocity-measurement (Section 
1.21). The small numerical value of both coefficients further suggests that the 
significant negative correlation which exists is being attenuated by some velocity 
perturbation uncorrelated with surface brightness. The dispersion of this 
velocity is presumably comparable with or larger than that of the velocity associated 
with surface brightness, since in one of the four spectra, that exposed at Oxford 
in 1950, the correlation is completely masked. 

2.2. Velocities of granules.—Table V summarizes the standard deviations of 
velocity and surface brightness and their correlation for the three Oxford spectra 
and for the Mount Wilson spectrum. The last column of the table gives the 
equation of the regression of V—V on kJ, that is, the least-squares solution on the 
assumption that the “errors” are wholly in V—V. In the three spectra where it 
exists the correlation is such that the bright photospheric matter has a negative 
sight-line component of velocity, or since cosy, —1 the bright photospheric 
matter is moving vertically upwards and the dark matter is moving vertically 
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downwards. This is what we should anticipate if the bright photospheric 
matter were composed of convection cells moving vertically upwards in a zone of 
instability. 

In an unpublished paper on the Physics of the Photosphere (14) it has been 
suggested that the photosphere may be divided into two zones. First, below an 
optical depth of unity (for 5485 A), but above Unsdld’s hydrogen-convection zone, 
lies a region of instability with a steep temperature gradient. It possibly arises 
from the onset of boundary-cooling. From the disappearance of granulation 
towards the limb of the Sun it may be shown that the upper boundary of the 
granules also occurs approximately at an optical depth of unity. Secondly, 
above the zone of instability and granulation lies a roughly isothermal region where 
the absorption lines are formed. 








TABLE V 
Sp. | o(V—T) o(kJ) —_ r Regression of V—V on kJ 
I 0°20 0°023 4°°2 —o-262 V— V = —2-359(kJ —1) 
2 0:26 0'026 ie —0°294 V—V=—2-960(kJ —o0-999) 
3 0°29 0'021 3° °6 _t 
4 0°26 0021 17-8 —0°304 V—V=—3°558(Am—o-002) 




















Sp. 1 is Oxford spectrum of 1949 Sept. 11-4306. 
Sp. 2 is Oxford spectrum of 1949 Sept. 11-4319. 
Sp. 3 is Oxford spectrum of 1950 May 11°7277. 
Sp. 4 is Mt. Wilson spectrum of 1949 May 24-6688. 


lf this suggestion be correct, then the Doppler displacement of the observed 
absorption lines, formed predominantly above the granules, will reveal only a 
fraction of the velocity of the granules themselves. This has indeed already been 
pointed out by de Jager and Pecker (5). In order to get a quantitative estimate 
of this effect, assume that the absorption lines are formed by pure absorption 
(Eddington’s «= 1) so that the equation of transfer becomes 


dl 
cos >* =(1+n)J, —(I+7)B,, (4) 


where 7 is the optical depth for the continuous spectrum and 7 is the ratio of 
line to continuous absorption coefficient—all quantities being defined for the 
wave-length A. The surface brightness of the Sun at this wave-length will 
therefore be 


I, (0, 0) = |. B, (r) exp [—sec@ \@ +) dt] (1+) sec@. dr. (5) 


Now let the top of the granulation zone occur at r=79, and let 
= (aconstant) for O<r<7p, 
n=, (aconstant) for ty9<7<0o. 
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Then the surface brightness may be written 


I,(0, 0) = | " Br) exp [ — (1 + y9)7 sec 4] (1 + yo) sec 8 dr 


+ exp[—(I +7 )7) sec 4] | ; B,(r) 
x exp[—(1+7,)(7 — 7) sec O](1 + 7;) sec 0 dr 


ager secO B, y 
” Ce *) expt —yldy 


“0 





'cos 6 
+exp|—(1 + y)7_ sec al| B, (=0+ on expt) dy 
0 I+” 


=B, [+4 = — exp[- (1+ 7)» ee 6] (—— =) peose |, (6) 


where the integration od been effected on the assumption that the source function 
is linear in the optical depth, or B,(r) = By(1 + Ar). 

Hence the dip in the absorption line, R,(A)=1—J,(0, 6)/1)(0,@), where 
I,(0,0) is the surface brightness in the adjacent continuous spectrum so that 
No =, =0, is given by 

age~ hae [ o_ + exp[—(1 + n)rpsec 8) ( —_— )] 


1+Pcos@| 1+ % I+ I+”, 


Taking @=0, corresponding to the centre of the disk, 7»=1 for the boundary 
of the granulation zone, and B=1-9 from the darkening to the limb in the 
continuous spectrum at 5900 A, we have 


Rd)= $3[ 7 +expt-(r+m)l (= - =) ]- (7 


TABLE VI 


























Actual V Obs. V Equiv. width 





km sec™! km sec~* 
° ° 48°5 mA 
I 0'213 504 mA 
2 0°472 54°4 mA 

















Since the average Rowland Number of the solar absorption lines in Table II 
is 3°5, we shall assume that the line absorption coefficient is determined by the 
Doppler core. We may therefore write 


no=Cexpl—A—A}*64, m=Cexpl-A-- =>)", 8) 


where V is the sight-line component of the velocity in the granulation zone. 
For T =6000 deg. K, A= 5910 A and atoms of Fe, b=0-0263A. To give approxi- 
mately the equivalent width of a line of Rowland No. 3°5 we put C= Io. 

From equation (7) we may calculate the profile of the absorption line for any 
required value of V, and from the profile may find the position of the centroid by 
numerical integration. The results of the calculation are given in Table VI. 
We note that even if the sight-line component of the velocity of the granular 





No. 2, 1954 Motions in the Sun at the photospheric level 265 


matter is 2km sec !, the observed line has a Doppler displacement corresponding 
to only 0-47kmsec'. Also possibly worthy of comment is the increase of 
equivalent width with increasing velocity of the granular matter. We thus have 
a mechanism for increasing the length of the Doppler part of the curve of growth 
without the necessity of assuming isotropic turbulence in the zone where the 
absorption lines are predominantly formed. 

From Fig. 2 the most marked fluctuations of surface brightness, due 
presumably to granulation, lie between 0-95 and 1-05 for the spectrum of 1949 
September 11-4319. From the equation of the regression line for this spectrum 
the corresponding observed velocities of descent and ascent are 0-145 and 0-151 
kmsec!, or from Table VI the actual velocities are 0-70 and 0-75 kmsec™! 
respectively. Since the model of absorption-line formation we have used is 
very crude, these “‘actual” velocities can only be regarded as estimates. 

2.3. Granulation and granules.—What now is the relation between the granu- 
lation shown in our spectra and the granules in Keenan’s remarkable photograph 
(6)? In this photograph the smallest granules have a diameter of approximately 
1”-1 of arc, a contrast of 0-15 and a density of one granule per four square seconds 
of arc. On account of poor seeing and imperfect telescopic definition, each 
luminous point on the solar surface may be regarded as being smeared out in the 
focal plane into an image with a flux distribution 


I x?+y 
2nae*P (- 2a" may, 


where a is the radius of the so-called tremor disk. According to measures by 
Wanders (16) on the solar limb as observed at Utrecht, a may lie between 2”-2 
and 5” of arc. The illumination at a point (x, y) on the solar image will therefore 
be given by 





O(%,y)= 3 [[Ménexpf-S— PVN Gedy, 9) 


2 
2G" . 


Assume each Keenan-granule to be a square of 1” side centrally situated in a 
square of 2” side (the sides of the granule and its surrounding square being parallel 
to the coordinate axes), with a contrast of granule to surrounding square of 0-15, 
and that the surface of the Sun is closely packed with such squares and their 
central granules. ‘Then equation (Q) may be integrated in terms of error functions 
with the result for O(x,0) shown in Fig. 4. Even when the seeing is so good that 
the tremor disk has a diameter of one second of arc, a condition rarely realized 
in solar seeing, the observed contrast is reduced to 0-037 from its true value of 
0-15. With a tremor disk one and a half seconds in diameter Keenan’s granules 
will have effectively disappeared, and the Sun will appear to have an almost 
uniformly bright surface. 

Nevertheless, observations, such as the present ones, made with a tremor 
disk larger than that applying to Keenan’s photograph show a granular structure. 
This must be a consequence of the fact that the Keenan-granules are not, as 
assumed in Fig. 4, uniforml? spaced over the surface of the Sun. From the 
figure we may take it that the smallest structure visible on the surface will be 
roughly fixed by the tremor disk, and that consequently the diameter of the 
tremor disk in the Oxford observations is equal to the diameter of the bright 
maxima, say 3"°7 of arc. Consider then an imaginary graticule on the surface 
of the Sun in the form of squares with sides of 3"°7 of arc. The average number 


18 
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of Keenan-granules in a square will be m=(3-7)?/4=3-4, since on the average 
according to Keenan there is one granule per four square seconds of arc. The 
chance of finding m Keenan-granules in a square will be given by a Poisson 
distribution, or 


P(n)= = exp(—m), 


where m= 3°4 is the average number, and \/m=1-84 is the standard deviation 
of the distribution. The squares with an excess over or a defect under m will 
be arranged at random over the surface, and we may take many spectra showing 
little or no fluctuation in surface brightness. However, the spectra we have 
measured (cf. Table V) show a dispersion of surface brightness with an average 


Surface brightness 
1-15 T T T | 





b a=0707 
1:10 


1-05 























1-00 | ! i | ale, 
pa ania” 0” +1” +2’ 





Fic. 4.—TIllumination of spectroscope slit when the tremor disk has a 
diameter, 2a, of © seconds of arc (rectangular lines), and 0707, 1*0 and 
1°414 seconds of arc (curves). 


standard deviation of 0-023. On our proposed interpretation a square with a 
surface brightness of 1-023 arises because of an excess, Am, of granules, where 
0°15(m + Am) + (3-7) 
o-15m-+ (3:7)? 

since each Keenan-granule has an area of one square second of arc and a contrast 
of 0-15. For m=3°4, Am=2-24, or 1-22 times the standard deviation of the 
Poisson distribution. Since the spectra were selected for the contrast of the 
granulation, this excess of Am over the standard deviation is not unreasonable. 
The granular structure shown in the Oxford spectra can at least be understood 

on this hypothesis of chance crowding of Keenan-granules. Unfortunately, 
since the upflow and particularly the downflow in regions containing an average 
number of granules on the one hand and an excess on the other are likely to be 
very different, it does nut seem possible to use the present observations to estimate 
the velocity of upflow in an individual Keenan-granule. All that can be said is 
that the velocity of upflow of 0-75 kmsec~!, deduced in Section 2-2 for a region 





= 1023, 
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of surface brightness 1-05, is likely to be a lower limit to the actual velocity of 
upflow in a Keenan-granule. 


3. Velocities of localized regions 

The smallness of the correlation between velocity and surface brightness 
suggests, as already remarked, that the velocity of granulation is being either 
partially, or in one spectrum completely, masked by some other velocity-disturb- 
ance. ‘To find the nature of this disturbance we may use the observed surface 
brightness of Table IV, and the equation of the regression of V—V on kJ in 
Table V, to calculate the granulation velocity, (V—V),ran, from the observed 
surface brightness. Subtracting this calculated quantity from the observed 
V—V in Table IV we have left the residual velocity, (V — V),.,, due to this other 
source of velocity-disturbance. In order somewhat to reduce the effect of errors 
in the velocity-measurement a moving average, three at a time, is taken of 
(V —V) res, So that we have values of this quantity for 114 instead of the original 
116 heights in each spectrum. ‘The upper two curves of Fig. 5 give (V—V)ru.s 
for the two Oxford spectra of 1949, and the fourth curve gives the similarly 
calculated (_V — V),., for the Mount Wilson spectrum. Since there is no corre- 
lation between surface brightness and velocity in the Oxford spectrum for 1950, 
the curve for this spectrum in Fig. 5 differs from that in Fig. 2 only in that the 
observed points are a moving average, three at a time, of the original velocity 
measures in Table IV. ‘This curve therefore represents a combination of effects 
due to granulation and to the other source of velocity-disturbance, while the other 
three curves are solely due to the latter. 

It is obvious from the appearance of the three curves for (V— V),,, that we 
have oscillatory velocity fields, and not purely random velocities. ‘This is over- 
whelmingly confirmed by the usual statistical tests for random series (8). Thus, 
for example, the number of turning points (peaks and troughs in the curves) differs 
from that predicted on chance, namely (2N—4)/3, where N is the number of 
observed points, by many times the standard deviation of the latter quantity. 
Incidentally it is satisfactory to note how closely the run of the Mount Wilson 
values of (V — V’),., confirms that found at Oxford. The only difference between 
them is the slightly finer structure of the velocity field shown by the higher- 
resolution Mount Wilson spectrum. 

These oscillatory velocity fields, showing the sight-line component of the 
motion, suggest the presence in the photosphere of large-scale, organized motions. 
Unfortunately the available evidence does not indicate how long such motions 
persist, nor at what level in the photosphere they occur. With regard to persis- 
tence it is true that the two Oxford spectra of 1949 were exposed within two 
minutes of each other but show no recognizable features in common. This, 
however, is as likely to be due to imperfect guiding as to short lifetime; the 
observer’s attention was concentrated on exposing at the instant of best definition 
rather than on accuracy of guiding. With regard to location we can only note 
that if these organized motions occur in the zone of instability the actual velocity 
component must from Table VI be of the order of 2kmsec"!. If this be the 
case it is perhaps a little remarkable that they are unaccompanied by any change in 
surface brightness. From this point of view it is perhaps slightly more probable 
that the motions are occurring in the upper photosphere so that the measured 
velocity is the actual velocity component. 
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If the present evidence, and that earlier obtained by Miss Hart (4), are to be 
trusted, there is possibly a hierarchy of organized motions in the photosphere. 
From her solar rotation spectra Miss Hart has measured common horizontal 
components of velocity extending over several degrees of longitude at the equator 
and persisting for at least three hours. Expressed in seconds of arc at the centre 
of the disk the observed extent of these motions lies between 30” and 80” of arc. 
These distances are comparable with the length of the slit in the present obser- 
vations, and suggest that possibly the values of V, in Table II are the vertical 
components of the organized motions of which Miss Hart has measured the 
horizontal component. At the lower end of the scale of size of organized motions 
come the oscillatory fields of Fig. 5 with an average size of some 5” ofarc. Whether 
in this possible hierarchy of organized motions we are for example observing the 
conversion of potential and internal energy into kinetic energy in accordance with 
the theorem by Margules, or observing some small part of a general photospheric 
circulation, can only be settled by the accumulation of much more evidence. 


Conclusion 


It is satisfactory that in a somewhat difficult observational problem there is 
such good agreement between the velocities and surface brightness measured 
at Mount Wilson by Richardson and Schwarzschild and at Oxford by myself. 
Where the two investigations part company is that Richardson and Schwarzschild 
believe the velocities to be purely random and regard them as a manifestation of 
homogeneous, isotropic turbulence, whereas I believe they originate in large-scale, 
organized motions. 

The similarity theory of turbulence tells us nothing about the large-scale, 
turbulence-producing processes. When these processes occur they will no 
doubt be accompanied by small-scale isotropic turbulence, but the astrophysically 
interesting problem lies in the occurrence and the origin of the large-scale processes 
which contain most of the energy. ‘The present investigation suggests that there 
are two categories of turbulence-producing processes. In the zone of instability 
in the lower photosphere convection cells—perhaps identifiable with Keenan- 
granules —move vertically upwards with a velocity of at least three-quarters of 
a kilometre per second. Also in this zone, or rather more probably in the upper, 
approximately isothermal zone occur regions of apparently common motion with 
dimensions ranging from 5” to perhaps 100” of arc. If there is isotropic turbu- 
lence in the photosphere, these may be among processes which produce it. 

In conclusion I wish to thank Mr R. J. Bray for permission to quote his new 
results on the apparatus function prior to publication, and Mr D. W. N. Stibbs 
for correcting an error in the application of Fisher’s z—{ distribution. 


University Observatory, . 
Oxford : 
1954 April 30. 
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PHOTOGRAPHIC AND PHOTOVISUAL MAGNITUDES OF 
7™rjom STARS IN THE +15° SELECTED AREAS* 


A. Beer, R. O. Redman and G. G. Yates 


Summary 


In order to provide a closer link between star magnitudes in the northern 
and southern skies, pg and pv magnitudes of 1268 stars in the twenty-four 
Selected Areas at + 15° have been determined on the international magnitude 
system of the NPS. The first part of the programme was carried out by 
in-focus photography, using the Cape photometric cameras, following approxi- 
mately the same procedure as that used at the Cape. The standard stars 
were taken from the Seares, Ross, Joyner (SRJ) catalogue. Stars in nine 
of the Selected Areas were also measured photoelectrically, using a pulse- 
counting photometer, again with comparison stars near the Pole. In 
addition 135 SRJ stars were measured photoelectrically to verify the relation 
of SRJ magnitudes to the NPS, while all NPS stars brighter than 10™-7 pv 
were measured photoelectrically, to examine their self-consistency of scale. 

After taking into account the scale and colour errors of SRJ, the inter- 
agreement between photographic and photoelectric measures at -+-15° is 
satisfactory. For nine Areas together the mean zero difference (photoelectric 
minus photographic) is -+-0o™-007+0™-o10 pg, and +o0™-019+0"-007 pv. 
Use of a comparatively short series of measures made earlier at the Cape has 
allowed a provisional comparison with the S magnitude system, the present 
basis of stellar photometry at the Cape. 

With single photoelectric readings in the order “‘ standard star S, field 
star F, standard star S, sky”’, the standard error of the field magnitude 
is oO™-06 pg, o™-045 pv. The standard error of one photographic measure of 
an individual star, based on four exposures in the order SFFS, is o™-06 pg, 
o™-065 pv, provided the plate zero error is known independently from photo- 
electric measures. If there is no such information with regard to zero error, 
the standard error for one star measured photographically on one plate only 
is O™-08 pg and o™-08 pv. There is good evidence that in England the errors 
of stellar photometry are determined chiefly by the sky, and depend heavily 
on the stringency with which observing nights are selected. The standard 
error of one catalogue magnitude varies, according to Area, over the ranges 
o™-930-0"-045 pg, o™-035-0™-045 pv. 

The paper concludes with a few comparisons with other catalogues and a 
short examination of the relation between spectral types and colour indices. 


* The full text of this paper is published in Memoirs of the R.A.S., 67, Part I, 1954. 
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